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The Hubbard model on a square lattice is one of the most studied condensed-matter quantum 

("^ ■ problems. Here we find evidence that for intermediate U/4t values and a hole-concentration range 

^Nj ' X € {xc, a;*) the ground state of the Hubbard model on the square lattice perturbed by weak three- 

j, , dimensional (3D) uniaxial anisotropy has long-range d-wave superconducting order. Here t is the 

rS ,1 effective nearest- neighbor transfer integral and U the effective on-site repulsion. The lower critical 

.^^ ' concentration Xc involves the Ginzburg number Gi and is approximately given by Xc ^ Gi + xo ~ 

I 0.05. Here xo < Gi is a small critical hole concentration that marks a sharp quantum phase transition 

00 . from a Mott-Hubbard insulator with long-range antiferromagnetic order for a; < a;o to an Anderson 

insulator with short-range incommensurate spiral spin order for x G (a;o, Xc). The value of the critical 

hole concentration a;* depends on U/4t and is given by a;, « 0.27 for U/At « 1.525. The long-range 

d-wave superconducting order emerges below a critical temperature Tc for a hole concentration 

f^ ' range centered at Xop = {xc + a;«)/2 ~ 0.16. It results from the effects of the residual interactions 

of the charge c and spin-neutral two-spinon si fermions of Ref.i, as a by-product of the short-range 

spin correlations. The spin subsystem provides through such interactions the energy needed for the 

effective pairing coupling between the c fermions of the virtual-electron pair configurations. 
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I. INTRODUCTION 

^ . The Hubbard model on a square lattice^" is one of the most studied condensed-matter quantum problems. However, 

it has no exact solution and many open questions about its properties remain unsolved. Besides being a many-electron 
problem with physical interest in its own right, there is some consensus that it is the simplest toy model for describing 
the effects of electronics correlations in the cuprate superconductors^"^*' and their Mott-Hubbard insulators parent 
vj [ compoundsiiiiii^. In addition, the model can be experimentally realized with unprecedented precision in systems of 
f-y~. • ultra-cold fermionic atoms on an optical lattice. One may expect very detailed experimental results over a wide range 
QQ I of parameters to be availablei^. The pairing of fermions lies at the heart of both superconductivity studied in this 
(T^ . paper and superfluidity. Recent studies of a system of ^Li ultra-cold atoms observed a zero-temperature quantum 
Cn| ' phase transition from a fully paired state to a partially polarized normal state^**. Hence, our studies are of interest 
.,-tI , for both cuprate superconductors and systems of ultra-cold fermionic atoms on an optical lattice. 
f~^ ■ The virtual-electron pairing mechanism studied in this paper is consistent with the evidence provided in RefsJ^ii^ 

00 , that unconventional superconductivity is in different classes of systems such as cuprate superconductors, heavy-fermion 
^D ■ superconductors, and iron arsenides mediated by magnetic fluctuations. Our investigations have as starting point the 
^ ' square-lattice quantum liquid of Refi, which refers to the Hubbard model in the one- and two-electron subspace-. It is 
spanned by the ground state and the excited states generated by application onto it of one- and two-electron operators. 
For such a square-lattice quantum liquid only the charge c fermions and spin-neutral two-spinon si fermions play an 
active role. The c and si fermion description and the related general rotated-electron description are consistent with 
the global 5*0(3) x SO{3) x U{1) = [80(4) x U{l)]/Z2 symmetry found recently in RefJI for the Hubbard model on 
any bipartite lattice and thus on the square lattice. Such a global symmetry is an extension of the 50(4) symmetry 
known to occur for the model on such lattices ^^. 

An extended presentation of the results of this manuscript, including their application to the study of the unusual 
properties of the hole-doped cuprate superconductors, can be found in Ref.™. 

The Hubbard model on a square lattice with torus periodic boundary conditions, spacing a, N^ = [-^a]^ sites, 
lattice edge length L = Na a, and Na ^ 1 even reads, 

H=-t Y. E [4„. cry,'^ + h.c] + U [N! ~ Q]/2 ; Q = ^ E '^-^..- (^ ^ "-~.— ) • (1) 

Here np.^a- = cp. cr'^^-.a and — cr =^ (and —a =1) for a =1 (and a ='[) counts the number of electron singly occupied 
sites. Consistently, the operator D = [N — Q]/2 where N = ^^ N^, and A^^r — X])=i ^^rj.o- counts that of electron 



doubly occupied sites. We denote the 77-spin (and spin) value of the energy eigenstates by S*,, (and Sg) and the 
corresponding projection by S^ — ^Wa ^ ^]/2 (and S^ = —[N^ — iV^]/2). We focus our attention on initial ground 
states with hole concentration x — [N^ — N]/N^ > and spin density m — [N^ — N^/N"^ — and their excited states 
belonging to the one- and two-electron subspace. 

As found in Ref."'^, one can perform an extended Jordan- Wigner transformation^^'^^ that maps the si bond particles 
onto si fermions with operators /I. ^-^ and fp ,si whose real-space coordinates r, where j — 1,...,N^ are those of 
the square si effective lattice with N^_.-^ w (1 — x)N^/2 sites. Such operators are related to the corresponding si 
bond-particle operators 5I ^^ and gf .si, respectively, defined in Ref3^ as follows, 

si = X! 4 , sifr,,,si (l>j',j,si ; (l)j',j,si = arctan ( -^ ^ _ -^^ ) . (2) 



I^J,SL 



Here the phase (j)j'j.si is defined in the range 4>j\j,si G (0, 27r) and Xj-^ and Xj^ (and Xji-^ and Xj'^) are the Cartesian 
components of the vector r*,- (and fj'). The phase (fij^si is that created by a gauge field whose fictitious magnetic field 
and corresponding effective vector potential read. 



Bsi{rj) = Vf^ X Asi{fj) = *o ^ nf^,.si S{fj' - tj) e. 



3:3 ' 
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respectively. In these expressions 62:3 is a unit vector perpendicular to the plane, $0 = 1 in our units, and np.gi = 
/f si fr .si is the si fcrnfion local density operator. In equation ([3]) and the remaining of this paper we denote by. 
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(4) 



an in-plane unit vector whose direction is defined by the angle (j). The effective potential ([3]) generates long-range 
interactions between the si fermions. In addition, each spin-neutral two-spinon si fermion has in average a fiux tube 
of one flux quantum attached to it. 

The components qj and qj of the discrete momenta qj of the si fermions are the eigenvalues of the two 

corresponding si translation generators in the presence of the fictitious magnetic field Bsi of Eq. (j3|). The si fermion 
operators are defined in and act onto subspaces spanned by mutually neutral states^. Therefore, such operators 
commute and the two components qj^ and qj^ can be simultaneously specified. A property that plays a central 
role in the studies of this paper and Refi^ is that for vanishing spin density m — ground states at finite hole 
concentrations a; > the si fermion momentum band is full and for one-electron and two-electron excited states 
displays a single hole and none or two holes, respectively. The si - si fermion interactions associated with the 
effective vector potential of Eq. ^ are stronger than those that arise between the emerging si fermions and pre- 
existing c fermions. In spite of that, the former do not lead to si - si fermion inelastic scattering. The obvious reason 
is that due to phase-space restrictions associated with the exclusion principle and energy and momentum conservation 
requirements there are no available momentum values in the si band for excited-state occupancy configurations. 

According to Mermin-Wagner-Berezinskii Theorem3i"— , in 2D quantum systems destruction of the fluctuations of 
long-range orders occurs at finite temperatures. To search for long-range superconducting orders at finite temperature 
we add a small three-dimensional (3D) uniaxial anisotropy perturbation to the square-lattice Hamiltonian ([l}, 

HsD^H + H^; H^ = -t^ Y. E [4 .- ^n' ,a + h.c.]; M=^. (5) 

(rjr,.,;_L)(T=T,; ^ 

Here the sum '^/fr.f^.,.j_\ runs over first-neighboring sites on nearest-neighboring planes, ij^ ^ t is a small transfer 
integral associated with electron hopping between such planes, and the mass Af is given in units of lattice spacing 
and Planck constant one. For the Hamiltonian iJs^i of Eq. ([5]) N denotes the expectation value of the number of 
electrons per plane. In addition to weak 3D uniaxial anisotropy, in this paper we find evidence that for the hole 
concentration range x g (xc,x^,) the effects of the cuprates intrinsic disorder are very weak and can be accounted for 
by an oversimplified scheme involving a single effective suppression coefficient 7^. Here Xc ~ 0.05 is a critical hole 



concentration introduced below and x* is the critical hole concentration above which there is no short-range spin 
order at zero temperaturei. 

In this paper we consider three completely different types of anisotropy: (i) The 3D uniaxial anisotropy associated 
with the small Hamiltonian term H± of Eq. ([SJ ; (ii) The in-plane anisotropy as measured by the values of the Fermi- 
velocity anisotropy coefficient ija — max V^-^/Vpc and Fermi-energy anisotropy coefficient ryo = jAj/W^'' introduced 
in Refji and given in Section II - It refers to the anisotropy of the Fermi line associated with the dependence of its 
finite energy width on the Fermi angle (j) given in Eq. (jA6p of Appendix A; (iii) Finally, the in-plane superfluid-density 
anisotropy, which refers to the in-plane anisotropy of the superfluid density of some cuprate superconductors. These 
three types of anisotropy are in the following in general called 3D uniaxial anisotropy, anisotropy, and superfluid- 
density anisotropy, respectively. 

The effects of the 3D uniaxial anisotropy perturbation H± on the square-lattice Hamiltonian ([T|) can for very 
small anisotropy parameter e^ = m*/M <^ 1 be ignored for most properties. Here m* is the c fermion mass of 
expression (IA11[) of Appendix A and M = l/2t± ^ l/2i the effective mass given in Eq. ([SJ associated with electron 
hopping between nearest-neighboring planes. The c fermion mass becomes infinite in the limit U/At -^ 0. Hence in 
spite of t±/t <C 1 for the original Hamiltonian H^d of Eq. ([S]), which is written in terms of electron operators, for 
small U/4:t <C uq « 1.3 the 3D uniaxial anisotropy parameter reaches values e^ — m*^/M > 1 larger than one upon 
decreasing C//4t. Indeed, rather than the bare electronic transfer integral t and corresponding mass m^ = l/2t, the 
relevant parameter for that quantum liquid is the c fermion mass m* — l/(2rct). That follows from the occupancy 
configurations of the c fermions generating the charge degrees of freedom of states close to the energy eigenstates of 
H3U. The parameter Tc = m'^ /m* in the c fermion mass expression is the mass ratio given in Eq. (jAip of Appendix 
A. According to that equation, its limiting values are Tc — )■ and Tc — )■ 1 for U/4t —>■ and U/4t — >■ 00, respectively. 
Our study refers to the range U/At > uq for which approximately Vc G (0.74, 1.00). For it that t±/t <C 1 for the 
original Hamiltonian H^u of Eq. ([5]) implies as well that £^ — m'^/M <C 1. Here uq « 0.13 is the U/At value at which 
the important energy scale Aq of Eq. (|A4[) of Appendix A reaches its maximum magnitude. 

The studies of this paper focus on the Hamiltonian H^d of Eq. (O in the one- and two-electron subspacel. For 
hole concentrations x G {xc,x^,) we call virtual- electron pair quantum liquid (VEP quantum liquid) that quantum 
problem under the weak suppression effects considered below. The main effect of the small perturbation H± of 
the VEP quantum liquid is for approximately U/At > uq the emergence of the 3D uniaxial anisotropy parameter 
e^ = m* /M <C 1 in the expressions of some physical quantities sensitive to the thermal and quantum fluctuationsi^. 
Specific examples are the critical hole concentration Xc ~ Gi-l-xo, the phases 9 of the virtual-electron pairs considered 
below, and related quantities. Here Gi is the Ginzburg number— and xq < Gi is a small critical hole concentration 
that marks a sharp quantum phase transition from a Mott-Hubbard insulator with long-range antiferromagnetic order 
for a; < xo to a Anderson insulator with short-range incommensurate spiral spin order for x S {xq,Xc)- As discussed 
in Appendix B, the strong intrinsic-disorder hole-trapping effects present in the hole-doped cuprate superconductors 
for x € (0, Xc) render xq a critical hole concentration. Moreover, such effects, which are active mainly for hole 
concentrations below Xc, shift the critical hole concentration from Xc ~ Gi to Xc ~ Gi + xq- This is their only effect 
on the physics for x e {xc,x^). 

Concerning the short-range spin order occurring for x € (xq, x*), for both x e [xq, Xc) and < (a; — Xc) ^ 1 it is a 
incommensurate spiral spin order as that of the square-lattice quantum liquid of Refi for small x. Further details on 
that order for instance for x near the optimal hole concentration Xop — {xc + x^,)/2 and f//4i £ {uo,ui) remains an 
interesting open problem. We find in this paper that for both < (a; — Xc) ^ 1 and < {x.^, — a;) ^ 1 the quantum 
fluctuations are large, so that the VEP quantum-liquid ground state refers to a liquid. In turn, it is found that at and 
near x sa Xop such fluctuations are smaller. Hence the ground state is expected to be intermediate between a liquid, 
where such fluctuations are large, and a crystal, where they are small. The results of Ref.^'^ seem to indicate that the 
VEP quantum-liquid physics is for the latter range of x values controlled by a quantum critical point. 

The possibility of using U/At as a tuning parameter plays a central role in our scheme. Indeed, the change in the 
U/At ratio is strongly linked to a change to the electron-lattice coupling, since increasing U/At also involves increasing 
the strength of the periodic potential. In contrast to large U/At values for which the electron-lattice coupling strength 
further increases and the important energy scale Aq of Eq. (jA4p of Appendix A becomes small, for intermediate 
U/At G (mq, Mi) values that strength remains much smaller, hence freeing the correlated VEP quantum liquid from its 
rigid link to the underlying lattice. Intermediate U/At values then refer to relatively high effective correlation strength 
without an increase in the effective strength of the periodic potential. The magnitude of the energy scale Ag provides 
a measure of such an correlation strength effectiveness, vanishing both for U/At — )■ and U/At —>■ 00 and being largest 
for U/At S {uq,ui). Such a correlation strength effectiveness is probably behind the exotic physics emerging for both 
intermediate U/At S (uq, ui) values and hole concentrations at and near Xop — {xc + x^)/2, which as mentioned above 
is according to the studies of Ref.^"^ controlled by a quantum critical point. 

Hence often our results focus on the intermediate- t//4i range U/At € (uq,ui) and hole concentrations x € {xc,x^,). 



As discussed in Ref.^-, the value of the parameter e^ — rnl/M ^ 1 is set so that Xc ~ G\ + Xq k, 0.05. In turn, 
that of the zero-temperature upper hole-concentration a;* is insensitive to e^ and rather depends on U /At. Alike 
for the square-lattice quantum liquid of Refj^, it changes smoothly from x» « 0.23 for uo « 1.3 to x* « 0.28 for 
Mo ~ 1.6. Specifically, the value cc* = 0.27 appropriate to the families of hole-doped cuprates considered in this paper 
corresponds to U /At « u* = 1.525. Provided that e^ = m*^/M <C 1, it is expected that the momentum occupancy 
configurations of the c and si fermions that generate the energy eigenstates of the e^ — square-lattice quantum 
liquid of Ref.i are for x S (cCcX*) close to energy eigenstates of the VEP quantum liquid. It then follows that the 
interactions of such objects remain being residual. 

The phases of the phase-coherent virtual-electron pair configurations studied in later sections have the form 9 = 
Oq + Oi. Their fiuctuations play an important role in the VEP quantum liquid physics. At zero temperature the 
fluctuations of the phases Oq and 9i become large for x ^ Xc and a; — >■ x* , respectively. Expression in terms of the 
rotated-electron operators^ of the effective microscopic Hamiltonian describing the quantum fluctuations of the VEP 
quantum liquid leads to a problem with basic similarities to that considered in Ref.^. Fortunately, the very involved 
problem of the derivation of an effective action for the phases 9 can be mapped onto a corresponding problem already 
investigated in that reference. That is achieved on replacing in the Hamiltonian (1) of Ref.~ electron operators by 
rotated electron operators^. The transformation that relates rotated electrons to electrons is unitary, so that the 
effective action for the phases 9 considered in this paper is valid for approximately U/At > uq and specifically at 
U/At ^ u^, = 1.525. Expression of the quantum problem that describes the fiuctuations of the phases 9 — 9o + 9i of 
the virtual-electron pairs in terms of c and si fermion operators reveals the central role plaid by the c - si fermion 
interactions in the emergence of a long-range superconducting order for x G {xc,x^,). 

Indeed, strong evidence is found below that for intermediate U/At values, vanishing spin density to = 0, and 
hole concentrations in the range x € {xc,Xi,) the ground state of the VEP quantum liquid has a d-wave long-range 
superconducting order, coexisting with its short-range spin order. Superconductivity emerges below a x dependent 
critical temperature Tc and results from the effects of the residual c - si fermion interactions, as a by-product of 
the short-range spin correlations. The spin subsystem provides through such interactions the energy needed for the 
effective pairing coupling between the c fermions of the virtual-electron pair configurations introduced in the following. 
The suppression effects due to intrinsic disorder or superfluid density anisotropy slightly lessen the magnitude of Tc 
and related physical quantities. Otherwise the jd = ^ physics remains unchanged under these small effects. 

Our scheme is used in Refs.»iS'— in the further understanding of the mechanisms behind the unusual properties of 
hole-doped cuprates. The studies of these references focus on five representative systems of the hole-doped cuprates 
for which Xc ~ 0.05 and x^ « 0.27: YBazCugOg+a (YBCO 123), BisSrsCaCusOg+a (Bi 2212), HgBa2Cu04+5 (Hg 
1201), Tl2Ba2Cu06+5 (Tl 2201), and La2-2:Sr^Cu04 (LSCO). In Sertion III-D it is found that for jd ~ I some of the 
VEP quantum liquid expressions are valid provided that Xc/x^ < 1/4. The critical hole concentration values Xc ~ 0.05 
and x^, « 0.27 are compatible with such a requirement. One finds {xc — xq) « Gi « 0.026 for the four representative 
hole-doped cuprates other than LSCO and {xc - xo) w Gi w 0.037 for LSCO. The relation Gi = G/e^ = G [M/to*] 
then implies that our scheme is not valid for the e^ = ml/M = 2D square-lattice quantum liquid studied in Refi 
Indeed, such a problem corresponds to vanishing values of e^ = to*/M, such that 1/e^ = M/to* -^ cx) and thus 
Xc/Xi, > 1/4. Fortunately, one finds that G is for the five representative systems very small, G sa 10^^ — 10~^. This 
allows that the ratio e^ — ml/M « 10~^ — 10^^ is smaller than Gi. Interestingly, our results provide evidence that 
the zero-temperature hole-concentration width (a;* — Xc) of the superconducting dome decreases upon decreasing the 
3D uniaxial anisotropy. This is consistent with the weak 3D uniaxial anisotropy effects associated with the small 
Hamiltonian term H± of Eq. ([5|) being needed for the occurrence of superconductivity in the quantum liquid under 
investigation in this paper. 

Our preliminary studies on the five representative systems focus mainly onto the hole concentration range x e 
{xc,x^.) where Xc ~ 0.05 and x* sa 0.27 at U/At ~ u* = 1.525 for which the strong hole trapped effects caused 
by intrinsic disorder discussed in Appendix B are not active. For that x range the suppression effects originated by 
intrinsic disorder or in-plane superfluid density anisotropy are very weak for the four representative hole-doped cuprates 
other than the random alloy LSCO. Indeed, the minimum magnitude 7™™ reached at Xop = (xc-|-x*)/2 « 0.16 of the 
corresponding suppression coefficient 7^ = Tc/Tc\aa=o is found in Reff^^ to belong to the range 'y™" g (0.94,0.98), 
whereas 7^ = 1 both for < (x — Xc) <^ 1 and < (x* — x) <C 1. Our investigations and those of Refi^ provide 
evidence that for x S (xc, x*) the interplay of the electronic correlations described in this paper by the c - si fermion 
residual interactions with the weak effects of the 3D uniaxial anisotropy perturbation H± of Eq. ([5]) is behind the 
unusual universal properties of the hole-doped cuprate superconductors^Ti^. 

For the VEP quantum liquid referring to the range x S (xc,x*) the x dependence of the fluctuations of the 
phases 6* = ^o + ^1 is for intermediate U/At values found to imply naturally a dome-like dependence on x for the 
critical temperature Tc- It reaches its maximum magnitude J"™"^ at the above-mentioned optimal hole concentration 
Xop = (xc + x^)/2 « 0.16 for Xc — 0.05 and x^ ~ 0.27. The energy parameter Aq appearing in the critical temperature 
and other energy scales expressions plays an important role in the square-lattice quantum liquid^. It vanishes both in 



the limits U /At — >■ and U /At -^ oo and goes through a maximum magnitude max{Ao} ~ t/ir at U /At = uq ~ 1.302. 
For instance, it is found in RefJ^ to control the dependence of T™°^ on pressure P. Under application of P it 
first increases until it reaches a maximum magnitude at a pressure value P = Pq. A generalized formula valid for 
arbitrary magnitudes of the critical hole concentrations Xc and x* and pressure P e (0, Pq) is introduced. At P = it 
becomes the empirical formula (1 — Tc/T/!^°'^) = 82.64 {x — 0.16)^ found by many authors to apply to several families 
of hole-doped cuprates for the range x g {xc,Xs/^^. Importantly, an universal ratio r™"^(Po)/r™°^(0) ~ 1.26 is 
predicted for hole-doped cuprates with P — Q and T = critical hole concentrations Xc ~ 0.05 and x* « 0.27. The 
validity of such a theoretical prediction is confirmed for two systems with completely different values of r'™°^(o) 30i3i . 

Moreover, the low-temperature incommensurate peaks in the inelastic neutron scattering of LSCO observed at 
momenta 5P = [tt ± 27ra;,7r] and 5P — [tTjTt ± 2'kx\ for the hole concentration range x G (0.05, 0.12)iii^"— follow 
within the VEP quantum liquid physics from a suitable generalization of the si fermion microscopic processes found 
in Ref.^ to lead at a: = to the sharpest features in the neutron scattering of LCO. The incommensurate character 
of the LSCO peaks results from the contraction of the si momentum band boundary line upon increasing the hole 
concentration x. 

The VEP quantum liquid is consistent with the coexisting two-gap scenario2i^i22i^: A dome-like superconducting 
energy scale 2 1 fi I = 4fcBTc/(l— a;cT'c/2;* T™"^"^) and pseudogap 2|A| = (l—x/x*) 2 Aq, over the whole dome x € [xc^x^,). 
The energy parameters 2|A| and |A| are the maximum magnitudes of the spinon pairing energy 2|Asi((fJj,]^)| = 
2|A|| cos2i?!)| and corresponding one-electron gap |A|| cos2(/)|, respectively. Here q^^-^ is a si band boundary hue 
momentum and d = ±1 the doublicityi. It is reached at the values = 0, 7r/2 of the Fermi angle </> of Eq. (|A6|) of 
Appendix A. The energy scale 2\^\ is the maximum magnitude reached ai 4> = 7r/4 of the superconducting virtual- 
electron pairing energy 2|r2si((j*Jgj^)| = 2|H|| sin2(^| introduced in this paper. Both the magnitudes and x dependences 
of the energy scales 2|A| and 2\Q,\ and of the superfluid density predicted by the VEP quantum liquid scheme agree 
with the corresponding experimental results of the five representative cuprate superconductors. That for x sa 0.04 < Xc 
only the d-wave gap ~ | A| | cos 20| is observed in the experiments of Ref."^^ is not inconsistent with our two-gap picture. 
Indeed, 2\U,\ emerges for x> Xc- 

The paper is organized as follows. The energy scales of the quantum liquid and other basic physical quantities are 
discussed in Section II. In Section III evidence is found that the ground state of the VEP quantum liquid has for 
intermediate U /At values and hole concentrations in the range x g [xc, a;*) a d-wave long-range superconducting order 
coexisting with its short-range spin order. A corresponding effective superconductivity theory is then introduced. The 
results presented in Section IV include the study of the VEP quantum liquid general energy functional, dependence 
on the hole concentration x of the superfluid density, and role of the c - si fermion residual interactions in the c 
fermion strong effective coupling. Finally, Section V contains the concluding remarks. 

II. ENERGY SCALES, CROSSOVER HOLE CONCENTRATIONS, AND FERMI MOMENTUM 

In this section and in Appendix A some of the results of Refs.-i'22, needed for our studies are shortly reviewed. 
That as considered in ReL-l the angle between the c Fermi hole momentum and si boundary line momentum of the 
c fermion and si fermion hole, respectively, created upon an one-electron addition excitation is exactly 7r/2 is a good 
approximation for the study of some properties. However, having in view the study of the one-electron scattering 
properties fulfilled in Refi^, in this paper we account for the small deviations of that angle from 7r/2. 

Concerning our short review of the results of Refsiii^, there occurs for the square-lattice quantum liquid a sharp 
quantum phase transition such that the m = ground state has a long-range antiferromagnetic order at a; = and a 
short-range spiral-incommensurate spin order with strong antiferromagnetic correlations for < a; ^ 1. As a result 
of it, the maximum magnitude 2| A| of the si fermion spinon-pairing energy has a singular behavior at a; = 0. Indeed, 
one has that 2A|a;=o = M° for a; = and m = whereas lim2,_j.o 2|A| = 2Ao. Here ^^ is the chemical potential /z in 
the limit /i" = lima;-i.o /i whose behaviors are given in Eq. (IA3|) of Appendix A. It is one half the Mott-Hubbard gap, 
consistently with for a; = and m = the chemical potential belonging to the range fj, € (— /i°,/i°). In spite of 2/x° 
being the charge Mott-Hubbard gap, it also refers to the spin degrees of freedom^: The energy scale 2/S\x=o = ijP is the 
excitation energy below which the long-range antiferromagnetic order survives for x — Q^ m — 0, and zero temperature 
T = 0. In turn, the energy scale limx^o 2| A| ~ 2Ao plays a key role in the x > and a: > a;o physics of the square- 
lattice quantum liquid of Ref.^ and the VEP quantum liquid, respectively. The small hole concentration Xq < Xc is 
given below. According to the analysis of Appendix B, the hole-traping effects caused for hole concentrations below 
Xc by the cuprates intrinsic disorder are strong and shift that sharp quantum phase transition from a; = to a; = xq- 

The approximate limiting behaviors of the energy scale 2Ao are provided in Eq. (|A4p of Appendix A. As a function 
of u = U/At, it is such that 9Ao(u)/9m = at m = U/At = uq. This is consistent with it reaching a maximum 
magnitude t/n at that U/At value. It can be expressed in terms of the si fermion energy dispersion bandwidth W^i 



as, 

Ao = r, 4W^0i = 4M^0i 6"^= ; W",^ = lim W^i = W^.i U=o , 

a:— J-0 

Wsi = [6,i(gtl.f)-e.i(0)]; A.-|ln(Ao/W°i)|. (6) 

Here eii(<f) is the si fermion energy dispersion derived for the square-lattice quantum liquid in Ref.^ whose generalized 
expression is for the VEP quantum Uquid introduced below and (f^gf is the si boundary line nodal momentum defined 
in that reference. 

For the Hubbard model on the square lattice at zero temperature T = a small concentration x of holes prevents 
the occurrence of long-range antiferromagnetic order—. At vanishing spin density m — G and both for small finite 
temperature T > and zero hole concentration a: = and for vanishing or small finite temperature T > and small 
finite hole concentration < a; <C 1 the system is driven into a renormalized classical regime where the T = and 
a; = long-range antiferromagnetic order is replaced by a quasi-long-range spiral-incommensurate spin order as that 
studied in Refi^ for simpler spin systems. For hole concentrations obeying the inequality < a; < x* the system is in 
a short-range spin ordered phase whose order parameter is the maximum magnitude of the si fermion spinon-pairing 
energy 2|A| ~ 2Ao (1 — x/a*). Following the analysis of Appendix B, for the quantum problem considered here such 
an order occurs for a smaller range x £ {xo,x^,). In either case the above 2|A| expression is valid approximately for 
the range U/4t e {uqjU-^) where u^ > ui is the U/At value at which x* = I/tt sa 0.32. For approximately U/At > ui 
and thus including at U/At sa u^ one has that r^ < re < irs and the approximate expression r^ w g-'^^'^o/u qJ gq_ 
(IA1[) of Appendix A is not valid. Often our results refer to the smaller range f//4i e {uq,ui) for which the relation 
Tc ~ 2rs ~ 2e~'^*"°^^ approximately holds. 

Alike the quasiparticle mass ratios in Fermi-liquid theory^, the [//4i-dependent spin ratio r^ = Aq/AW^i and charge 
mass ratio Tc — m'^ /ml, both given in Eq. (|Aip of Appendix A, control the effects of the electronic correlations in 
many physical quantities. These ratios play an important role in the physics of both the square-lattice quantum liquid 
of Refii and VEP quantum liquid. In contrast to a Fermi liquid, here the mass m^ = limt//4i_j.oo "t* — l/2t, which 
refers to the limit of infinite on-site interaction, plays the role of bare mass. In turn, the VEP quantum liquid mass 
ratio 1/e^ = M/m*^ involving the mass M of Eq. ([5]) and the c fermion mass m* of expression (|A11[) of Appendix A 
controls the effects of the weak 3D uniaxial anisotropy. 

Unlike the Fermi-liquid quasiparticles, the c and si fermions do not evolve into electrons upon adiabatically turning 
off the interaction U . Instead, upon turning off adiabatically the parameter At^ /U the c fermions evolve into the 
spinless fermions that describe the charge degrees of freedom of the electrons that singly occupy sites within the 
energy-eigenstate configurations of the state basis used in the studies of Refi^. Furthermore, for At/U — )■ their 
spin degrees of freedom are described by the spin- neutral two-spinon si fermions whose energy-dispersion bandwidth 
vanishes in that limit. Within a mean- field approximation for the fictitious magnetic field Bsi{fj) of Eq. ^ the 
si fermion occupancy configurations that generate the to = ground states are in that limit those of a full lowest 
Landau level with iV^i = N^^_^ = N/2 one-sl-fermion degenerate states of the 2D quantum Hall effect^. Here N^^_^ is 
the number of both sites of the si effective square lattice and si band discrete momentum values. 

In turn, for At'^/U > the si fermion dispersion acquires a finite energy bandwidth. Both the square-lattice 
quantum liquid of Refii and the VEP quantum liquid are defined in the one- and two-electron subspace. For finite 
U/At values the c and si fermions describe the charge and spin degrees of freedom, respectively, of the occupancy 
configurations of the rotated electrons that generate the states that span that subspace. In it only the c and si 
fermions play an active role. The c momentum band has the same shape and area {2'!r/L)'^ N^ as the first Brillouin 
zone. The si momentum band is exotic and such that its momentum area and shape are subspace dependent. For 
ground states corresponding to hole concentrations x € (XcX*) and spin density m = the c fermions have for the 
square-lattice quantum liquid of Ref.^ an isotropic hole like c Fermi line. The hole momenta qp^ belonging to that 
line, defined in Eq. (jlOp . are centered at the momentum —n. Out of the c band N^ discrete hole momentum values 
Qj^ where j — 1, ..., N^, Nc = 2Sc are filled and N^! = [N^ — 2Sc] unfilled and correspond to c fermions and c fermion 
holes, respectively. For the ground state one has that Nc — 2Sc — N = (1 — x) N^ and the hole c Fermi line encloses 
a momentum area [A^^ _ ^yj [2n/L]'^ = [N^ - 2Sc] [2TT / Lf = At:'^x. 

For ground states corresponding to vanishing spin density m = Q the si band is full. Hence the number of holes in 
the si band vanishes, N^^ — 0, and the number of si fermions reads iV^i — N/2 = (1 — a) N^/2. The si momentum 
band is particle like so that its boundary line is centered at g = and encloses a momentum area 27r'^(l — x). The 
momenta ^Jsi "-"^ ^^^ ^^ band boundary line are given in the following. For the states that span the one- and two- 
electron subspace the number of holes in the si momentum band is given by Nl/'i = 1 and N^^ = or N!,\ — 2, 
respectively^. 

For each [A^-^, N^] — [N/2, N/2 — l]-electron ground state there is a corresponding zero-spin-density and [N^, N^] — 
[N/2, 7V/2]-electron ground state with one more spin-down electron. We say that such a ground state is the m — 
generating ground state of both the [N/2, N/2 — l]-electron ground state and its [N/2, N/2 — l]-clectron excited states 



of small momentum and low energy. The reduced one-electron subspace is spanned by a [N /2, A^/2— l]-electron ground 
state and its [N/2, N/2 — l]-electron excited states of small momentum and low energy. The doublicity d = ±1 labels 
the si fermions whose momentum occupancy configurations generate the states that span that subspace. The si band 
momenta (f of the reduced one-electron subspace si band are related to the si band momenta go of the corresponding 
TO = generating ground state as^, 

q = Ai^q^- limAf, =1. (7) 

a:— vO 

Here we neglect terms of order 1/N^, which vanish in the thermodynamic limit. I denotes the 2x2 unit matrix and 
for X G {Xc,x^) and momenta q at or near the si boundary line one has that A'^i = Ap. The F rotation matrix Ap 
is given below. 

The Fermi-velocity anisotropy coefficient ?7a = raax V^gj^/Vpc, Fermi-energy anisotropy coefficient 770 — \A\/W^, 
and hole concentrations x^ and Xq at which 1]^ — 1 and tjq = 1, respectively, play an important role. For approximately 
C//4i S (mojWi) such quantities read^, 



_ Ao f X, \ _ Ao 1 _ 1 

^° " t {i2r,r + Ao/tJ^ t (2^)2x, ' ^'^i-g- ^' 

Here Ao is the energy parameter of Eq. ((B]), ^c ~ 2rs « 2e~^*'"°/'^ where Uq ~ 1.302, and Wj} ~ x2Tr/m* is the 
unfilled c fermion sea energy bandwidth. Following the discussions of Appendix B, the hole trapping effects caused by 
strong intrinsic disorder render xo a critical hole concentration. It marks a sharp transition between zero-temperature 
states with long-range antiferromagnetic order and short-range spin order. In Ref.— Xci was defined as the hole 
concentration at which the equality 77A — 2xo is satisfied. For f//4i g (wo,wi) that formula leads to Xci ~ 1/8. Since 
Xci is a mere crossover hole concentration in this paper we use the magnitude Xci — 1/8, valid for the interaction 
range [//4i e (uq, ui). The absolute value V^^^ of the si fermion velocity Vg^{q) of Eq. (JAISP of Appendix A at the 
si boundary line and that of the c fermion velocity at hole momenta belonging to the c Fermi line also appearing in 
the expressions of Eq. ([8]) read, 

T^i^.i-Kt('fJ.i) = ^|Bin20|; ^.c = |f-^, (9) 

respectively. The dependence on the Fermi angle (j) G (0, 2tt) defined in Eq. (jA6[) of Appendix A of the absolute value 
^Bsi ^^ ^^^ *1 fermion velocity confirms the anisotropic character of the si boundary line. In turn, that for hole 
concentrations x € (xc, a;*) the c Fermi line is isotropic follows from the independence of the angle (f> of the c fermion 
velocity. 

For X € {xc,x^,) and U/At values approximately in the range U/At € {uq,ui) the Fermi-line level of anisotropy 
is controlled by the interplay of the si boundary line anisotropy and c Fermi line isotropy. The Fermi energy has 
the form Ep — n + SEp where /i w /i° -I- W^, pP is given in Eq. (JB1[) of Appendix B, and below it is found that 
the square-lattice quantum liquid expression 5Ep w |A||cos2(/)| acquires an extra term due to the superconducting 
fluctuations of the VEP quantum liquid studied in this paper. 

According to the criteria of Refji, for the approximate range U /At £ {uq, ui) of intermediate U/At values the Fermi 
line is strongly anisotropic for hole concentrations below Xci = 1/8, has some Fermi-energy anisotropy yet the Fermi 
velocity is nearly isotropic for the hole concentration range x S (a;ci,Xc2), and is nearly isotropic for the x range 
X G (a;c2,a;*). Here the hole concentration Xc2 ~ ([270 + l]/[37o + 1])2;« where 70 = (1 — Xo/xt) is introduced below 
in Section III-E. The Fermi line is hole and particle like for x < Xh and x > Xh, respectively. The value of the hole 
concentration Xh > x^i is not accurately knowni. Likely it belongs to the range Xh G (a;c2, a^*)- Provided that x^ < x*, 
the angle ^an of Eq. (|A6p of Appendix A vanishes for hole concentrations x < Xh and is small for x G {xh,x^,). 
Hence for hole concentrations in the range x S (xh,x^,) we consider only corrections up to first order in (J)an, so that 
miix{5Ep} w |A| and minjy^^J w [^aat V2]\A\. 

The hole Fermi momentum kp = kp + if given in Eq. (IA6|) of Appendix A where 7? = [tt, it] can be expressed in 
terms of the corresponding c hole Fermi momentum qp^ and si boundary-line momentum qj^^i as follows^, 

fc^ = 4c - qisi = 4 W 60 ; qpc = 4c + ^ , 

T*hd ^h I ±d\ r* , T^d 



(iFc 



= q%A(t>i) e^^ ; qLi = IBsM) V, ' d = ±l. (10) 
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These general expressions are for one-electron excited states valid for x G {xc,x^,). For x G {xci,Xc2) the specific 
general dependences on the Fermi angle cj) of the angles (j)f S (0, 2tt) and 0^]^ S (0, 27r) are, 



dTTcsfi - 



4] 



(/.fl = [0,1 + 0' 



4] 



±1. 



(11) 



In the studies of Ref.^ it is considered that for the Fermi velocity isotropic x range x € {xci,Xc2) the c Fermi hole 
momentum qj}^ and si boundary line momentum Qg^i of the c fermion and si fermion hole, respectively, created upon 
an one-electron addition excitation are perpendicular, so that qj}^ ■ q^si — 0. Since ttcs/'^ is the angle between the 
momentum-space directions of qp^r and q^^i, this is equivalent to considering that tTcs/'^ = ''■/2 in the 0J? expression of 
Eq. (in]) . Hence, in spite of the c momentum band remaining unaltered upon such an excitation, for states belonging 
to the reduced one-electron subspace there are two alternative c Fermi momenta. Those are associated with the two 
doublicity values d = ±1, respectively. Their expression is given in Eq. ([TU|. Specifically, the doublicity d — ±1 of 
a c fermion created or annihilated under an one-electron excitation equals that of the corresponding si fermion hole 
created under the same excitation. 

The concept of doublicity remains valid when the angle tTcs/'^ between the momentum-space directions of qp^ and 
'is si appearing in the (j)^ expression of Eq. (|11|) is (p dependent, provided that the following integrals involving the 
corresponding function tTcs = t^cs{4') and the c Fermi line hole momentum and si boundary line momentum absolute 
values q%c{4>) and qBsi(0), respectively, vanish. 



7r/4 



d(j)qp^{(j))qBsi{4>)cos 







T^cs{4>) 



7r/2 



7r/4 



d<t>q^pM<lBsi{4>)cos(^^\=Q 



(12) 



This ensures that the Fermi line centered at —7? encloses the correct momentum area, 

rg7r[4(</>)]^ = £^g,(^[4^(^)]%[,^,,(0)]V2g^,(0)te,i(0)cos(^^ 
= X 47r2 + (1 - x) l-K^ = (1 + x) 2-K^ . 



(13) 



Here x47r^ and (1 — a;) 27r^ are the momentum areas enclosed by the c Fermi line and si boundary line when centered 
at — TT and = [0, 0], respectivelyi. In this paper we consider the general case for which the angle 7rcs/2 = 7rcs(0)/2 
is for X g [xcIt'^c'i) a periodic function of </> of period 7r/2. In addition, for e (0, 7r/4) it obeys the equality 
7rcs((/))/2 — 'Kcsij^l'^ — '/>)/2 where 7r/2 — € (7r/4,7r/2). It follows that for x £ {xci,Xc2) the one-electron F angle 
4>p{(j)) appearing in Eq. ([TT]) reads. 



6f^((/)) = darctan 



g^,(0)sin(2^) 
^9B.i(</') + g^,(0)cos(^) 



d = ±l. 



(14) 



It can have two values. 



^-/ 



and 



The two corresponding rotations refer to the doublicity d = — 1 and d = +1, 



respectively. For tTcs/2 = tt/2 this general expression recovers that found in Rcf.^. The F angle 4>p{<j>) is associated 



with the one-electron F rotation matrix Ap such that, 



QBsl 



A%qBs 



Ad 



cos (j)p — sin (j)p 
sin (j>p cos (j)p 



d = ±l 



(15) 



Alike in Refii, we denote the auxiliary momentum of the si boundary momentum qg^^ by qBsi- For x € {xci,Xc2) 
and g at or near the si boundary line the matrix A'^^ of Eq. ([7]) is orthogonal and equals the F rotation matrix Ap 
given here. 

If follows from the above expressions that for the hole-concentration range x £ {xci,Xc2) the absolute value kp{(f>) 
of the hole Fermi momentum given in Eq. (jlOp and Eq. (|A6|) of Appendix A reads, 



kFW^\[<lkW? + kBsim' 



2q'^^{(p)qB si {(I)) COS 



.w 



(16) 



Again, for tTcs/^ = tt/2 this recovers the corresponding expression of Ref.^. The hole Fermi momentum kp of Eq. 



([TOl) for any value of the Fermi angle 
as follows. 



and the nodal Fermi momentum kp 



ZhN 



7? for (j) = 7r/4 can be expressed 



k^{^) 



7Tt9b^i(W: 

qBsiW 



kp 



I.N 
t^p 

yssi 



?li 



qBsl[(f>) = qBsl{TT/2- 



(17) 



Here kp{(t)) and kp are centered at —7? and 0, respectively, and kp is the absolute value of the vector kp . Note 
that the auxiliary si boundary momentum qBsi of Eq. (jlSp has been constructed to inherently pointing in the 
same direction as kp. This justifies the validity of the equalities (J17p . which involve the auxiliary si boundary-line 
momentum rather than the corresponding si boundary- line momentum. 

The ratios Ts and re of Eq. ljAip of Appendix A, which for C//4t e (uo, ui) are related as re = ttx* w 2rs, control the 
effects of electronic correlations in many quantities. For instance, for such a U /At range they control the range width 
max(7rcs/2)-min(7rcs/2) = r^ of the angle tTcsI'^- Indeed, 7rcs/2 S ([1 — xa\ 7r/2, [1 -|- xa] t^I'^)- Here xa ~ a:*/2 = 0.135 
is a crossover hole concentration slightly larger than Xd — 1/8 = 0.125. According to the studies of Ref<^ it marks 
the emergence of a scale invariance in the VEP quantum liquid, which dominates the physics in a hole concentration 
range x € {xa-,Xc2)- The smallest and largest 7rcs/2 magnitudes are reached for hole Fermi momenta pointing in the 
nodal and anti-nodal directions, 

TTes{TT/A)/2=[l-XA]TT/2; TT,, (0)/2 = TT,, ( ^2) = [1 + -T^] 7r/2 ; XA^Y^ (1^) 

respectively. For U/At « 1.525 this gives 7rcs(7r/4)/2 « 0.43 tt and 7rcs(0)/2 = 7rcs(7r/2)/2 « 0.57 tt. Hence 7rcs/2 S 
(0.43 TT, 0.57 tt) has indeed values near 7r/2. Both this and that its average value is 7r/2 justifies the approximation of 
Ref.^ that it reads 7r/2. 

Except for the angle tTcs/^ = tt/2 being replaced by a dependent angle tTcs/2 G ([1 — xa] 7r/2, [1 -|- xa] 7r/2), the 
expressions and physics reported in Refi^ remain valid. Moreover, the c fermion energy dispersion edq ) provided in 
Eq. (JAIOP of Appendix A is not changed by the VEP quantum liquid superconducting fluctuations studied in this 
paper. In contrast, the si fermion energy dispersion esiio) is, as found in the following. 

III. LONG-RANGE d-WAVE SUPERCONDUCTING ORDER 

Here evidence is provided that the ground state of the VEP quantum liquid is superconducting for zero spin density 
m = 0, hole concentrations in the range x G {Xc,x^,), and intermediate U/At values. We recall that for approximately 
U/4:t > uq the VEP quantum liquid refers to the Hamiltonian H^jj of Eq. ([5|) in the one- and two-electron subspacei 
for t±/t <C 1 plus the very small suppression effects considered below. For such a U/At range of values, that the 
inequality t±/t <^ 1 holds assures that the anisotropy parameter e'^ = m*^/M is very small as well. 

We start by considering the ground state of the 7d = 1 and t^/t = square-lattice quantum liquid of Refi at 
vanishing spin density to = 0. We find that a necessary condition for it being that of a d-wave superconductor is 
fulfilled for finite hole concentrations below x* . Such a result is valid for the VEP quantum liquid as well provided that 
the only significant effects of the weak 3D uniaxial anisotropy and suppression effects are the emergence of the small 
anisotropic parameter e^ = to* jM and parameter 7^ « 1 in the expression of several physical quantities associated 
with quantum and thermal fluctuations. For the five representative systems the parameter e^ — m*/Al is found in 
Refji^ to belong to the range e^ € (3 x 10^^, 1 x 10^^). For such a small values of that parameter and 7d ~ 1 the 
energy eigenstates of the 7^ = 1 and t±/t — square-lattice quantum liquid studied in Ref^i are expected to be a 
good approximation for those of the VEP quantum liquid. 

Consistently with the Mermin-Wagner-Berezinskii Theorem^ir— , the introduction of the weak 3D uniaxial 
anisotropy prevents the destruction of the fluctuations of long-range orders at finite temperature. For t±/t <^ 1 
very small and approximately U /At > uq we find indeed strong evidence of virtual-electron-pair phase coherence 
below some critical temperature Tc- (The concept of a virtual electron pair is introduced below.) At zero temperature 
such a virtual-electron-pair phase coherence occurs for the hole-concentration range x € (xccc*). Here the lower 
critical concentration Xc ^ Gi + xq where Gi — G [AZ/to*] = G/e^ is fully determined by the anisotropic parameter 
£^ = ml/M ^ 1 and xq emerges due to the hole-trapping effects of Appendix B. The magnitude of the proportionality 
constant G is found in Ref.^^ to be in the range G G (1 x 10~^, 3 x lO^"') for the five representative systems. Hence it is 
so small that Gi « 10"^ in spite of m^/M ^ 1 being much smaller than Gi. The weak 3D uniaxial anisotropy allows 
considering the pairing energy within a small coherent volume. Such an elementary volume controls the magnitude 
of parameters associated with the thermal and quantum fluctuations of the square-lattice quantum liquid perturbed 
by weak 3D uniaxial anisotropy. The suppression effects slightly lessen the Tc magnitude, which is linear in 7^ ~ 1. 

The physical picture that emerges from our study is that superconductivity arises in the VEP quantum liquid as a 
by-product of the short-range spin correlations. We then construct a consistent scheme concerning the d-wave pairing 
mechanism, phase-coherent-pair superconducting order, and corresponding order parameter, which follows from the 
properties of that quantum liquid. Our results are inconclusive on whether the ground state of the 2D Hubbard model 
on the square lattice is superconducting. They seem to indicate that some small 3D uniaxial anisotropy is needed for 
the emergence of superconductivity. 
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A. A necessary condition for the ground state at zero spin density being that of a d-wave superconductor 

Let us consider the t^/t — and 7d = 1 Hubbard model on the square lattice ([IJ. The one- and two-electron 
subspace is spanned by states whose deviation 5N^ in the number of c band holes and number N^-^ of si band holes 
read", 

5Nl' - -6N = 0, Tl, T2 ; Ni\ = ±{5N^ - 6N^) + 21,,^^^ + 2Ns2 - L^, _i/2 + L,^ +1/2 + 2iV,2 = 0,1,2. (19) 

Here SN = [SN^ + 6N^] is the deviation in the number of electrons, SN^ and SN^ those in the number of spin-projection 
t and 4, electrons, respectively, Ns2 the number of the excited-state s2 fermions, and is, ±1/2 that of independent 
spinous of spin projection ±1/2. 

For finite hole concentrations below x* and intermediate C//4t values the use of the energy functional introduced 
in Refji reveals that excited states with si band hole numbers iV^^ — and A^^^ = 1,2 and involving addition to or 
removal from the c Fermi line of c fermions refer to a gapless excitation branch and have an energy gap, respectively. 
For the latter excited states the energy gap vanishes only if the auxiliary momentum of the si band hole {N^i = 1) 
or both si holes {NJ^i = 2) points in the nodal directions. Use of Eq. (TTOl) then reveals that addition or removal 
of (i) one electron and (ii) two electrons with the same spin projection to or from the hole-like Fermi line whose 
hole Fermi momenta are given in Eq. (|10p involves creation of (i) one and (ii) two holes, respectively, in the si 
momentum band. Except for the excitation-momentum nodal directions these excited states have a finite energy gap. 
In contrast, note that addition or removal of two electrons of opposite spin projection to or from that Fermi line leads 
to a final excited state whose si band is full alike that of the initial ground state. Indeed such excitations correspond 
to Ns2 = is,=Fi/2 = in Eq- (HH), so that Nj}^ = for SN^ = SN^ = ±1. Hence the latter processes refer to a gapless 
branch of two-electron excitations. Specifically, for excitations whose si holes are created at the si boundary line the 
energies of such processes read to first order in the c and si hole momentum distribution-function deviations, 

6E = SEpiip) = I All cos 201; 5N^ =±1, 5N^„^Q, 

5E = 5EF{(t)) + 5EF{(t)')^\l^\[\cos2(l)\ + \cos2(J3'W- 5N^=±2, 5N^„^Q, 

5E = 0; 5N^^ SN^ = ±1 . (20) 

Such general spectra refer to vanishing spin densities m — and approximately U/At e (uojUtt). (For the square- 
lattice quantum liquid of Ref."'^ and the VEP quantum liquid they refer to a; G (0, x*) and x e (xc, x*), respectively.) 
In the expressions provided in Eq. ((20)) . SEp — \Asi{qj^gi)\ is the anisotropic Fermi-energy term and |Asi((fJgj^)| is 
the si fermion pairing energy per spinon given in Eq. (JA12[) of Appendix A. The Fermi energy reads Ep = M + SEp 
where /z is the zero-temperature chemical potential. The anisotropic Fermi energy SEp vanishes for x > x^, since then 
TaaxISEp} = |A| — 0. It follows that for the hole concentration range x S {x^,, 1) the ground state is a disordered 
state without short-range spin order. Consistently, for that range of x values the si fermion spinon-pairing energy 
2|Asi(gg^j^)| vanishes for all momentum values. 

The numbers of si fermions and si fermion holes in the si momentum band equal those of occupied and unoccupied 
sites, respectively, in the si effective lattice. Hence in the following we discuss the structure of the different energy 
spectra SE of Eq. (|20p in terms of the numbers of sites, occupied sites, and unoccupied sites of that lattice. Those 
equal the corresponding numbers of discrete momentum values, filled discrete momentum values, and unfilled discrete 
momentum values, respectively, of the si fermion band. That the energy SE of Eq. (j20p vanishes for excitations 
involving creation of two electrons of opposite spin projection and except for the nodal directions is for x € (xc,a;*) 
finite and given by SE = [SEp{4>) + SEp{(f>')] for those involving creation of two electrons with the same spin projection 
is not a trivial result. (The additional term to 6Ep = \Asi{qQgi)\ found below in Section IV-A as a result of the 
superconducting fluctuations does not change the basic property that SE = and SE > for excitations involving 
creation of two electrons of opposite spin projection and the same spin projection, respectively.) That behavior 
follows from the number N^ of sites of the si effective lattice being for the Hubbard model on the square lattice a 
subspace-dependent functional^ . For the one- and two-electron subspace referring both to the square-lattice quantum 
liquid of Ref.^ and the VEP quantum liquid, the expressions of the number deviations SN^^ and SNsi of sites and 
occupied sites, respectively, of the si effective lattice read, 

SNl^ = SN^ + L,^ _i/2 ; SN^i = SN^ - L,^ _i/2 - 2iV,2 . (21) 

Such expressions are consistent with that of iV]\ = [SN'^^_^ — SNsi]^ given in Eq. (|19p. For the deviation numbers 
and numbers SN^ = SN^ = ±1 and is, ±1/2 = Ns2 = of the SN = ±2 excited states for which the two added 
or removed electrons are in a spin-singlet configuration one finds from the use of the expressions provided in Eqs. 
P^ and (pi]) that the corresponding deviations in the numbers of sites and occupied sites of the si effective lattice 
read SN^^_^ — ±1 and i^A^si = ±1, respectively. Therefore, the deviation in the number of unoccupied sites vanishes, 
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(5iVJ\ = [SN^^_^ — SNsi] = 0. Indeed, under such excitations the creation or annihilation of one si fermion is exactly 
cancelled by an increase or decrease, respectively, in the number of sites of the si effective lattice. As a result the 
number of unoccupied sites of the corresponding excited states remains being zero, as for the initial ground state. 
Use of the energy functional introduced in Ref.^ and given below in Section IV- A then leads to the gapless spectrum 
SE = OoiEq. 1^. 

In turn, for the deviation numbers and numbers 6N^ = ±2, SN^ = 0, Ls±i/2 = 2, and is,ipi/2 = Ns2 = or 
SNi = ±2, 5N^ = 0, Lg^ ;pi/2 = 2, and L^ ±1/2 = A^s2 = of the 5N = ±2 excited states for which the two added or 
removed electrons are in a spin-triplet configuration one finds from the use of the expressions provided in Eqs. p^ 
and ([2T|) that the corresponding deviations in the numbers of sites and occupied sites of the si effective lattice read 
SN^ — [±1 + 1] and SNgi = [±1 — 1], respectively. It then follows that the deviation in the number of unoccupied 
sites are given by SNl^i — [SN^_,^ — SNsi] — 2. In contrast to the above excitations involving a spin-singlet electron 
pair, the creation or annihilation of one si fermion is not cancelled under the latter excitations by an increase or 
decrease, respectively, in the number of sites of the si effective lattice. As a result, its number of unoccupied sites 
increases from zero to two for the excited states. Use of the general energy functional of Ref.>^ then leads to the 
gapped spectrum given in Eq. ([20]) for these excitations. (Use of the modified VEP quantum liquid functional of 
Section IV-A leads to a gapped spectrum including an additional term, which as mentioned above does not change 
the physics discussed here.) 

A similar analysis can be carried out for the 6N = ±1 electron and 6N = spin excitations. Their energy spectrum 
is in general gapped. The d-wave-like structure of the one- and two-electron spectra of the Hubbard model on the 
square lattice given in Eq. (j20p follows from the momentum dependence of the si fermion dispersion studied in Refji. 
Whether the general energy spectrum (j20p is gapless or displays a gap is fully controlled by the interplay between the 
deviations in the numbers of sites and occupied sites, respectively, of the si effective lattice. That lattice is exotic in 
that its number of sites is subspace dependent. The spectrum (pOl) refers to a gapless branch of excitations whenever 
the corresponding creation or annihilation of si fermions is exactly cancelled by an increase or decrease, respectively, 
in the number of sites of the si effective lattice so that its number of unoccupied sites remains being zero as for the 
initial ground state. Above we gave the example of creation or annihilation of a spin-singlet electron pair. Such a 
canceling occurs for creation or annihilation of any finite number of such spin-singlet electron pairs. 

In turn, for one-electron excitations and creation or annihilation of spin-triplet electron pairs the spectrum of 
Eq. (I20p is in general gapped except for some momentum directions. The d-wave-like structure of that spectrum is a 
necessary condition for the ground state of the model being for vanishing spin density to = 0, finite hole concentrations 
X € {xc,Xi,), and thus finite 2|A| that of a d-wave superconductor. However, it is not a sufficient condition for 
the occurrence of coherent pairing needed for the macroscopic condensate. Below we find strong evidence that for 
intermediate U/it values, vanishing spin density to. = 0, and the hole concentration range x G {xc,x^.) the sufficient 
condition of phase coherence is met by the VEP quantum liquid. At zero temperature and both for hole concentrations 
in the ranges x £ (xq^Xc) and x € {xc,x^,) there is short-range spin order. For x € {xo,Xc) pairing correlations occur 
that due to strong phase fluctuations do not lead to coherent pairing and superconductivity. 

Finally, we emphasize that the number deviation and number expressions provided in Eqs. (|19p and (j2ip also 
hold for I? = 1 spatial dimensions with N^^_^ replaced by Na^^^- However, for the ID Hubbard model there is no 
short-range spin order for any range of x values so that | A| = in Eq. (PU]) and SE = for all one- and two-electron 
excitations under consideration. Consistently, the ground state of the ID model is not superconducting. Such an 
analysis reveals that a necessary condition for the occurrence of a superconductivity order in the VEP quantum liquid 
is the occurrence of short-range spin order associated with the finite energy parameter 2|A| > 0. That strongly 
suggests that the occurrence of superconductivity in such a system is a by-product of the short-range spin correlations 
associated with the energy scale 2|A|. However, our results are inconclusive on whether for the Hubbard model on 
the square lattice the phase coherence needed for the occurrence of long-range superconducting order does occur. 

B. Short-range spin order, the pseudogap energy scale, and the pseudogap temperature T* 

Here we provide some basic information about the energy scales associated with the short-range spin order, which 
for vanishing spin density to = and hole concentrations in the range x € {xq , x* ) refers to finite temperatures below 
a pseudogap temperature T* . Indeed, the VEP quantum liquid scheme accounts for the hole trapping effects reported 
in Appendix B, so that the short-range spin order occurs for the range x G (a;o,x«) rather than x S {0,x■^,) for the 
square-lattice quantum liquid of Refj^. Evidence is provided in Ref^iS that for the hole concentration range x € (xc, a;*) 
the suppression effects due to intrinsic disorder or superfluid density anisotropy are for the four representative systems 
other than LSCO very small. In contrast, the LSCO cation-randomness effects introduced in RefJ^ are not small. 
Fortunately, they are merely accounted for by multiplying the energy scale Aq and related magnitudes by a factor of 
1/2, leaving the pseudogap temperature T* considered in the following unaltered. Its magnitude remains unaltered as 
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well under the suppression effects. The results of this paper and Ref.^'^ confirm that our oversimplified description of 
the effects intrinsic disorder or superfluid density anisotropy and LSCO randomness in terms of suppression effects and 
cation-randomness effects, respectively, leads to agreement between theory and experiments on the five representative 
systems. 

At zero temperature the energy order parameter 2|A| of the short-range spin correlations is the maximum magnitude 
of the si fermion spinon-pairing energy of Eq. (jA13|) of Appendix Aa^. Upon increasing x the residual c - si 
fermion interactions tend to suppress the short-range spin correlations. For approximately U/At € (uo,U7r) and hole 
concentrations x £ {xo,x^) such an effect leads to the linear decreasing 2| A| « (1 — x/x^,) 2Ao of the zero-temperature 
order parameter of the corresponding short-range spin order. That energy parameter plays the role of pseudogap 
energy scale, as it controls the magnitude of the pseudogap temperature T* above which there is no short-range spin 
order, associated with si fermion spinon pairing. However, there are several definitions of T* corresponding to effects 
of the pseudogap associated with 2|A| on different physical quantities. Indeed, such effects may appear at different 
temperatures. 

For the square-lattice quantum liquid of Rcf.^, for which the short-range spin order occurs for < x < x^, we 
identify in the limit x — >■ the temperature T* with the temperature Tx of Ref.— . In that limit it is related to the 
energy parameter Aq = lim2,_i.o |A| as T^ ~ Ao/fcs^. For the quantum problem considered here and approximately 
U/4:t € (uo,u^) the corresponding expression valid for finite hole concentrations in the range x € {xo,x^) is then, 

T*^(l-^]^. (22) 

This is actually the maximum magnitude that T* can achieve. Its magnitudes as obtained from different properties 
should obey the inequality T* < {1 — x/x^,)[Ao/kB]- Hence the zero-temperature magnitude of the short-range 
spin order parameter 2|A| « (1 — x/x^.) 2Ao controls the range of the pseudogap temperature T* « 2|A|/2fcB. For 
temperatures above T* the system is at zero spin density m = and for finite hole concentrations in the range 
X £ (xo,x*) driven into a spin disordered state without short-range spin order. 

C. Selected Hamiltonian terms: The VEP quantum-liquid microscopic Hamiltonian 

For the hole concentration ranges x € {0,Xc) and x G {xc,x*) the quantum problem considered in this paper is 
qualitatively different. For x G (0,Xc) it refers to that reported in Appendix B, for which the effects of intrinsic 
disorder are strong. For x € (xc, x*) the VEP quantum-liquid suppression effects are found in Refj^ to be very small. 
That under the LSCO cation-randomness effects introduced in that section the energy scale Aq is lessened by a factor 
two is behind the hole concentration xq of Eq. ([5]) being different for the parameters appropriate to LSCO and the 
remaining four representative systems: It reads xq « 0.013 and xq « 0.024, respectively. 

Expression of the t±/t <^ 1 Hamiltonian H^jj of Eq. ([5|) in terms of rotated-electron creation and annihilation 
operators leads to an infinite number of termsi. Although in terms of electron operators the Hamiltonian has only on- 
site interactions, in terms of rotated electron operators there emerge effective interactions involving rotated electrons 
on different sites. Only a small number of such Hamiltonian terms are relevant to the physics of the present quantum 
problem. However, the appropriate selection of the latter Hamiltonian terms is a problem of huge complexity. 

Here we use as criterion for that selection the general strongly correlated microscopic Hamiltonian that, according 
to the analysis of Ref<^, almost certainly underlies the essential physics of the representative hole-doped cuprates. For 
U/it > uo this leads to a simplified effective Hamiltonian with the same general form as that given in Eq. (1) of Refj^, 
but with the electron creation and annihilation operators replaced by corresponding rotated electron operators. Its 
kinetic-energy terms include those generated from the operator H± given in Eq. ([5]). In Section IV- A we introduce 
the general energy functional corresponding to such a microscopic Hamiltonian. The small hole concentration denoted 
by Xq in Ref."*, which reads xq ~ 0.01 for LSCO, is identified here with the hole concentration xq of Eq. ([5]). For 
the VEP quantum liquid a:o is the critical hole-concentration at which there occurs a sharp quantum phase transition 
from the Mott-Hubbard insulator with long-range antiferromagnetic order to a short-range incommensurate-spiral 
spin ordered state. That in RefJ^ it is found to read xq ~ 0.013 for LSCO and as given in Eq. ^ is proportional to 
Ao/t is consistent with the results of Ref.''. 

The selected Hamiltonian terms of our microscopic Hamiltonian are basically the same as those of the microscopic 
Hamiltonian (1) of Ref 4, with the electron operators replaced by rotated-electron operators and without the term 
containing the parameter Up. Indeed, after expression of the rotated-electron operators in terms of c and si fermion 
operators the one- and two-electron subspace with no rotated-electron double occupancy is well-defined: The gener- 
ators of the states that span such a subspace have simple expressions in terms of c and si fermion operators^. As 
a result there is no need of introducing artificial Hamiltonian terms to impose the lack of double occupancy. In our 
case that is achieved by defining the microscopic Hamiltonian in the one- and two-electron subspace. Indeed, the 
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corresponding VEP quantum liquid is defined in that subspace. Below we express the terms of the microscopic Hamil- 
tonian that control the thermal and quantum fluctuations in terms of c and si fermion operators. The kinetic-energy 
terms of the VEP quantum-liquid microscopic Hamiltonian involve both the C//4i-dependent c fermion mass m* and 
the much larger mass M of Eq. ([5]) . 

The advantage of the rotated-electron operator description of Ref.^ is that it has been constructed to inherently single 
and double rotated-electron occupancies being good quantum numbers for U/At > 0. Hence for such operators the lack 
of rotated-electron double occupancy is exact for the range U/At > uq of the VEP quantum liquid. The transformation 
that relates rotated electrons to electrons is unitary, so that the effective action for the phases 9 considered in this 
paper is valid for approximately t//4t > uq and specifically at U/At fn u^, — 1.525. This is important in view of the 
results of RefSfi^i?'^. Such results provide evidence that within the description of the properties of several classes of 
hole-doped cuprates with Xc ~ 0.05 and x^ sa 0.27 by the VEP quantum liquid the appropriate value of U/At is not 
large and reads C//4t « u* = 1.525 € (mo,wi). This includes the five representative systems. The intermediate value 
U/At w 1.525 is also that found in Reffi to be appropriated for the description of the spin- wave spectrum of the a; = 
parent compound La2Cu04 (LCO)^^. 

Most of the VEP quantum-liquid microscopic Hamiltonian terms refer to in-plane processes. The very small t±- 
dependent terms generated from the operator H± given in Eq. ([5]) control the magnitude of the hole concentration 
Xc and thus have effects on the virtual-electron pair phases 9jfl and 9j^i introduced below. In the superconducting 
phase considered in the following such terms allow the Josephson tunneling through nearest-neighboring planes of 
a very small density of vanishing-energy spin-singlet electron pairs. For each square-lattice plane such excitations 
correspond to the gapless excited states of the general spectrum (j20p . which refer to creation or annihilation of spin- 
singlet electron pairs. The small average numbers of such pairs that leave and arrive to a given square-lattice plane are 
identical. The main role of that tunneling is to allow for thermal and quantum fluctuations associated with in-plane 
long-range superconducting order at finite temperatures. 

The general energy functional introduced in Section IV-A corresponds to the VEP quantum-liquid microscopic 
Hamiltonian, yet it has an in-plane character. It includes implicitly the needed 3D uniaxial anisotropy effects. 
Indeed, part of such effects are effectively described by the use of a suitable mean-field approximation for the square- 
lattice physics. Specifically and alike in the microscopic Hamiltonian of Ref."*, in some of the selected Hamiltonian 
terms studied below a complex gap function Aj/j// replaces the corresponding pairing operator. Indeed within mean- 
field theory finite-temperature long-range orders may occur in the 2D system. Although this is excluded by the exact 
Mermin-Wagner-Berezinskii Theorero^l^— , mean-field theory refers to an additional effective way to indirectly account 
for the effects of the very small spin-singlet electron pair Josephson tunneling through nearest-neighboring planes. 
In turn, in the pseudogap state considered below the small i_L-dependent terms are behind very small one-electron 
transfer between first-neigbhboring planes. For each square-lattice plane such excitations correspond to excited states 
of the general spectrum (PH)) with a finite energy gap. Those refer to creation or annihilation of a single electron. 
As discussed in Ref-^^.^ such excitations are behind the energy gap of the normal-state conductivity in the direction 
perpendicular to the planes. 

Before introduction in Section IV of the general energy functional corresponding to the VEP quantum-liquid 
microscopic Hamiltonian, in the following we express the creation and annihilation rotated-electron operators of 
the terms of such a Hamiltonian that control the thermal and quantum fluctuations in terms of c and si fermion 
operators. Analysis of the form of the obtained effective Hamiltonian terms reveals that the spin bonds and the 
charge-2e sector of Refi^ correspond to the spin-singlet two-spinon si fermions and c fermion pairs, respectively. For 
the hole concentration range x € (xc,a;*) of the VEP quantum liquid both the effects of 3D uniaxial anisotropy 
and intrinsic disorder are small. Hence for such a x range the interactions of these objects are within our description 
residual. This results from the residual character of such interactions within the starting square-lattice quantum liquid 
of Reffi. That property greatly simplifies the derivation of the one-electron scattering rate carried out in Refj^. 

The d-wave-likc structure of the Fermi-line one- and two-electron energy spectrum given in Eq. (j20p is only a 
necessary condition for at zero spin density m = and finite hole concentrations in the range x G {xq, a;*) for which 
|A| > the ground state of the present quantum problem to be that of a d-wave superconductor. That energy 
spectrum refers to e^ = m*^/M — >■ 0. It also applies here provided that the 3D anisotropic parameter e^ = m*^/M is 
very small and 7d ~ 1- In the remaining of this section we access the x dependence of the T = parameters associated 
with the VEP quantum-liquid thermal and quantum fluctuations for the approximate range C//4i € {uq^Ut^). Indeed 
some of our expressions are not valid for U /At > u^. 

In the following evidence from the study of pairing phase fluctuations is found that for zero spin density m = Q and 
hole concentrations in the range x € {xc, a;*) the ground state of the VEP quantum liquid has phase coherence, which 
is associated with a long-range superconducting order. While our results refer to U/At e (uojUtt), this may hold for 
U /At > Utt as well. Our results also indicate that for vanishing spin density to = and finite hole concentrations in 
the range x £ {xq, Xc) there is short-range spin order and pairing correlations yet the lack of phase coherence prevents 
long-range superconducting order. Finally, it is expected that alike for the system of Ref.^, for m — and x £ (0, xq) 
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monopole-antimonopole pairs of the type considered in that reference unbind and proliferate, leading to long-range 
antifcrromagnetic order, consistently with the analysis of Appendix B. 

1. Hamiltonian terms that control the quantum fluctuations of the phases associated with competing orders 

Out of the selected terms of the VEP quantum-liquid microscopic Hamiltonian, here we consider those that control 
the quantum fluctuations of the phases associated with the competing orders. It is assumed that the c and si fermion 
occupancy configurations that generate the energy eigenstates of the t±/t = and 7^ = 1 problem generate states 
close to energy eigenstates for the VEP quantum liquid at very small tj^/t <C 1 values and 7d ~ 1- Such terms are 
given by, 

N/2 
H^onas ^ ^ ^ ^^.,^.„ [gt^^^ ^ ~^^^__ ^ _ ,t^^^ ^ ,t^^^^ ^^] -, (h.C.) , (23) 

consistently with the maximum number of in-plane independent bonds being iV/2. Indeed, for hole concentrations 
below Xtf and zero spin density m = the ground states are spin-singlet stateai. For x < xq and x e (xo,^*) 
such states have long-range and short range spin order, respectively. They contain TV^i = iV/2 two-spinon si bond 
particlesi^. The two spinous of such spin-neutral objects describe the spin degrees of freedom of rotated electrons 
that singly occupy sites in the ground-state configurations. Therefore, in the one- and two-electron subspace without 
rotated-electron doubly occupancy where such ground states are contained there is an energetic preference for the 
formation of spin-singlet rotated-electron bonds. 

The summation ^ -, „//r,_constl "^^ ^^"^ right-hand side of Eq. (P5|) is over all in-plane iV/2 spin-singlet rotated- 
electron bonds considered in Ref.— . Alike in Ref.— , to arrive to the Hamiltonian terms given here we used the 
usual mean-field approximation within which the complex gap function ^y y has replaced the corresponding pairing 
operator. The weak effects of the very small 3D uniaxial anisotropy perturbation iJj^ of Eq. ([5]) occur on the 
quantum liquid in-plane Hamiltonian terms (|23p through the phase of the complex gap function ^y j" , as discussed 
below. Each spin-singlet rotated-electron bond is centered at a real-space point [r*,-/ -I- rj"]/2, near that of coordinate 
"^3 ~ i^j' + ^j"]/2 ~ ^ [as/2] e^d- Here the two-site bond indices I ± 1 and d — 1,2 are those used within the notation 
of Ref.^^ discussed below and Qs — a/y/l — x is the spacing of the spin effective lattice. (The bond index d = 1, 2 is 
unrelated to the doublicity d — ±1 considered in Section II and Appendix A.) 

The rotated-electron creation and annihilation operators of the Hamiltonian terms (I23p act onto the one- and two- 
electron subspace with zero rotated-electron double occupancy^. For finite hole concentrations in the range x G {xq , x* ) 
such terms refer to energy scales below and around the energy Aq of Eq. (jB)). The complex gap function is defined on 
the two-site bonds in the spin effective lattice studied in Ref.^^.. The energy scale 2Ao equals the absolute maximum 
excitation energy below which the short-range spin order survives. Therefore, Aq is the absolute maximum magnitude 
of the pairing energy per spinon of the two-site and two-spinon bond associated with the complex function Aj/jn—. 
The amplitude of Aj'j" is for finite hole concentrations in the range x G (a;o,a;*) frozen below the energy Ao/v2. 
Here l/v2 is a suitable normalization factor, which is absorbed in the expressions of the two-site and two-spinon 
operators of the si bond-particle operators considered below. What remains are the fluctuations of the phases 9j'j" 
of rotated-electron pairs. The complex gap function reads. 






^3' 3" = ^ '''" T^^O ; Xo = ' V C'o ; Co = I ,„ ,, , \ ,, I « 77:77^ , Xo<X<X^ 



^3' 3" = ^''''' -^^ , < X < xo . (24) 

The hole concentration xo of Eq. ([8]) is here expressed in terms of |Aj'j//| = Ao/-\/2. It corresponds to the hole 
concentration also called xq in Ref.''. The studies of that reference estimated it to be proportional to |Aj/j//|/t and 
such that xo — [|Aj/j"|/t]Co <C 1. Its expressions given in Eqs. (jH]) and (p4|) are consistent with such an estimation. 

Following the d-wave character of the spinon pairing of the si bond particles, one has in the expressions of Eq. 
(1211) that d = 1 and d = 2 for the families of spin-singlet two-site bonds whose primary bonds considered in Refi^ 
are horizontal and vertical, respectively. The different magnitudes Aq and /i''/2 that the energy scale |Aj/j//| has for 
finite hole concentrations in the ranges x € (xq, x*) and x e (0, xg), respectively, are due to the sharp quantum phase 
transition occurring at x = xq. (We recall that f/' is the energy scale given in Eq. (|A3p of Appendix A.) 

Usually one restricts the in-plane summations ^ , „ of Eq. (P5)) to nearest neighboring sites. In turn, the terms 
given in that equation involve contributions from all possible in-plane bonds whose two rotated electrons are located 
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at arbitrarily distant sites. In order to capture the physics of the quantum problem studied here one must start by 
taking into account all such contributions. After some algebra one then arrives to an effective Hamiltonian, which is 
equivalent to restricting the in-plane summations ^ -, ■„ of Eq. ([23| to nearest neighboring sites. 

The importance of the selected Hamiltonian terms (|23|) is that they contain the phases Oj'j" whose fluctuations 
control the physics of the VEP quantum liquid. It is useful to express (j23p in terms of c fermion operators and 
two-site and two-spinon bond operators. This is straightforwardly achieved by direct use of the expressions of the 
rotated-electron creation and annihilation operators in terms of c fermion operators and two-site and two-spinon bond 
operators^, with the result, 

A//' [ci, ,, £t ,„, - £t „, £t „,,] = (- 1) i-"^ 72 A,,,. fl,jl,A,,,,si,,,,,- (25) 

Here "Tj'j" = [rj' +rj"]/2 = tj + I [as/2] e-c^, a^ = a/\/l — x for x < x*, fot^-^ ^ , is the in-plane two-site bond operator 
defined in Ref.^'^, and the index g — Q,\, ..., [Nsi/4: — 1] refers to the link or bond type also defined in that reference. 
(The index g appearing here is not the amplitude g introduced below in Section III-D.) 

A si fermion operator of real-space coordinate rj is defined in terms of a superposition of such bond operators, 

Nsi/4:-l 2 

4,.i = e*- E EE^;^I-.+.o,..M,...' (26) 

g = d=ll = ±l 



where (j)j,si is the operator phase of Eq. ([2|). The absolute value |/ig| = \/h*hg of the coefficients appearing in this 
expression decreases for increasing magnitude of the two-site bond length ^g = |2rJJ. The minimum and maximum 
values of the length ^g are ^q = a^ and (,g = ^/2as (A^si/4 — 1) for g = [Nsi/4: — 1], respectively. These coefficients 
obey the normalization sum-rule X]g=o I'^sP = ^/'^- Each of the Ngi spin-singlet two-spinon bonds of a si 

bond particle of real-space coordinate rj involves two sites of coordinates r ~ r^i and r -\- f^i, respectively, where 
rj = r — r^i ^-nd r^i = I [as/2] e^j^. It follows that the two-site bond centre f = rj + r^^ is the middle point located 
half-way between the two sites. There are four families of bonds labeled by the numbers d = 1,2 and Z = ±1. For 
each family there are Nsi/A link vectors r^i of different link type g = 0, 1, ..., [A^si/4 — 1]. 

For simplicity, in Eq. (j25p we call rji and rjn the two real-space coordinates of the sites of a two-spinon bond. 
According to that equation, those are also the real-space coordinates of the two c fermions, respectively, involved in 
the corresponding spin-singlet rotated-electron pair. However, we recall that the two real-space coordinates of such a 
bond and corresponding c fermion pair refer to well-defined indices d, I, and g. The g = primary bonds and primary 
c fermion pairs have most of the corresponding rotated-electron pair spectral weight. For each of the 7V/2 in-plane 
real-space coordinates fj = [fji + fjii]/2 — r^ i there are four primary c fermion pairs corresponding to d = 1 horizontal 
and d = 2 vertical bonds and I = —1 left or lowest and / = +1 right or upper bonds. 

^j'j" 



The phase factor e i'^" of the complex gap function ^j'j" appearing in Eq. (|25p can be written as, 



g»^'," =g«w,',".oe'^.'.".i . (27) 

Here dj'j"fi corresponds to the overall centre-of-mass phase. In turn, Oj/jn^i is the part of Oj'jn, which cannot be 
reduced to arbitrary configurations of site phases and regulates the relative motion of a pair. Hence the phases ^j'j",i 
are related to the internal pairing degrees of freedom. Within charge excitations the phases Ojijn are associated with 
the c fermion pair of Eq. (1251) . For one-electron and spin excitations they are associated with the virtual-electron pair 
defined below. 



Use of expression (|24p in Eq. ((25|) leads to, 

A,'." [£t , ^, 5^„ ,, - £^, ,, 5^„ ,,] = e^^. v Ao IIjI^„,A,,,^sia- (28) 



Hence one finds from use of this expression in Eq. (j23p . 

^,onds ^J2 E ^''''" Ao 4, ,c 4„ ,c bi,,, ,.M,,, + (h.c.) . (29) 

i = l j' ,j"lJ-const] 

The c effective lattice is identical to the original lattice-. For the particular set of rotated-electron pairs and 
corresponding c fermion pairs of real-space coordinates r}' and r*," associated with the same point of coordinate 
^j ~ R' + ^i"]/2 ^ ^di the phases 6j'j" of Eq. (P^ are nearly equal. The reason is that such a set of pairs interact 
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with the two-site and two-spinon bonds of the same si fermion of real-space coordinate fj whose operator is given in 
Eq. (p6)) . Thus we introduce the fohowing notation for the phase factors, 

The phases Oj^ and 9j^i are identical to the phases 9j'j"fi and Oj'j",i, respectively, of the primary c fermion pairs 
associated with the four primary two-site bonds referring to nearest-neighboring sites of the in-plane spin effective 
lattice whose centre of mass is near the real-space point of coordinate "Tj = [fj> + fji']/2 — r^i- 

In order to describe the physics of the quantum liquid in terms of c and si fermions, in Appendix C it is shown 
that the Hamiltonian terms (1291) are approximately equivalent to, 

^e/r = E E ^l^ofl, 4, 4,, + (h.C.) . (31) 

Here the summation X](i' j") runs over c fermion pairs associated with rotated-electron pairs on singly occupied 
sites. In the spin effective lattice the centre [rji -|-r\,//]/2 of the corresponding in-plane nearest-neighboring real-space 
coordinates rji and ry of the spins of such rotated electrons is near the real-space coordinate f\, = [fji + fjii]/2 — 
I [as/2] Cxj of the si fermion. In the present N^ — > oo thermodynamic limit it coincides with it. Moreover, 



Here 4>j'.j,si are the phases defined in Eq. ([2]). It follows from the derivation of Appendix C that the primary two-site 
bonds, i.e. those whose length ^o == ^s = a/vl— 1e is minimum, have most of the spectral weight of a si bond 
particlc^-^ and corresponding si fermion. The coefficient l/|/io| of expression ([5T|) compensates the weight associated 
with bonds of larger length. 

The phases 9cp — 9cp{rj) appearing in the Hamiltonian terms (PT|) are associated with the rotated-electron pair 
and corresponding c fermion pair whose centre of mass [fji + fjir]/2 is approximately at the real-space coordinate 
^j ~ l^j' ~^ '^j"]/2 ~ I [as/2] ex^ of the si fermion of operator /I ^^. Such Hamiltonian terms describe the interaction 
of the c fermion pair with the si fermion. That pair feels the latter object through such a phase. Hence the real- 
space coordinate Vj in the argument of the phases 6cp = ^cp(?^j) corresponds both approximately to the centre of 
mass of the c fermion pair and to the real-space coordinate of the si fermion that the two c fermions interact with. 
For the construction of an effective action for these important phases the contributions from the Hamiltonian terms 
pip involving nearest-neighboring sites at rji and fji> are sufficient. In turn, concerning some other aspects of the 
virtual-electron pairing mechanism the contributions from rotated-electron pairs and corresponding c fermion pairs 
at larger distances cannot be ignored. 

Since 9cp = [Oj — 0° ^i], the c fermion pairs of Eq. ([5T|) feel the si fermion effective vector potential ^ through the 

phases 4>^j^si- Those are obtained from the operator phase t/ij^si also given in that equation by replacing fl ^ sifr-,.si by 
its average (/I , fp., si) ~ 1. Indeed, for the si fermion occupancies of the states that span the one- and two-electron 
subspace there are none, one, or two unoccupied sites in the si effective lattice. In addition, the total number of such 
sites N^^_^ « A''/2 = (1 — x)N^/2 is of the order of N^. Therefore, within the present thermodynamic limit such a 
replacement is a good approximation^. Furthermore, for the states that span such a subspace the fluctuations of the 
phases (p^ ^i are very small. Hence the fluctuations of the phases 9cp of Eq. l\'S2\i are fully controlled by those of the 
phases 9jfi and 9j^i of Eq. (l30l) . 

The effective Hamiltonian terms (j3ip are consistent with considering only the Hamiltonian terms, 

N/2 

^" = E E ^V'[%.t%„,-%,,4,„.t] + (h-c.). (33) 

i=iO'j") 

Those correspond to the spin-singlet rotated-electron pairs of Eq. (j23p involving nearest-neighboring sites. Indeed, 
the spin-singlet rotated-electron pairs of Eq. ([33]) involve the same sites as the c fermion pairs of Eq. ([3T|) . 



2. The fluctuations of the phases 9j^o and 9 



j,i 



A next step is the construction of an effective action for the phases 9cp. The exact calculation of some of the physical 
quantities involved in that derivation is an involved open problem. Nevertheless, one can construct by means of several 
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approximations an effective action whose key features faithfully reflect general properties of the VEP quantum-liquid 
microscopic Hamiltonian. To follow the fluctuations in 6j,q and 9j^i, we should integrate the c and si fermions in 
the expression of Eq. ([31]) about a suitable saddle point. Such a procedure can as well be fulfilled by integrating the 
rotated-electrons in the expression of Eq. p3p . Those are however extremely complex problems. 

Fortunately, since the VEP quantum-liquid microscopic Hamiltonian has the same general form as the microscopic 
Hamiltonian (1) of Ref.— , the corresponding effective action for the phases 9cp has also basic similarities to that 
constructed in that reference. Given the nearly equivalence between the two actions, we omit here the details that are 
common to both investigations. Those can be found in Refsi^i^. The main results are summarized in the following. 
They are used and further developed in Section HI-D and following sections. Alike in the problem of Refsi^i^, 
one arrives to an effective continuum Lagrangian. In spite of being a simplification of the selected rotated-electron 
Hamiltonian terms, such a Lagrangian is expected to faithfully refiect the general properties of the microscopic problem 
under consideration. 

First one finds that our phases 6j,q and 6j^i correspond to the phases df,¥, and OL,,,, respectively, of Refs<^i^, 
where the site indices //' are called jk. Except that the rotated electrons play here the role plaid by the electrons 
in that reference, the physics is very similar. For instance, the charge— 2e sector of Ref<^ refers to the charge— 2e 
c fermion pairs. Indeed, upon expression of the rotated-electron operators in terms of c fermion and two-site and 
two-spinon bond operators one finds that the Cooper pairs of that reference are mapped onto the virtual-electron pairs 
considered below. The charge degrees of freedom of such virtual-electron pairs refer to the c fermion pairs. Within 
charge excitations the latter pairs behave independently of the si fermion that contains the two-spinon spin-singlet 
configuration of each virtual-electron pair. 

Second it is found that the fiuctuations in the phases 9j_o and 9j_i proliferate in different regions of the phase 
diagram. For finite hole concentrations in the range x £ {xq , x* ) and vanishing temperature T — the fiuctuations 
of the phases Oj^ increase for x — >■ Xc and remain large for finite hole concentrations in the range x e {xq, Xc). Those 
of Oj^i increase for x — >■ x*. It follows that the fiuctuations of the phases 9j increase both for x — >■ Xc and x — ?> x*. 
Otherwise, the zero-temperturc fiuctuations of the phases Oj^q and Oj^i remain small for finite hole concentrations 
in the ranges x S (xc,x*) and x G (xo,x*), respectively. Here x* is the critical hole concentration of Eq. (IA2[) of 
Appendix A introduced in Ref.^. Above it there is no short-range spin order at zero temperature. As further discussed 
below, the critical hole concentration Xc that emerges from such studies is directly related to the Ginzburg number— 
and critical hole concentration xq, Xc « Gi -|- xq. The expression of Gi suitable to the present quantum problem is 
provided below in Section HI-D. That relation of Xc to the Ginzburg number and xq is valid for very small values of 
the 3D uniaxial anisotropy parameter e^ = m*/M and 7^ « 1. For those Xc ~ 0.05 and we find below that Gi w 0.037 
for LSCO and Gi « 0.026 for the remaining four representative hole-doped cuprates. 

Third one finds that at zero temperature and hole concentrations in the range x S (xc,x*) and for temperatures 
T lower than a critical temperature Tc < T* and a smaller temperature-dependent hole concentration range centered 
at the optimal hole concentration x — Xop ~ (xc -I- x*)/2 the phases Oj and thus 9cp all line up into a phase-coherent 
superconductor. For these x and T ranges both the c fermion kinetic energy and the energy order parameter 2|A| 
of the short-range spin correlations are finite. Combination of the results on the fluctuations of the phases 9jfi and 
9j^i with those of Section HI- A on the general d-wave spectrum ([20)1 then provides strong evidence that for the VEP 
quantum liquid there is a d-wave long-range superconducting order for such hole concentration and temperature 
ranges. In that phase the fluctuations of both the phases dj^ and Oj^i are small. This ensures that the average phase 
factors (e*^^°) ^ and (e*^^-') ^ are non vanishing. It then follows that the average phase factor {e^^"^) 7^ is also 
non vanishing. This is consistent with the phases O^p all lining up into a phase- coherent superconductor. 

Fourth it is found that for finite hole concentrations in the range x £ (xq, Xc) at zero temperature and a temperature- 
dependent hole concentration range for temperatures belonging to the range T G {Tc,T*) where Tc = for x < Xc 
and Tc > for x G (xc,x*) and the pseudogap temperature T* is approximately given by Eq. (j22p the fiuctuations 
of the phases 9j^i are small yet those of the phases Oj^ are large. As a result, one finds that the average phase factor 
^e'^3.0) = vanishes and the average phase factor (e*^^-^) ^ remains non vanishing. This then implies that the 
average phase factor {e^^"^) = vanishes. In the corresponding pseudogap state there is no phase coherence and thus 
no long-range superconducting order. Nonetheless there remain strong pairing correlations. Indeed, virtual-electron 
pairs remain existing in such a normal state yet their phases 9cp do dot line up. Consistently, the virtual electron pair 
- VEP quantum liquid refers both to the normal and superconducting states. In addition the short-range spin order 
prevails and infinite vorticity loops proliferate. Since the fiuctuations of the phases Oj^i remain small, the monopole- 
antimonopole pairs of the type considered in Ref4 remain bound and the vorticity is conserved in the low-energy 
limit. 

Finally one finds that for the hole concentration range x G (0, xo) at zero temperature T = and a small temperature 
dependent hole concentration range for finite temperatures the fiuctuations of the phases 9jfl and Oj^i are large. It 
follows that the average phase factors (e'^^'°) — 0, (e**^'^) = 0, and (e*^=p) — all vanish. As a result vorticity is 
not properly conserved. Then the above monopole-antimonopole pairs unbind and proliferate. At T = this is a 
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Mott-Hubbard insulator quantum phase. The corresponding ground state has long-range antiferromagnetic order. 

Alike in the problem of Ref.^, a sharp distinction is drawn between the microscopic properties of the quantum liquids 
referring to hole concentrations in the range x £ (0,xo) and finite hole concentrations x G (xojX*), respectively, in 
terms of the global symmetry of the effective action for the phases 9j . For the latter hole concentrations it is a global 
C/(l) symmetry. It is commonly understood that Hamiltonian symmetries by themselves are not sufficient to prove 
that a particular symmetry is broken in the ground state. However, for vanishing spin density m = and finite 
hole concentrations in the range x G (a;o,a;*) the global symmetry of that action is related to the hidden global U{1) 
symmetry of the Hubbard model on the square lattice found recentlyii. The latter f/(l) symmetry is contained in the 
model global 50(3) x SO{i) x C/(l) = [5*0(4) x U{l)]/Z2 symmetry, beyond its previously well-known global 50(4) 
symmetry-!^. Indeed the representations of the model hidden global U{\) symmetry correspond to the occupancy 
configurations of the c fermions in their c momentum band^iii. 

Consistently, for the present quantum problem in the normal state at zero temperature and finite hole concentrations 
in the range x g {xt^^Xc) and in the superconducting state at zero temperature and hole concentrations in the range 
X € {xc,x^,) the c fermions configurations are organized in terms of zero-momentum c fermion pairs. Those play a 
key role in the non-coherent and coherent virtual-electron pairing, respectively, of these two states. In turn, at m = 
and a; = the c band is for the absolute ground state of the Hubbard model on a square lattice full so that there is a 
single occupancy configuration. Consistently, the symmetry of the effective action for the phases dj is instead a local 
compact gauge symmetry. Here that holds for x e (0,xo) as well. 

Our results are inconclusive on whether the ground state of the t±/t — square-lattice Hubbard model ([T]) is 
superconducting for intermediate U/At values and hole concentrations in the range x € (xc,a;*). If it is, the zero- 
temperature lower hole concentration Xc must have an expression different from that found here, Xc ~ Gi-l-xo. Indeed, 
this expression gives a;c — > oo for t±/t —^ and thus M/ml — > oo for [//4i > 0. However the expression Xc ^ Gi + xq 
applies provided that Xc remains small. It indicates though that the superconducting dome hole-concentration width 
(x* — Xc) decreases upon decreasing tj_/t in the regime where (xc — xq) = Gi — G [Al/m*] + xq remains small. 

For small hole concentrations obeying the inequality < {x — xo) ^ 1 the motion of the c fermion pairs and the 
associated kinetic energy play the role of symmetry-breaking Higgs terms. The presence of such terms supresses free 
monopoles and is behind the replacement at T = of the long-range antiferromagnetic order for x € (0,a:o) by the 
short-range incommensurate-spiral spin order for < (x — xq) <^ 1—. Indeed, independently of the occurrence of phase 
coherence virtual-electron pairs and corresponding c fermion pairs prevail for x > xq and only for x € (0,a;o) are 
replaced by the monopoles and antimonopoles. The Higgs terms associated with the motion of the c fermion pairs 
leaves behind the global U{1) symmetry of the short-range spin ordered phase. This is alike in the related problem 
of Refi^, where the motion of the c fermion pairs corresponds to that of the — 2e charges. 

D. Relation of virtual-electron pairing to spinon and c fermion pairings and the superconducting order 

parameter 

In this and following sections further evidence is provided that for intermediate (and probably large) f//4i values 
a ground-state long-range superconducting order occurs in the VEP quantum liquid at zero spin density m = and 
hole concentrations in the range x S (xc,x*), as a by-product of the short-range spin order. 

1. Small suppression effects and relation of virtual- electron pairing to spinon and c fermion pairings 

For the hole-concentration range x G {xc,x^,) at zero temperature and a smaller temperature dependent x range 
centered at x = Xop ~ {xc+x^,)/2 for finite temperatures below the critical temperature Tc given below the fluctuations 
of the phases 0jfi and 9j,i remain small. Hence the amplitudes, 

ff = l(e''^)l = l(e^''-°+'-^l)l=ffo5i; 5o = |(e'^-°)|; gi ^ \{e^'-'}\ , (34) 

remain finite. Such amplitudes play an important role in the physics of the VEP quantum liquid. 

According to our scheme, the effects of the cuprates intrinsic disorder or superfluid density anisotropy^^i^, called 
here suppression effects, are behind the lessening of the experimental critical temperature relative to its magnitude 
predicted by the 7d = 1 VEP quantum liquid scheme. Within our oversimplified description of such effects, they are 
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accounted for by a single suppression coefficient, 



7c/ Tc\aa=Q \ X 

- = (' - -"'•') ^ {"■ S) '« I - (»• l) ^ « <! - !• S) '- S <!' 5) ■ <'^' 

Here g denotes the amplitude g of Eq. ([Ml) at T = 0, whose maximum magnitude ^™''^ — 7c/4 is found below to 
be reached at a: = Xop — {xc + a;*)/2, Tc\aa=o denotes the critical temperature in the absence of suppression effects 
introduced in the following, and Tc is its actual magnitude as obtained from experiments. 

For g — ^ the very strong phase fluctuations are the main mechanism for depressing T^, so that 7^; — > 1. This is 
found below to occur both for < (a; — Xc) ^ 1 and < (x^ — a;) <C 1. In turn, 7^ reaches its minimum magnitude 
7™*" = (1 — ctd) at the optimal hole concentration x ss Xop — {xc + x* )/2, when such fluctuations are smallest. For the 
five representative systems Xc/x^ < 1/4 so that the strength of their suppression effects is measured by the magnitude 
of ad & (0, Xc/x^) and thus of 7™"' g (7^, 1). The a^ = (1 — 7™™) ranges given in Eq. ([55]) are justified below. 

The validity of our scheme requires that the experimental T™"^ magnitude is smaller than T™°-^\a^=o and larger 
than 7cT™''^|„,=o where 7^ « 0.81 and thus xjx* « 0.19 < 1/4 for x^ = 0.05 and x^ = 0.27. In Ref.^^ it is 
found that 7^™ e (0.94, 0.98) and thus ad S (0.02, 0.06) for the four representative systems other than LSCO. The 
corresponding weakness of the suppression effects is consistent with the experimental results of Refsj^i^. In turn, 
it is found that 77*" = 0.82 and thus aa ^ 0.18 for LSCO. The inequalities 7cT™°^|„^=o < T™"^ < T^-^^la^^o, 
^mm ^ ^^^ g^j^i^ Q,^ ^ Xc/xt for Xc/x* ~ 0.19 < 1/4 are then met for all representative systems. 

For temperatures below the pseudogap temperature T* the short-range spin order parameter 2|A| is identified with 
2|A| = gi 2Ao. For a large range of finite temperatures below T* such that [T* —T) > is not too small the amplitude 
gi = |(e'^^'^)| is expected to be approximately given by gi « gi where gi is its zero-temperature value. It follows 
that provided that (T* — T) > is not too small the energy scale 2|A| is for the range U/At E {uo,Ut^) for which 
X* e (0.23, 0.32) approximately given by, 

2|A| = 51 2Ao « ffi 2Ao « 2Ao (l - ^) ■ (36) 

The results of Ref.- about the dependence of the energy scale |A| on the hole concentration x then indicate that 
for the square-lattice quantum liquid the amplitude gi reads gi « (1 — x/x*) for small hole concentrations < x <C 1. 
This is expected to hold for the VEP quantum liquid for < (x — xq) ^ 1. Interestingly, this expression has for the 
whole range x S (xo,x*) the correct behavior gi < 1, which is imposed by the inequality gi = \{e^^^'^)\ < 1. Indeed 
its maximum magnitude reached for (x — xq) -^ reads ^i = 70 = (1 — xq/x*) < 1. Provided that U/At G (uo,U7r) 
and thus the expression 2|A| sa 2Ao(l — x/x^,) is valid^, we then consider that 51 w (1 — x/x*) for hole concentrations 
in the range x G (xo,x*). 

Concerning properties such as the fluctuations of the phases 9j^o and 0j^i, the Hamiltonian terms p3p involving 
spin-singlet rotated-electron pairs at nearest-neighboring sites contain the relevant information. In turn, in order to 
capture the extra information needed for the further study of the coherent virtual-electron pairing and clarification 
of its relation to c fermion pairing and si fermion spinon pairing one must take into account the more general 
Hamiltonian terms provided in Eq. (|23|) . When expressed in terms of c fermion operators and two-site and two- 
spinon bond operators these Hamiltonian terms are of the form given in Eq. (I29p . 

The amplitude fluctuations near a given real-space point can be neglected. Nevertheless the x dependence of the 
average over the whole system |(e j'^")! Aq of the quantity e^^'j" Aq appearing in the Hamiltonian terms provided 
in Eq. (l29l) plays an important role. It follows from Eqs. (p4)) and (|36|) that |(e*^j'3")| Aq can at zero temperature 
be written as go |A| = g Aq. Here go and g denote the T = values of the amplitudes go and g, respectively. The 
physical meaning of such terms is that the two-site and two-spinon bonds of the si fermion of Eq. (I26p provide 
through the residual interactions with the c fermion pairs the energy needed for effective pairing coupling between 
the c fermions. The source of such an energy transfer is the energy scale |A|, which is one half the order parameter 
associated with the short-range spin correlations. Hence within the VEP quantum-liquid scheme the short-range spin 
correlations provide the energy needed for the c fermion effective coupling involved in the pairing mechanism behind 
the long-range superconducting order. And this is consistent with the short-range spin order parameter 2|A| having 
the form 2|A| = gi 2Ao w gi 2Ao, as given in Eq. ((36l) . 

The expression in Eq. (23) of the VEP quantum-liquid effective microscopic Hamiltonian terms of Eq. (P^ in 
terms of the local c fermion operators and the two-site and two-spinon bond operators contains important physical 
information. Indeed it reveals that the set of c fermion pairs whose centre of mass [rj' + rj'i]/2 nearly coincides 
with the real-space coordinate rj = [rj> + fj>>]/2 — r^i w [rj> + fj>>]/2 of a given local si fermion interact with the 
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two-site and two-spinon bonds contained in that si fermion, as given in Eq. (pS)) . An important concept is then 
that of a local virtual-electron pair of real-space coordinate fj. Such an object has the same real-space coordinate 
as a given local si fermion. Indeed there is a local virtual-electron pair for each local si fermion. In addition 
to that si fermion, a local virtual-electron pair of real-space coordinate fj involves a superposition of quantum 
configurations of all c fermion pairs whose centre of mass [fji + fjii]/2 nearly coincides with its real-space coordinate 



Ti = ri' +ri 



/2 — f2i Pi [fji + fj"]/2. Each of its configurations includes the charge — 2e of a c fermion pair of centre 
of mass [fji + fj"]/2 and the spin-singlet configuration of two spin-1/2 spinous of the local si fermion of real-space 
coordinate fj — [fj' + fj"]/2 — f^^ « [fji + fjii]/2. In each of these configurations the local si fermion of real-space 
coordinate fj = [fji +fjii]/2 — f^i « [fji +fjii]/2 interacts with the c fermions of a pair of centre of mass [fji +fjii]/2. 
As discussed below, the c fermion effective coupling results from the energy supplied to each pair by the si fermion 
under consideration. 

In the pseudogap state whose x and T range was mentioned above the amplitude g vanishes. However, there remain 
c fermion pairing correlations associated with the amplitude gi, which remains finite. Indeed, the fluctuations of the 
phases 9i remain small. Therefore, the concepts of a virtual-electron pair and a VEP quantum liquid hold for the 
pseudogap state as well, yet for it such a pairing has no phase coherence. 

The c fermion discrete momenta qj — [q^ — f] where j = 1,...,N^ and si fermion discrete momenta qj where 
j = 1,...,A^^^^ are the conjugate of the real-space coordinates of the c and si effective lattices of such objects, 
respectively. Such momenta are good quantum numbers for the square-lattice quantum liquid of Ref.^. They are close 
to good quantum numbers for the VEP quantum liquid at hole concentrations in the range x € (xc, x*). Alike a local 
virtual-electron pair has the same spatial coordinate fj as the corresponding local si fermion, a virtual-electron pair 
configuration of momentum qj involves one si fermion carrying that momentum. Therefore, the virtual-electron pair 
configuration momenta belong to the si band. Such a virtual-electon pair configuration consists of a zero-momentum 
pair of c fermions of hole momenta <f'' and — (f'', which interact with a si fermion of momentum q. Within a local 
virtual-electron pair the well-defined set of c fermion pairs with the same centre of mass [fji + fjii]/2 interact with 
the si fermion of real-space coordinate fj = [fji + fjii]/2 — f^i ^ [^j' + '^j"]/2- Consistently, the effective pairing 
coupling between c fermions of hole momenta q^ and — g'' results from elementary processes where such objects have 
interactions with a si fermion of well-defined momentum q. 

We emphasize that a pair of c fermions of hole momenta q^ and — g'' is a superposition of local c fermion pairs 
whose centre of mass [fji -\-fjii\/2 is different and thus interact with different si fermions. Therefore, in the occupancy 
configurations of the superconducting ground state there is no one-to-one correspondence between a coherent c fermion 
pair and a given si fermion. The occupancies of such a ground state involve superpositions of phase-coherent virtual- 
electron pair configurations of different momentum. The same pair of c fermions of hole momenta q^ and —q^ 
participates in different virtual-electron pair configurations where it interacts with a si fermion of different momentum 
q. And vice versa, each si fermion of momentum g participates in different virtual-electron pair configurations where 
it interacts with different pairs of c fermions. 

The short-range spin correlations provide through the c - si fermion interactions the energy needed for the effective 
coupling between the two c fermions of a pair. These zero-momentum c fermion pairs couple to charge probes 
independently of the si fermions. The superconducting macroscopic condensate refers to such pairs. In turn, the 
virtual-electron pair configurations are coupled to spin probes through the si fermions. It follows that concerning 
charge excitations the zero-momentum coherent c fermion pairs behave as independent objects relative to the virtual- 
electron pair configurations. In turn, one-electron and spin excitations break phase-coherent virtual-electron pairs. 
Indeed, as discussed below virtual-electron pairs exist in intermediate virtual states generated by pair breaking 
processes upon such excitations. A virtual electron pair configuration involves the charge — 2e of its c fermion pair 
and the spin-singlet configuration of the two spin-1/2 spinous of its composite si fermion. 

A virtual-electron pairing configuration involves two types of pairing: (i) The zero-momentum c fermion pairing 
whose effective coupling between the c fermions results from residual interactions with the si fermion; (ii) The spin- 
singlet spinon pairing of the composite two-spinon si fermion. As reported below in Section III-E, at zero temperature 
and hole concentrations in the range x £ (xc,a;*) one may reach the pseudogap state by applying a magnetic field 
H aligned perpendicular to the square-lattice planes. The occupancies of the gi > Q and g — normal ground 
state involve again superpositions of virtual-electron pair configurations of different momentum q. However such pair 
configurations and corresponding c fermion pairs have no phase coherence. The same applies to the normal ground 
state for vanishing magnetic field H = and finite hole concentrations in the range x £ {xo,Xc)- 

Strong effective coupling between c fermions of a pair is defined below as that whose breaking under one-electron 
excitations leads to sharp spectral features in the corresponding (fc, u;)-space distribution. One then classifies the c - si 
interactions into two classes: Those that lead and do not lead to strong effective coupling, respectively. Importantly, 
only c fermion pairs with strong effective coupling can participate in phase-coherent virtual electron pair configurations. 
Indeed, c fermion strong effective coupling is a necessary but not sufficient condition for coherent pairing. 
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According to the analysis presented below in Section IV-C, the momentum g of strongly coupled virtual-electron 
pairs may belong to different si — sc lines labelled by the line normal momentum absolute value g^^, S {q^ciIbsi)- 
Here sc stands for strong coupling. The si — sc line of the momentum of a si fermion contributing to the strong 
effective coupling of two c fermions is determined by the absolute value q^ of their hole momenta g'' and —q^. 
Indeed, there is an one-to-one correspondence between each g^^ G (g^, q^si) value labeling a given si — sc line and 
qh _ |^/i| g {q%i.,q^c)- The set of c fermions with the same hole momentum absolute value q^ = \q^\ G {qpcjlec) s-^^ 
involved in virtual-electron pairs whose momentum g belongs to the same si — sc line. Hence the si fermions whose 
momentum belongs to a given si — sc line interact with the c fermions whose hole momenta q^ and —q^ belong to 
an approximately circular c — sc line of radius g'' £ {Qfc' lee) centered at — tt = — [tt, tt]. And vice versa, g^ and g^^ 
are the minimum momentum and corresponding maximum hole momentum, respectively, for which there is c fermion 
strong effective coupling. (Such properties refer only to c - si fermion interactions leading to c fermion strong effective 
coupling.) 

Combination of the expression of the number of holes in the si band N^^ given in Eq. (jl9|) with that provided 
in Eq. (PT|) of its deviation 6Nsi due to a ground-state - excited-state transition gives N^^ = [SN — 26Nsi — '2Ns2]- 
Here Ns2 = 0, 1 is the number of four-spinon s2 fermions^. This is the form that the exact relation between the 
excited-states number of holes in the si band and the corresponding deviation 5N from the ground-state electron- 
number value has for the one- and two-electron subspace. It implies that the matrix elements between one-electron 
(and two-electron) excitations and excited states with an even (and an odd) number of si fermion holes vanish. Such 
an exact selection rule applies when the initial ground state refers to vanishing spin density m — 0. It reveals that 
excited states with a single hole in the si band are generated by application onto that ground state of one-electron 
operators. In turn, those with none or two holes in that band are generated by application of two-electron operators. 
Furthermore, the transformation laws under the electron - rotated-electron transformation of the quantum objects 
whose occupancy configurations generate the energy eigenstates of the Hubbard model can be used to show that nearly 
the whole one-electron (and two-electron) spectral weight corresponds to excited states with a single si fermion hole 
(and none or two si fermion holes)^. The same is expected to hold for the VEP quantum liquid. 

In conventional superconductivity the objects that pair are Fermi-liquid quasiparticles^Sil^. Here virtual-electron 
pairing involves a pair of charge c fermions and one spin- neutral two-spinon si fermion. The occurrence of si fermion 
spin-singlet spinon pairing is needed for coherent virtual-electron pairing. In contrast, the former pairing can occur 
independently of the latter. Indeed, the si fermion spin-singlet spinon pairing occurs both in the superconducting and 
pseudogap states. Consistently, coherent virtual-electron pairing is a by-product of the si fermion spin-singlet spinon 
pairing. It occurs provided that the amplitude g — |(e'^^)| is finite, owing to a macroscopic quantum phase-coherent 
virtual-electron pair state. 

The coherent virtual-electron pairing involves phase coherence. The si fermion spinon pairing does not. Concerning 
charge excitations the coherent c fermion pairs of charge — 2e behave as independent entities, which carry the superfluid 
current. The flux quantization is ~hc/2e. Hence in the superconducting state in the presence of a vector potential, the 
c band dispersion whose expression is provided in Eqs. (jA10[) and (JA11[) of Appendix A is shifted and approximately 
given by. 
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Here Ocp is the phase of Eq. ([5^ and jdq'^) the c fermion current of Eq. (JA16I) of Appendix A. It reads jdq'^) ~ 
qp^/rnP for hole momenta q^ at or near the c Fermi line. Its absolute value is for x G (xcX*) and U/At S (uoi^ir) 
approximately given by Jc{q^) ~ y^xTT2/mP. The renormalized transport mass toJ? is provided in Eq. (JA16I) of 
Appendix A. In the units used here it reads trP = rc/2t. 

2. Virtual- electron pairs in virtual intermediate states of one- and two-electron excitations 

As discussed in Section III-A, removal of two electrons in a spin-singlet configuration leads to an excited state with 
no holes in the si band. Upon such an excitation, before the two involved c fermions recombine with one si fermion 
giving rise to the removed electron pair under consideration, the system is driven into an intermediate virtual state. 
In it the c fermion pair under consideration stops interacting with all si fermions other than that si fermion. Hence, 
in that intermediate virtual state both the c fermion pair and the si fermion participate in a single virtual-electron 
pair. In it such a pair exists as an individual object. The occurrence of that intermediate virtual state is behind the 
designation virtual-electron pair. This holds both for the superconducting and pseudogap states. Such a two-electron 
process corresponds to a gapless excitation branch. It refers to the last excitation energy 5E = given in Eq. ((20|) . 

In turn, under one-electron removal and spin excitations there emerges one and two holes in the si band, respectively. 
Also within such excitations the system is driven into an intermediate virtual-electron state. In it one c fermion pair 
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and one si fermion contribute to a single virtual electron pair: that broken under the excitation. Such a discontinuous 
change is brought about by application onto the ground state of the corresponding one-electron and spin two-electron 
operator, respectively. Within these processes the c fermion pair ceases to interact with si fermions other than that 
under consideration. Furthermore, that si fermion ceases to interact with all remaining c fermion pairs. The virtual 
electron pair and the corresponding c fermion pair and si fermion spinon pair are then broken under the one-electron 
removal or spin excitation. Finally, under the former excitation one c fermion and one spinon are removed along with 
the electron and one hole emerges in the si band. In turn, under a spin-triplet excitation a spinon spin-flip occurs. It 
generates two independent spinous. Under a spin-singlet excitation the two spinous recombine with those of another 
si fermion giving rise to a spin-singlet four-spinon s2 fermiorJ-. Such spin two-electron processes are followed by 
the emergence of two holes in the si band. In addition, the c fermion pair left over under these spin excitations 
restarts interacting within virtual-electron pairs with other si fermions. This is how one-electron and spin quantum 
measurements destroy the ground-state wavefunction and create a new wavefunction. 

3. Critical temperature, superconducting order parameter, and coherent length 

The order parameter associated with the short-range spin correlations is the maximum magnitude of the si fermion 
spinon-pairing energy 2|A|. It is expected that the absolute value 2|51| of the superconducting order parameter is the 
maximum magnitude of the coherent virtual-electron pairing energy. According to the pairing mechanism emerging 
from the above analysis, the short-range spin correlations supply the latter pairing energy. It is then expected that 
for hole concentrations in the range x G (xc,a;*) and zero temperature both the inequality 2|J7| < 2|A| and the 
relation 2|51| ex 2|A| hold. The pseudogap temperature ([^^ is given by T* w giAo/fc^. It is controlled by the 
zero-temperature amplitude gi, which is finite for hole concentrations in the range x € (a;o,x*) where there is short- 
range spin order. Within the present approach the superconductivity critical temperature reads Tc ~ j^g AQ/2kB, 
so that Tc\aa=o ~ gAo/2kB in Eq. ([55)1 . It is controlled by the amplitude g, which is finite for hole concentrations 
in the range x e {xc,x*). For it a long-range superconducting order coexists with the short-range spin order, as its 
by-product. 

Physically, the energy scales (i) 2|r2||T=o and (ii) AksTc are the maximum magnitude of the coherent c fermion 
pairing energy within pair breaking (i) at zero-temperature under spin-triplet excitations and (ii) upon increasing 
the temperature. The relation 2|A||j-=o ~ 2fcB T* holds. The corresponding expected relation 2|J7||r=o ~ ^ksT^. 
holds both in the limits < (a; — Xc) <C 1 and < (x* — x) <C 1. Hence in these limits when the phase fluctuations 
are strong the superconducting energy scales 2|J7||t=o and AksTc have the same magnitude. In turn, at cc = Xop = 
{xc + a;*)/2 such fluctuations are smallest at zero temperature and these two energy scales have slightly different 
magnitudes. Due to correlations the energy scales considered in the following are not additive, so that our discussion 
is merely qualitative. The maximum energy that the short-range spin correlations associated with the energy scale 
2|A|(xop) can supply to coherent pairing is at x = Xop approximately given by « 2|A|(a;op)|T=o/[27c] = Ao/2. 
In the absence of suppression effects the above pairing energies read 4:kBTc{xop) ~ 50 2|A|(xop)|t=o = 7c Ao/2 and 
2|J7|(a;op)|T=o ~ ^o/2. For pair breaking under spin excitations the pairing energy 2|51|(xop)|t=o equals that delivered 
by the short-range spin correlations « Ao/2. In contrast, in the case of Aks Tc{Xop) ~ (1 — Xc/x^,) Ao/2 only a fraction 
of the latter energy is supplied to pairing, the energy sa [xc/x^,) Ao/2 left over being lost due to phase and thermal 
fluctuations. If one accounts for the suppression effects the above two related pairing energies then approximately 
read ^kBT,{xop) « 7r"5o 2|A|(xop)|t=o = 7^™ 7c Ao/2 and 2\^\{xop)\t=o ~ 7^™ Ao/2. Hence the relation, 

^ksT^ixop) _ -}_u ^^) f38) 

holds. The following expression satisfles such limiting behaviors, 

I 11^ ^^ /3c V x*r,"-v ^ ' 

The inequalities ad < Xc/x* and Xc/x* < 1/2 in the range of the parameter a^ = (1 — 7™*") of Eq. ([55]) are justifled 
by the physical range and inequality 2|ri||y°Q e (g2Ao, Ao/2) and 2|A||t=o/[27c] = Ao/2 < 2|A||t=o, respectively, 
corresponding to x = Xgp. The latter follows from the physical requirement that the maximum energy supplied by 
the short-range spin correlations being smaller than their own energy scale 2|A|. In turn, concerning the above range 
only a fraction (1 — Q;rf)Ao/2 of the available maximum energy Ao/2 is supplied, whereas the energy a^ Ao/2 left 
over is lost due to the suppression effects. Our oversimplified scheme in terms of such effects is valid provided that 
2\Vl\\'^^1 = (1 - arf)Ao/2 G (52A0, Ao/2). At ad = x^/x* the equality 7^ = /3c holds, so that 2|J7||5?°g = 52A0. 
Otherwise, 7^ < $c and g2Ao < 2|r2||™°g < Ao/2 for ad < Xc/x* and 'y™™ > -y^. Below it is confirmed that 
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rj^max _ Xc(xop) and 2|r2||™°Q = 2|r2|(xop)|T=o- For the critical hole concentration values Xc ~ 0.05 and x* « 0.27 of 
the five representative systems the optimal hole concentration of Eq. (j38p is given by Xop ~ 0.16, in agreement with 
experiments. 

It follows from the above analysis that the order parameter 2f2 and the energy scale 4fcs Tc read, 

2n = e^'^-2\n\; 9,p ^ 6^ - ^l,,; ^, = ^,- o + ^j,i ; 4^^ T, « 7^ 5 2 Aq . (40) 

Here 2|f2| is for T — provided in Eq. ([39]), 9cp — Ocp{rj) are the phases of Eq. ([32]) . and g = gofji is the 
zero-temperature value of the overall amplitude given in Eq. (1341) . 

Generalization to finite temperature T > of the zero-temperature 2|ri| expression provided in Eq. (|39p leads for 
a temperature-dependent x range centered at a; = Xop to, 

2N = ^^; P.^fi-'-^'-^). (41) 

Here 7c is the parameter given in Eq. (PSJ) . Note that (3c\t=o = $c where Pc is defined in Eq. ([55]) . 

The 2fl expression of Eq. (HUl) is such that 2n ex e^^"" |(e*^'=p)| 2Ao « 6*^'='' 52A0. However, the amplitude fiuctua- 
tions are not accounted for twice. The point is that |(e*^'=p)| refers to an average over the whole system. Hence the 
amplitude |(e*^'^p)| is independent of the spatial coordinate r*,-. In turn, the fj dependence of the phase factor e^^"'' 
corresponds to a small real-space region around fj where the phases 9cp change little and smoothly. Then amplitude 
fluctuations can be neglected in that small region around fj and averaging e*^'=p over the local virtual-electron pairs 
contained in it gives |(e'^=i')„ear-j| ~ 1- This justifies why such fiuctuations are not accounted for twice. 

That local normalization is not fulfilled when the phase fiuctuations become large. Nevertheless, then 5 — >■ and 
thus 2U — >■ 0. Hence the expressions provided in Eqs. (|10)) and (PT|) remain valid. One then concludes that the range 
of validity of the expression 2J7 = e*^'^p2|ri| whose fj dependence occurs through the phases 9cp = [9j — (j>^ ^i] refers to 
a small real-space region around fj . We recall that the fluctuations of the phases 9cp result mostly from those of the 
phases 9j = [9jfi + 9j,i]. Indeed, the fiuctuations of the phases (fi^ ^i are very small for the whole hole-concentration 
range x G {xo,x^,). 

At ad = and thus 7d = 1 the superconducting energy scale Aks Tc of Eq. (|40p equals the zero-temperature 
average over the whole system |(e*^j'j")| Ao of the basic quantity c'^j'j" Aq appearing in the Hamiltonian terms (|29p . 
The expressions provided in Eqs. (IMl) . (15^1) . (PO]) . and (PT|) confirm that the short-range spin order parameter 2|A|, 
energy scale Aks Tc, and superconducting order parameter 2il are closely related. This is consistent with the short- 
range spin order and long-range superconducting order being closely related as well. The presence of the amplitude 
gi = |(e*^J'^)| within the overall amplitude g = go gi of the superconducting order parameter (|4ip is consistent with the 
short-range spin correlations providing the energy needed for the coherent virtual-electron pairing. Simultaneously, 
such a supplying of energy by the short-range spin correlations suppresses them through the increasing fluctuations 
of the phases 9j^i. The latter phases refer to the internal degrees of freedom of the virtual-electron pairs. 

The phases 9jfi and 0_,_i are such that go = |(e*^J")|T=o — and 51 — |(e*^J'i)|T=o — Ic for < (x — Xc) — > 
and go = \{e'^^^-°)\T=o = 1 and 51 = |(e'^^'i)|T=o = for < (a;* — x) — )■ 0. These behaviors are consistent with the 
critical temperature Tc being small and given by Tc ~ go Ao/2kB and Tc « gi Ao/2kB for < (a; — Xc) ^ 1 and 
< (x* — x) ^ 1, respectively. In these limits g — gogi ^ 0, so that 7^ — >■ 1, the suppression effects play no role, 
and the physics is controlled by the very strong phase fluctuations. Symmetry arguments associated with the physics 
specific to these limits for which the critical temperature Tc is controlled only by strong phase fluctuations imply that 
go ~ {x/x^y and g\ « ([a;* — a:]/a:*)^'' are for < (a: — Xc) ^ 1 and < (a:* — x) ^ 1, respectively, controlled by the 
same dynamical exponent z — 1 and unknown exponent v. 

Treatments involving the use of the effective action for the phases 9j without incorporating the Berry phased lead 
to I' « 2/3. More detailed treatments, incorporating the latter phase, lead to different values for that exponent^-. 
Symmetry arguments suggest that 51 « ([x* — a;]/a;*)^'' = (1 — x/x^.y has the same overall exponent zv for hole 
concentrations obeying both the inequalities < (a;* — a;) <C 1 and < (a; — a;o) <C 1. On combining the expression 
2A = gi 2Ao of Eq. ([55)1 with the square-lattice quantum-liquid results of Ref.^ for < a; ^ 1, which in our case 
hold for < (a; — Xc) ^ 1, we then find z^ = 1 so that v ~ 1. Within our scheme the suppression effects lessen the 
critical temperature and related energy scales, yet do not affect the phases and corresponding amplitude g — gogi- 
We then consistently consider that for approximately the range C//4i S (uoj^tt) and hole concentrations x G {xo,x^) 
the zero-temperature amplitude gi is given by 51 ss (1 — x/x^,). Furthermore, the expression go ~ {x — Xc)/{xsf — Xc) 
obeys for x e {xc,x^) the inequality go = |(e*^J'')|T=o < 1- It also obeys the expected boundary condition go — > 1 as 
X — >■ a:*. Therefore, we consider that go ^ (x — Xc)/(x* — Xc) for hole concentrations in the range x G {xc, x*). Hence 
within the present approach the expressions of the amplitudes go and gi of Eq. (p4)) are for intermediate interaction 
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values U /At E {uq, u-,^) and zero temperature approximately given by, 

g = |(e^'^)|T=o=5offi«/^^ f 1 " -) ^ ffo = |(e'''-'')|T=o « /^^ for x e (x„:z:.) , 

91 = |(e'^^'OlT=o« (l- — ) forxG (a;o,a;*); max{gi} = 70 = ( 1 - — ) for < (x - 2:0) < 1 , 

5o = g = for < a; < Xc and x G (a;*, 1) ; go = .9i = .9 = for a; e (0, a;o) and x G (a:*, 1) . (42) 

Due to the hole-trapping effects reported in Appendix B, the fluctuations of the phases Oj^i are very strong for 
X e (0,a;o) and the zero-temperature amplitude gi — |(e* J'^)|t=o vanishes. For x G (a;o,a;c) the hole-trapping 
effects remain active yet are weaker than for x € (0,a:o)- In contrast to the range a: S (0,a;o), the short-range 
spiral incommensurate spin order of the related square-lattice quantum liquid of Ref ji' survives for hole concentrations 
X € (xo^Xc), coexisting with the Anderson insulating behavior brought about by the hole-trapping effects. 

The behaviors of the zero-temperature amplitudes go and gi follow from the fluctuations of the corresponding 
phases 9o and 611 of Eqs. ([50)) and ([5^ becoming large for < (a: — Xc) ^ 1 and < (a;* — a;) <C 1, respectively. The 
singular behavior gi -^ 70 for < (a; — a:o) ^ and gi = at (a; — a;o) = is due to a sharp quantum phase transition. 
It marks the onset of the long-range antiferromagnetic order for x < xq. In turn, the singular behavior go — >■ 1 for 
< (a;* — a;) — > and go = at (a;* — a;) = is also due to a sharp quantum phase transition marking the onset to a 
disordered state without short-range spin order for x > x*. Due to such sharp quantum-phase transitions, the phases 
61 and ^0 have also a singular behavior at a; = xo and x = x* , respectively. The fluctuations of these phases are small 
for < (x — xo) ^ 1 and < (x* — x) <C 1 and large for x < xq and x > x*, respectively. 

Both in the Anderson insulator phase with short-range spiral incommensurate spin order occurring in the hole 
concentration range x £ (xo,Xc) due to the hole-trapping effects and in the phase corresponding to the range x S 
(xc,x*) for which the short-range spin order coexists with a long-range superconducting order the symmetry of the 
phases 6* = ^o + ^1 action is a global U{1). Such a global U{1) symmetry is related to the Hubbard-model global 
[5*0(4) X U{l)]/Z2 symmetry found recently beyond its known 50(4) symmetrjsii. 

On combining Eqs. (|40p and (HH one finds the following expression for the critical temperature, 

rr ^ rr., ^f, 4.9 \ Aq ^.Afl (x^xj/ x\Ao 

The related energy scale 2|0||7'=o is then given by expression ([M]) with T^ as provided here. Hence the energy scales 
2|A||t=o of Eq. ([5S)) at zero temperature and 2|r2||T=o associated with the short-range spin order and long-range 
superconducting order read, 

oiAii " oA fl,' \ oioii ^^fcsTc 7dg2Ao , , 

2|A||t=o =5i2Ao6'(x-xo); 2|S2||t=o « — ^ — = ^ , (44) 

Pc Pc 

respectively. Analysis of these expressions reveals that the physical requirement that 2|51| < 2|A| is met for the whole 
hole concentration range x € (xc,x*) provided that the a^ = (1 — 7™™) ranges given in Eq. (1551) are fulfilled. 

It follows from the expressions of the energy scale 2|ri||T=o and critical temperature Tc provided in Eqs. ([5^ and 
(|43p . respectively, that the maximum magnitudes of these quantities read, 

|ii|lT=o = 7d ~Y ' " ^^'^ ^"8^' ^'^ = (l-"d) G (7c,l)- (45) 



Those are achieved at the optimal hole concentration given in Eq. 

Complementarily to the complex parameter 2il of Eq. PO)) . one can consider the order parameter (j)cp — 
•y/n^Jp/2 e'^'=p . Here Ucp denotes the density of paired c fermions contributing to phase-coherent virtual-electron pair 
configurations. The x dependence of Ucp is studied below in Section IV. The complex parameter (pcpir) = yncpj^e"^^"'' 
describes the macroscopic properties of the zero-momentum c fermion pairs superfluid condensate. Indeed, |<^cpP = 
ncp/2 is a measure of the local superfiuid density of coherent c fermion pairs. The average c fermion distance or length 
^1 of the set of c fermion pairs contributing to a local virtual-electron pair is related to the coherent length ^ associ- 
ated with coherent virtual-electron pairing. For the isotropic Fermi- velocity range x G (xci,Xc2) and approximately 
U/At G {uo,UTr) such lengths read, 
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In the expression of ^i we included explicitly here the lattice spacing a and the factor ^d^g appearing in the ^ 
expression assures that its minimum magnitude ^™™ equals ^i. It is achieved at zero temperature and x — Xop- 
Furthermore, ^g is the distance between the c fermions of each pair contributing to the same local virtual-electron 
pair. For the present isotropic Fermi-velocity range the Fermi velocity Vf equals approximately the c fermion velocity 
Vpc ~ -y/icTr 2/?7i* of Eq. ^. The quantity \2hg\ — y/2h*2hg is the absolute value of the coefficients controlling 

the ratios {fl,„^jl„„Jg/\{fl,,jl„Jo\ ~ [2h*g/\2ho\] of Eq. ^ of Appendix C. Those are associated with a c 
fermion pair with the same real-space coordinates fj' and fj" as a two-site bond of a given rotated-electron pair. The 
relation between such c fermion and rotated-electron pair is given in Eq. (|28p. The coefficients 2hg obey the sum-rule 
provided in Eq. (jC3[) of Appendix C. They are twice those appearing in the expression of the annihilation operator 
of the local si fermion of a local virtual-electron pair. Such an expression is obtained from Eq. ((26|) . The absolute 
value \2hg\ = y/2h*2hg decreases upon increasing the magnitude of the length ^g^^- \2hg\ is largest a.t g — 0, which 

corresponds to ^o ~ Os = a/^/T~~x. It falls rapidly upon increasing g. Therefore, the average length ^i is typically 
very small as compared to the penetration depth considered below in Section IV-E. 

The overall amplitude g — go9i associated with the phases 6j = Oj^ + 9j^i whose dome-like x dependence is for 
zero temperature given in Eq. (|42|) gives a measure of the strength of the c fermion effective coupling leading to 
the superconducting state coherent pairing. Such an effective coupling also occurs when gi > and go = yet for 
the corresponding pseudogap state does not lead to coherent virtual-electron pairing. In either state it is due to 
the energy provided to the c fermions of such pairs by the short-range spin correlations. This occurs through the 
residual interactions of the two c fermions with the si fermion within each virtual-electron pair. The energy lost by 
the spin-subsystem to supply the c fermions of a pair with the energy needed for their effective coupling weakens the 
short-range spin order. This occurs through the enhancement of the fluctuations of the phases 6j^i associated with the 
amplitude gi. For the c fermions gi refers to their effective pair coupling. Complementarily, for the si fermions it is 
behind the weakening of the short-range spin correlations upon increasing x. Consistently, the value of gi provides a 
measure of the ability and power of the short-range spin correlations supplying the c fermions with the energy needed 
for the occurrence of virtual-electron pairing. 

The physical picture that emerges is that the coherent-pair superconducting order rather than competing with 
the short-range spin order is a by-product of it. It coexists with the latter order for the hole concentration range 
X € (xc, Xi.) at zero temperature and a smaller x range centered at a: = Xop for finite temperatures below Tc- Moreover, 
the type of d-wave long-range superconducting order considered here cannot occur without the simultaneous occurrence 
of short-range spin correlations. As the limit < (cc* — x) ^ 1 is reached at zero temperature, the short-range spin 
correlations use up all their energy. This is consistent with the behaviors gi —> 0, g = go 9i = 9i —^ 0, 2|A||t=o —^ 0, 
2|ri||T=o —> 0, and 2|r2||T=o/2|A||T=o/ -^ 1 occurring in that limit. Indeed, 7^ — >■ 1 as g ^- for < (x* — cc) <C 1 
when the physics is controlled by strong phase fluctuations. The order parameter of the long-range superconducting 
order vanishes upon the vanishing of that of the short-range spin order. Such parameters vanish upon the simultaneous 
disappearance of the corresponding orders. 

Except in the limits < (x — Xc) <S^ 1 and < (x* — a;) <C 1 when the phase fluctuations are very strong, for 
X G {xc,x*) the suppression effects lessen both the critical temperature Tc of Eq. pS]) and superconducting energy 
scale 2|i7|| of Eq. (pi)) and enhance the coherent length ^ of Eq. (|46p . relative to their a^ = magnitudes. In turn, the 
c fermion energy dispersion £c{q'^) given in Eqs. (lAlOp and (JA11[) of Appendix A, the form of the si fermion energy 
dispersion esi{q) given below, and the magnitudes of the pseudogap temperature T* of Eq. ((22|) and the critical hole 
concentration x* of Eq. (|47p remain unaltered. Although here we consider intrinsic disorder, this is consistent with 
the experimental results of Refs4^i^, according to which the disorder induced by Zn substitution does not lead to a 
clear change in the magnitude of the pseudogap temperature T* . 



4- The ratios controlling the effects weak 3D uniaxial anisotropy and electronic correlations 

In the present limit of very weak 3D uniaxial anisotropy its effects occur mainly through the dependence of the 
critical hole concentration Xc on the parameter e^ . It follows that the amplitudes g and go of Eq. (|^^ , whose expres- 
sions involve Xc, also depend on the latter parameter. For approximately U/At > uo the critical hole concentration 
deviation (xc — xo) and the critical hole concentration x* are proportional to two important ratios. 
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U/At > uo . (47) 



That {xc — Xo) equals approximately the Ginzburg number Gi is a result obtained from the study of the effective 
action of the phases 9j. In the G expression provided here Xab denotes the in-plane penetration depth introduced 
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below in Section IV-E, T|7"^=^ is the maximum critical temperature of Eq. (H5t reached at a: = Xop^ and ^'"™ is the 
corresponding x — Xop minimum magnitude of the coherent length given in Eq. (j46l) . The ratio e^ = to*/M appearing 
in the {xc — xq) expression of Eq. (|47|) controls the effects of weak 3D uniaxial anisotropy. Here the U/At dependent 
mass m* is that of the c fermion energy dispersion of Eq. (|All|) of Appendix A and M ^ m* is the effective mass of 
Eq. ([5]) associated with electron hopping between the planes. 

The critical hole concentration x* expression given in Eq. (j47p is derived in Ref.^. The spin ratio Vg appearing 
in that expression and the a;* limiting values are provided in Eqs. (JA1[) and (jA2[) of Appendix A, respectively. For 
the smaller range of intermediate interaction values U/4t e {uq,ui) one approximately has that r^ = m'^/m* ~ 2rs. 
Hence for such a range the upper critical hole concentration can be expressed as x■^, « Tc/tt, as found in Ref.-'^. The 
ratios r^ = Ao/4VF°^ and r,, = m'^ /m'^ control the effects of electronic correlations. Here m^ = lirau/^^i^^ m* = l/2i 
and Aq and 4:Wgi are the energy scales of Eq. © . 

Our scheme is valid for very small values of the 3D uniaxial anisotropy mass ratio m*/M. That the relation 
{xc — Xq) « Gi holds is fully consistent with the related results of Ref.*. The studies of Ref.^^ reveal that a critical 
hole concentration Xc ~ 0.05 corresponds to e^ = m*^/M w 10^'' for LSCO and e^ — ml/AI w lO^'^ — 10^^ for 
the other four systems. This follows from G being very small and given by G w 10^^ and G w 10^* for LSCO 
and the other four systems, respectively. The smallness of G allows that Gi w {xc — xq) w 10^'^ is small in spite of 
1/e^ = M/ml being large and {xc — Xq) oc 1/e^ = M/m*^. Hence the value Xc ~ 0.05 is set by choosing an appropriate 
small magnitude for the 3D uniaxial anisotropy parameter e^ = m*^/M . In turn, the value x* « 0.27 is set by choosing 
U/4i ~ 1.525. Often one considers the range U/it G (wq,"^) for which according to Eq. (|A2p of Appendix A the 
X* values belong to the domain a;* G {2e~^/7T, I/tt) « (0.23, 0.32). Some of our expressions refer to the smaller range 
U/4t € {uo,ui) of intermediate interaction values for which Vc ~ 2rs « ttx* « 2e~''*"''''^ and x* G (0.23,0.28). 

Our results reveal that for constant values of C//4i the superconducting-dome hole-concentration width (x* — Xc) 
decreases upon further decreasing the ratio mass m*/M. That could mean that such a width vanishes in the 2D limit 
m*^/M — >■ so that the ground state of the Hubbard model on the square lattice is not superconducting. However, 
note that our results do not apply to the 2D limit e^ = ?7i*/A/ — > 0. Then Gi becomes large, so that the expression 
Xc ~ (Gi -I- xo) obtained from the effective action of the phases 9j may not be valid. 

5. The rate equations for suppression o/ 2| A| and variation of 4fc_B Tc 

The dependence on the hole concentration x of the zero-temperature short-range spin order parameter 2|A||t=o of 
Eq. p6p and critical temperature Tc provided in Eq. (j43p is described by the two rate differential equations under 
suitable and physical boundary conditions given in the following. Those are valid for vanishing spin density m — Q 
and finite hole concentrations in the ranges x G (xo,x*) and x G (xc,x*), respectively. 

The rate equation for suppression of the order parameter of the short-range spin correlations 2|A|(x)|t=o upon 
increasing the value of x and its boundary condition are for the approximate range U/At G {uq^Ut^) and x G (xo,x*) 
given by, 

i92|A|(x)|t=o Ao / x^\ 

^ = -TT — 6'(x*-x); 2|A|(xo)|t=o =7o2Ao; 70 = I 1 I- (48) 

Here 9{z) is a theta function such that 0{z) = 1 for z > Q and 0{z) = for z < 0. In turn, the rate equation for 
variation of the superconducting energy scale 4^^ T^ upon increasing x and the corresponding boundary condition 
read, 

a4fcB7U^^^27rfj:^^p|A|(x)-2|A|(xop)]|T=oe(a;*-x); 4fcBT,(xe)-0. (49) 

ox 7c Ts 

That the source term of the rate equation (f49l) is given by [77(27^ — l)/7crs][2|A|(x) — 2|A|(xop)]|t=o is consistent 
with the superconducting order being a by-product of the short-range spin correlations. The source of the energy 
provided by such spin correlations to sustain the superconducting energy Aks Tc is the pseudogap energy 2|A| in that 
source term. Consistently, the latter energy is the order parameter associated with short-range spin correlations. For 
hole concentrations below and above the optimal hole concentration Xop of Eq. (H5t the short-range spin correlations 
can be considered strong and weak, respectively. Indeed, for x G (xc,Xop) and x G {xop,x^) the above source term is 
positive and negative, respectively. 
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E. Zero-temperature magnetic-field c fermion diamagnetic orbital coupling and si fermion Zeeman coupling 

It is useful for the study of the pseudogap state to consider a uniform magnetic field H aligned perpendicular to 
the square-lattice plane, to suppress superconductivity. The usual type II superconductor field Hd is very small for 
the VEP quantum liquid at intermediate 11/4,1 e [uq^u.,^) values, so that here we consider that H^ sa 0. Our study 
refers to such a U /At range. 

It is well known that for the type-II superconductors the Pauli field Hp for which the Zeeman splitting spin alignment 
starts to be profitable is given by Hp « fc^ Tc/ ^.39- Here /is is the Bohr magneton and g ~ 2. For conventional 
type-II superconductors such a field is found to lie above the Hc2 line. In that case the Pauli pair breaking will become 
at zero temperature an important issue at high magnetic fields H ^ Hp > Hc2- In contrast, we find in the following 
that here Hc2 > Hp. Important properties are then that: 

The diamagnetic orbital coupling of the magnetic field is to the charge c fermions. It occurs through the c 
band hole momenta, q'^ — eA, where A is the vector potential associated with the magnetic field H = W x A. Hence 
diamagnetic pair breaking refers here to the zero-momentum charge c fermion pairing. 

The Zeeman coupling of the magnetic field is to the spin-singlet two-spinon si fermions. Hence Pauli pair 
breaking refers here to the spin-singlet spinon pairing of the si fermions. 

Since virtual-electron pairing involves both zero-momentum charge c fermion pairing and spin-singlet spinon pairing 
of the si fermions, diamagnetic orbital coupling and Zeeman coupling are direct probes of the two types of pairing 
involved in the virtual-electron pairing. 

That the energy needed for coherent c fermion pairing is supplied by the short-range spin correlations associated 
with si fermion spinon pairing implies an interplay between Pauli spinon pair breaking and diamagnetic c fermion 
pair breaking. Below a critical field Hq the energy supplied to the c fermions by the spin correlations associated with 
si fermion spinon pairing is at zero temperature enough to prevent that diamagnetic c fermion pair breaking stops 
phase coherence. For H < Hp « k^ T^/ ^b9 there are no significant effects on the spin-singlet spinon pairing from 
Zeeman coupling to the si fermions and only diamagnetic orbital coupling to the c fermions plays an active role in 
the physics. In turn, for H > Hp the si fermion spinon pairing starts to be affected by the Zeeman coupling. This 
lessens the energy supplied by the spin sub-system to the c fermion strong effective coupling. As a result, the effects 
of diamagnetic c fermion pair breaking become stronger for H > Hp and the critical field Hq > Hp above which there 
is no c fermion pairing phase coherence is at zero temperature proportional to Hp. However, for H € (Hp^ Hq) the 
diamagnetic orbital coupling to the c fermions remains much stronger than the Zeeman coupling to the si fermions. 
Within our scheme the critical field Hq is that above which both the diamagnetic orbital coupling to the c fermions 
and the Zeeman coupling to the si fermions play an active role. It is approximately given by Hq « 2Hp. Since 
Hp « ks Tc/fisg where 5 ~ 2 one then finds that at zero temperature Hq « fc^ Tc/ ^b- 

Hence within our picture, at zero temperature the fluctuations of the phases 9j remain small and long-range 
superconducting order prevails provided that the uniform magnetic field H is below a critical field Hq w ksTc/ fiB- 
Above that critical field, diamagnetic c fermion pair breaking prevents phase coherence yet there remain c fermion 
pairing correlations. At zero temperature such a critical field expression approximately reads, 

Hq « kB TJ^lB ~ Id i^^^ (i^^\^^ a; e [x,, x,) ; i/o"'"" = iT^'lc ^ , x^x^p. (50) 

(a;* - Xc) V a;* / 2^_b 8^s 

Here the hole concentration Xop is given in Eq. (p8| . Interestingly, the expressions of the energy scale |2ri||T=o given 
in Eq. ([55]) . critical temperature Tc in Eq. P5|) . and magnetic field Hq in Eq. ((5(1)) have the dome- like x dependence 
observed in the five representative hole-doped cuprates^"^. Within our scheme such a type of x dependence follows 
directly from that of the overall zero-temperature amplitude g = go 5i of Eq. ()42|) . The x dependence of such an 
amplitude is fully controlled by the interplay of the x dependences of the fluctuations of the phases 9jfl and Oj^i such 
that 9j = Ojfl + Oj^i, as discussed above in Section III-D. 

At zero temperature, magnetic field H in the range H S {Hq, Hc2), and finite hole concentrations in a iJ-dependent 
range below a hole concentration x = Xc2 the phases 9jfi have large fluctuations. This refers to the pseudogap state 
for which both the diamagnetic orbital coupling to the c fermions and the Zeeman coupling to the si fermions play 
an active role. In it the phase-factor averages {e^^^-°) and {e^^'^) and the amplitudes go and g vanish. However, 
the fluctuations of the phases 9j,i remain small. As a result, the amplitude gi — |(e*^J'i)| remains finite and the 
short-range spin order prevails in the pseudogap state. The same applies to the pseudogap state reached at zero 
temperature, magnetic field in the range H € {Hq, H*), and a H-dependent hole concentration range above x — Xc2 
and below x — x^,. Here Hc2 and H* are the fields above which the phases 9j^i have large fluctuations, so that the 
amplitude gi = |(e'^^'^)| vanishes and there is no short-range spin order. 
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While the phase transition occurring a,t H = Hq results mostly from c fermion diamagnetic pair breaking, the 
disappearance of short-range spin order occurring at the Hc2 and H* lines results both from c fermion diamagnetic 
and Zeeman spin-singlet spinon si fermion pairs breaking. Consistently, at zero temperature the fields H = Hc2 for 
X G {xo,Xc2) and H — H* for x G {xc2,x*) are those at which Zeeman spin-singlet spinon si fermion pair breaking 
fully suppresses the short-range spin order and thus the diamagnetic c fermion pair breaking fully suppresses the 
corresponding incoherent c fermion pairing correlations. 

For hole concentrations in the range x G (xo,Xci) we use the method of Ref.** to derive the x dependence of the 
magnitude iJc2 of the uniform magnetic field H aligned perpendicular to the square-lattice plane. This involves 
expanding the quadratic terms of the continuum Lagrangian representing the charge sector associated with the c 
fermion pairs. The quadratic terms reflect the centre of mass motion of such pairs. They are expanded in charge 
— 2e Landau levels resulting from the diamagnetic orbital coupling of the magnetic field to the charge c fermions. In 
turn, we did not derive an expression valid for hole concentrations in the range x G {xci-,Xc2)- For that x range we 
could only access an inequality obeyed by Hc2 and the hole concentration Xc2^ respectively. In turn, the x dependence 
of the upper field H* is a simpler problem. Alike Hq ~ fcs Tc/^ib is proportional to the superconducting critical 
temperature Tc, for x G {xc2, x*) it is proportional to the pseudogap temperature T* of Eq. ([22|) . For the interaction 
range f//4i G (wojMit) our results for the fields Hc2 and H* are, 

Hc2 ~ J^ ^ —, :-- , X G {X0,Xcl); TT-, 7 -F < Hc2 < J^ , X G [Xcl,Xc2) , 

(X^JO) SlJ'B 9l[Xcl) 

k,Q 1 ^ Zaq 

H « «5i — , xG{xc2,x*). (51) 

Mb ^^B 

Here the amplitude gi = gi{x) is given in Eq. P^ . We call iJ*2 the field magnitude reached at the hole x = Xc2 at 
which Hc2 = H* . Hence HI2 = Hc2{xc2) = H*{xc2)- The inequalities obeyed by H*2 and Xc2 are then given by. 



[3"fo + l)nB V 370 + 1 



< H:2 < ,^ ', " ; X™" < x,2 < X, ; X™" = ( ^^^^ ] X, , (52) 



respectively. Here H*2 — 7oAo/[(37o -I- 1)/^b] and H*2 = for a; = x™" and x = x"^^ — x*, respectively, and 70 
is the parameter given in Eq. (H51) . It is expected that H*2 is closer to « 7oAo/[(37o + l)(Us] than to zero, so that 
Xc2 ~ a^™" and Hc2 ~ ^ both for x G {xq,Xci) and x G (xci,Xc2)- 

That the expression H* « ^i[Ao//zb] involves the zero-temperature and zero-field value of the amplitude gi = 
|(e*^ji)| is consistent with it marking the full suppression of the short-range spin correlations. Indeed, for H > H* 
and hole concentrations in the range x G {xc2 , x^, ) the system is driven into a disordered state without short-range 
spin order. In turn, for hole concentrations x G {xo,Xci) the field Hc2 grows linearly with x. If such a behavior is a 
good approximation for x G {xci,Xc2) as well, for x* = 0.27 one finds Xc2 ~ a;™" « 0.20 for the five representative 
systems. For x G {xi,Xc2) the actual H(.2 x dependence may slightly deviate to below the linear-x behavior. If so, 
the hole concentration Xc2 ~ 0.20 may increase to « 0.21 — 0.22. 

The pseudogap state occurring for fields below Hc2 and H* and above Hq for a iJ-dependent hole-concentration 
range centered at x = Xc2 is a quantum vortex liquid. It shows strong vortex fluctuations and enhanced diamagnetism. 
Along the line H — Hq, this quantum vortex liquid freezes into a solid. Then any small deviation —5H < of the 
field H = [Hq — SH] near Hq leads to superconductivity. For magnetic fields above the field H* corresponding to the 
range x G {xc2,Xt.) the system is driven into a disordered spin state. It is similar to that reached at zero field and 
zero temperature for hole concentrations larger than x*. Upon lowering H from above H* , the field H* marks the 
onset of the short-range spin order. Therefore, it marks a crossover rather than a sharp transition. The same applies 
to the pseudogap temperature T* of Eq. (P^ for vanishing magnetic field H — and finite hole concentrations in 
the range x G {xo,x^,). Within a mean- field treatment, the upper field Hc2 refers to a crossover as well. For finite 
fields H and vanishing temperature the true quantum transition to the superconducting phase takes place at the line 
Hq = Hq(x). Along it the quantum vortex liquid freezes into a solid. Whether beyond mean-field theory the line 
associated with Hc2 {x) marks a sharp quantum phase transition and the onset of a new order for fields H > Hc2 and 
finite hole concentrations in the range x G {xq^Xc2) remains an open question. 

IV. THE QUANTUM-LIQUID ENERGY FUNCTIONAL, VIRTUAL-ELECTRON STRONG 
EFFECTIVE COUPLING PAIRING MECHANISM, AND SUPERFLUID DENSITY 

Here we study the general energy functional associated with the YEP quantum-liquid microscopic Hamiltonian 
considered in Section III-C. It refers to that Hamiltonian in normal order relative to the to = initial ground state. 
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The 2D energy functional considered here accounts for the weak 3D uniaxial anisotropy effects through the dependence 
of the critical hole concentration Xc and pairing phases on the parameter e^ — m*^/M <^ 1. 

Such a functional is obtained from that constructed in Ref j^ on introducing in it the effects of the coherent virtual- 
electron pairing phase fluctuations, accounting for the small suppression effects. For x G (xcX*) this refers to the 
occurrence in the VEP quantum liquid of a long-range superconducting order. That program involves the evaluation 
of the general superconducting virtual-electron pairing energy, which emerges in the spectrum of such an energy 
functional. The derivation of such a pairing energy requires the study of the virtual-electron strong effective coupling 
pairing mechanism. This involves as well investigations on the energy range of the c - si fermion interactions behind 
it and signatures of virtual-electron pairs in the one-electron removal spectral- weight distribution. In addition, in this 
section we study the hole-concentration dependence of the zero-temperature superfluid density. 

A. The quantum-liquid energy functional 



The phase and amplitude of the complex gap function (j24l) appearing in the VEP quantum-liquid Hamiltonian terms 
([23|l show up in the superconducting order parameter 2il of Eq. (|40|1. The absolute value of this order parameter is 
the maximum magnitude of the superconducting virtual-electron pairing energy 2|f2si('z)| evaluated below in Section 
IV-B for si band momenta q « <1bsI ^^ *-'^ near the si boundary line and in Section IV-F for momenta belonging 
to the general coherent si — sc lines introduced in that section. (Here the index sc stands for strong coupling.) For 
temperatures below T^ and a T dependent x range centered at x = Xop, which is largest and given by a: £ (XcX*) 
at T = 0, the superconducting pairing energy per electron |nsi(g)| derived in these sections emerges in the si 
fermion energy dispersion £51(0). This dispersion appears in the general energy functional introduced in the following. 
That functional involves the ground-state normal-ordered c and si fermion hole momentum distribution function 
deviations^, 

5N^^{q^) = [N'^{q^) - N^,^\q^)] ; 5N':,%) = [^i\('S) " ^'l°(9.)] • (53) 

Here N^'^{q^) and N^-l (qj) are the corresponding initial-ground-state values. According to the results of Ref ji., the 
c and si fermion discrete hole momenta (f. and momenta qj, respectively, are good quantum numbers for the e^ = 0, 
ad = 0, and Xq = square- lattice quantum liquid of Ref^i. Under the very weak 3D uniaxial anisotropy effects and 
suppression effects they are for x £ (xc, x*) close to good quantum numbers of the VEP quantum liquid. This is in 
contrast to the range x € (0, Xc), in which the strong hole tapping effects reported in Appendix B change qualitatively 
the square-lattice quantum liquid physics of Ref.^. 

It follows that alike for that quantum liquid, for the the VEP quantum liquid considered here the hole momentum 
distribution functions N^{qj^) and iV^i(q^) approximately read 1 and for unfilled and filled, respectively, discrete- 
hole- momentum values qj^ or discrete- momentum values qj . Furthermore, the first-order energy terms of the energy 
functional given in the following correspond to the dominant contributions. The c — c fermion interactions vanish or 
are extremely small and can be ignored, whereas the si — si fermion interactions do not lead to inelastic scatterin g^i^*^ . 
In turn, the effects of the c— si fermion interactions of interest for the studies of this paper are implicitly accounted for 
in the pairing energy involved in the si band energy dispersion esiio) appearing in the first-order terms of the energy 
functional. This justifies why its second-order terms in the deviations (j53l) are not studied here. An approximate 
expression of the matrix element limp'_j.o Wc.si{q''',QtP) of the c - si fermion effective interaction between the initial 
and final states needed for the derivation of the one-electron inverse lifetime and related to the second-order forward- 
scattering fc,si{<i^:<t} function is for transfer momentum p — >■ evaluated in Ref.^. 

For the interaction range 11/4,1 G {uq^Ut^), hole concentrations x S (xc,x*), and vanishing spin density rn = the 
VEP quantum liquid general energy functional is to first order in the deviations (1531) given by, 

6E = -Y, e.{q-'^)SNi\q-'^) - J^ ^^M^N'^iq) . (54) 

g'» q 

The energy dispersion edciP") appearing here is given in Eqs. (|A10I) and (JA11|) of Appendix A. The dispersion esi{Qj) 
reads, 



e.i{q) = -^\eO,M\' + \^v-el{q)\^; \A,-el{q)\ ^ \A,im + l^sMl . (55) 

Here the energy dispersion e^i{q) and pairing energy per spinon |A5i(q)| are given in Eq. (|A12p of Appendix A, 
|^t>-ei(9)| is the virtual-electron pairing energy per electron, and |r25i((7)| is its coherent part introduced below in 
Sections IV-B and IV-F. The short-range spin order and its order parameter 2|A| refer to the range x € (xo,a;*). 
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However, due to the effects of the Anderson insulator behavior coexisting with that order for x G {xq,Xc), the si 
fcrmion energy dispersion esi{q) of Eq. (j55p and the corresponding momentum-dependent energies |Asi((7)| and 
1 0^1 (g) I are weh defined only for the range x € {xc,x^). 

As discussed in Section III-D, removal of two electrons in a spin-singlet configuration leads to an excited state with 
no holes in the si band. Hence it follows from the form of the energy functional (|54l) that no virtual-electron pair 
is broken under that excitation. Instead, such a pair becomes the removed electron pair. In turn, one-electron and 
spin excitations lead to the emergence of one and two holes in the si band, respectively. Thus their energy spectrum 
involves a finite pairing energy associated with virtual-electron pair breaking. 

We recall that the virtual-electron pairs carry the same momentum q as the corresponding si fermion. Indeed, 
their two c fermions have hole momentum g'' and —q'^, respectively. Therefore, the c fermion pair net momentum 
contribution vanishes. This justifies why the energy scale |A^_e/((7)| given in Eq. (j55p appears in the g dependent 

energy dispersion esi{q) — —\/\^^i{(l)\'^ + \^v-ei{q)\^ ■ Indeed, that it and thus |r2si('?)| appear in such a dispersion 
follows from the virtual-electron pair associated with these pairing energies per electron carrying momentum q. As 
discussed in Section III-D, such virtual-electron pairs of momentum q exist as individual objects only in intermediate 
virtual states of one-electron and spin excitations. Nevertheless, the virtual-electron pair energy is related to the 
spectrum of one-electron excitations and has physical significance, as confirmed below in Section IV-D. It reads, 

E,^,l{q) = 2\e,{±q'')\ + |e,i(g)| . (56) 

Note that such an energy is not in general a pairing energy and thus is different from the virtual-electron pairing 
energy per electron \Ay-ei{q)\ of Eq. ((55|) . It is found below that for momenta g « q^^^ at or near the si boundary 
line coherent pairing is associated with c fermions of hole momenta q^ and —q^ near the c Fermi line. In this case 
Ey-ei{(iBsi) ~ l^t)-ei(9Bsi)li ^^ that the virtual-electron pair energy becomes the virtual-electron pairing energy per 
electron. 

At zero temperature the energy |J7si((7)| is finite for si band momenta g corresponding to si fermions that mediate 
coherent c fermion pairing. These momenta belong to the range q G Q^^ defined in Section IV-F. It refers to a set 
of coherent si — sc lines centered at zero si band momentum. Moreover, si fermions whose momentum belongs to 
a given coherent si — sc line only interact with c fermions whose hole momenta belong to a related nearly circular 
coherent c— sc line centered at the c band momentum — tt = — [tt, tt]. There is a one-to-one correspondence between 
such si — sc lines and c— sc lines. 

The non-coherent |Asi(g)| and coherent |r2si(g)| parts of the virtual-electron pairing energy \/S.y-ei{q)\ of Eq. ([55)) 
are associated with si fermion spin-singlet spinon pairing and coherent zero-momentum c fermion pairing, respectively. 
That they appear in the spin si fermion spectrum and their maximum magnitudes | A| and |r2| = 7^ go I A| are closely 
related and as confirmed in Ref.^^ correspond to the superconducting- and normal-state maximum one-electron gap 
|A| and superconducting-state low-temperature magnetic resonance energy 2|ri|, respectively, is consistent with the 
interplay between magnetic fiuctuations and unconventional superconductivity discussed in Ref^. 

The corresponding general excitation momentum functional is linear in the deviations (j53p and reads, 

5P = Sq^ - Y,[q'^ - Tf] SN^^iq'^) - ^^ qSN'^iq) . (57) 

Here 6q^ is the subspace-dependent small momentum deviation considered in Refii. 

Except for the superconducting pairing energy per electron |J7si((7)| appearing in the expression of the dispersion 
esi(q} of Eq. (|55p. the first-order terms of the energy functional (j54l) are those derived in Refji for approximately the 
range U/At > mq ~ 1.302. Indeed and as found in that reference, for [//4i < uq the energy scale 2Ao whose limiting 
behaviors are given in Eq. (jA4l) of Appendix A becomes very small. It follows that the amplitude fluctuations of the 
order parameter of the short-range spin correlations cannot be ignored. As discussed in Ref.^, for small U/t values 
the si fermion energy dispersion esi{qj) is then not expected to have the form given in Eq. ([55]) . 

B. d-wave virtual-electron pairing mechanism near the Fermi line 

Here we evaluate the superconducting pairing energy per electron |J7si(9)| in the si band energy spectrum of 
Eq. (|55|) for momenta q ~ QbsI ^* '^^ near the si boundary line. As given in Eq. (IA13P of Appendix A, the 
order parameter 2| A| associated with the short-range spin correlations is the maximum magnitude of the si fermion 
spinon-pairing energy 2|Asi(g)|. Similarly, the absolute value 2 1 $711^=0 of the order parameter associated with phase- 
coherent-pair superconducting order is the maximum magnitude of the superconducting virtual-electron pairing energy 
. The maximum magnitude of the si fermion spinon-pairing energy 2|A<,i(g)| corresponds according to Eq. (IA13[) of 
Appendix A to si band momentum values q — qsti^ ■ Those belong to the si boundary line and their corresponding 
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auxiliary momentum values qo — QbsI given in Eq. (IA7p of that Appendix point in the anti-nodal directions. (The 
equality gj^^ = QbsI holds for ground states at vanishing spin density m — and their two-electron excited states. 
In turn, one has that q^fi^ ^ g^iT ^^^ their one-electron excited states^.) 

The VEP quantum-liquid d-wave long-range superconducting order occurring at zero temperature for the x range 
X e {xc,x^,) is a by-product of the short-range spin correlations. For small < (a; — Xc) <C 1 this shows up in the 
relation of the zero temperature energy scale |f2si(<?i3si)| to the si fermion pairing energy per spinon associated with 
such short-range correlations. Indeed for < (x — Xc) <^ 1 the energy scale |51si(q)| is the deviation of the si fermion 
pairing energy per spinon |Asi((f|^;^ - ggj| from its value |Asi(gJ^J| at gj^ = 0. Here qg^ = [go/v^Je^A. (The 
physical meaning of the vector qgg is discussed below.) Expanding |Asi((fJsi ~ Qgo)! gives, 

l^sliqLl - «§o)l « \^sl{qLl)\ + Vj^iQLl) ■ % = \^sliqisl)\ + \^sliqLl)\ . T^O. (58) 

Hence for < (a; — Xc) <C 1 the zero-temperature superconducting virtual-electron pairing energy per electron reads, 

1^^.1(9^.1)1 = Vj^iqLi) ■ qso = 9oV2V^,, ^\n\\T=o\sm2cf>\, T = . (59) 

Here the si fermion velocity V^lq^^i) is given in Eq. (IA15|) of Appendix A, V^^i denotes its absolute value of Eq. 
(|n]), and |r2||T=o is the superconducting energy scale of Eq. (P^)) . 

The expansion of Eq. ((55| refers to small < (a; — Xc) <C 1 and thus to small absolute values go [l/\/2] ~ 
[{x — Xc)/{x^ — Xc)] [l/v2] of the vector (fg,, = [go/v2] c^a . As mentioned above, that for such a range of x values 
near and above the critical hole concentration Xc the coherent virtual-electron pairing energy per electron |^si(9jsi)| 
associated with the phase-coherent-pair superconducting order is the energy deviation Vj^{qBsi) ■ qgo of the expansion 
([58|) is consistent with such an order being a by-product of the short-range spin correlations. Indeed, the latter are 
associated with the si fermion pairing energy per spinon |Aii((f^j,]^)|, which corresponds to the zeroth-order term of 
that expansion. 

For larger hole concentration values in the range x £ (a;c,x*) the expansion of Eq. (j58p is not valid anymore. 
However, within our scheme the expression for the zero-temperature coherent virtual-electron pairing energy per 
electron |r2si((fjj,j^)| w 7d5o V^Vg^i of Eq. ([59l) multiplied by the suppression coefficient 7^ of Eq. (|35|) is. Indeed, 
that equation refers to the limit < (x — Xc) <C 1 for which 7^ = 1. For virtual-electron pairs of momentum g « qBsi 
the full virtual-electron pairing energy per electron has the following form, 

\K-ei{qisi)\ - \^si{qLi)\ + 1^.1(9^.1)1 - |A|| COS 201 + l^ll sin20| . (60) 

For small < (a; — Xc) ^ 1 this energy reads \Asi{q^gi — qgo)\, as given in Eq. ([55]) . The form of the pairing energy 
(|60p is consistent with virtual-electron pairing having a d-wave character and involving both spin-singlet spinon si 
fermion pairing and c fermion pairing. For the normal-state ground states at finite magnetic field H considered in 
Section III-E the virtual-electron pairs exist yet loose phase coherence. Then the virtual-electron pairing energy per 
electron equals for virtual-electron pairs of momentum (f w qBsi the spin-singlet spinon si fermion pairing energy for 
the same momentum, \Ay_ei{qBsi)\ = I'^siiqisiM- 

It follows from the above analysis that for hole concentrations in the range x € {xc,x^) the zero-temperature 
absolute value 2|ri||r=o defined by Eqs. (|5^ and ([33]) of the order parameter introduced in Eq. PO)) associated 
with the phase-coherent-pair superconducting order reads 2|0||t=o = '^l^siiqBsi )l- It corresponds to the maximum 
magnitude of the coherent virtual-electron pairing energy. Such a maximum magnitude refers to virtual-electron pairs 
of momentum at the si boundary line and whose auxiliary momentum points in the nodal directions. We recall that 
the maximum magnitude of the zero-temperature spin-singlet spinon si fermion pairing energy max2|Asi((fj5]^)| = 
2|Asi(q^g]^ )| = 2|A| corresponds to momenta belonging to the si boundary line and whose auxiliary momenta (Jb^ 
point instead in the anti-nodal directions. 

For hole concentrations in the range x e (a;c,a;*), momenta q ~ q^^i at or near the si boundary line, and 
temperatures T smaller than Tc the coherent part |J7si((7)| of the pairing energy per electron reads, 

\^siiq)\ = \^sii-q)\~7dV,^iq)-qgo=7dgQV2V,^{q); qg^=goe^A; q-qLi- (61) 



It involves the general amplitude go rather than its T = magnitude go = {x — Xc)/{xi, — Xc) of Eq. (1421) . 

The c - si fermion residual interactions are behind the effective pairing coupling of the two c fermions of hole 
momentum (f'* and — g'' such that (f'' « qp^ is at or near the c Fermi line. Such c fermion pairs are associated 
with the pairing energy per electron |f2si((7)| « 7^ V^{q) ■ qg^ of phase-coherent virtual-electron pair configurations of 
momentum q k, q^^i at or near the si boundary line. The short-range spin correlations provide the energy |f2si((7)| 
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through the interactions of the si fermion of momentum (f « QbsI^ energy esi{qBsi) ~ — |Asi((fBsi)|, and velocity 
V^{q) ~ — Vq-|Asi(g)| with the two c fcrmions. 

The vector qg^ plays the role of coupling constant of the c fermion- pair degrees of freedom to the si fermion of 
velocity Vj^{q). It is through it that the pairing energy |51si(q)| w TdKaC*?) ' 9sn needed for the c fermion effective 
coupling and corresponding coherent virtual-electron pairing is supplied to the two c fermions. The absolute value 
go/V^ of the vector qg^ involves the amplitude go = g/gi of Eq. ([M|) . It measures the relative strength of the coherent 
virtual-electron pairing and si fermion spinon pairing. Indeed, the magnitudes of the amplitudes g and gi also given 
in Eq. p4p provide a measure of the strengths of the coherent virtual-electron pairing and si fermion spinon pairing, 
respectively. 

As confirmed in Section IV-D, breaking under one-electron excitations of the virtual-electron pairs considered 
here refers to removal or addition of electrons from and to the Fermi line. Therefore, virtual-electrons of si band 
momentum at or near the si boundary line are within the electron representation at or near the Fermi line, respectively. 
Consistently, the superconducting fluctuations lead to an additional contribution |i^si('?Bsi)l ~ |i^|| sin2(/)| to the 
anisotropic term 5Ep — \Asi{qg^^)\ of the Fermi energy introduced in Refji, which appears in the spectra given in 
Eq. (I20[) . Hence for the present VEP quantum liquid the Fermi energy reads, 

Ef = fi + SEpicj)); n^il° + w!^, 
5Ef{4>) = -esiiqisi) = \A,^ei{qLi)\ = \^si{qLi)\ + \^si{qLi)\ ~ |A|| cos20| + \n\\ sin2</>| , (62) 

where the superconducting energy scale |51| is provided in Eq. (HT)) and VF^ is the ground-state c fermion- hole 
energy bandwidth given in Eq. (lAlip of Appendix A. For intermediate U/At values approximately in the range 
C//4i e {uq,ui) the W^ expression provided in that equation is a very good approximation for hole concentrations 
X e {xc,Xci) and a reasonably good approximation for the x range x € {xci,x»). Rather than the square-lattice 
quantum-liquid expression /i « /io + W^!', the expression ij, ra fiP + W^ given here for the chemical potential refers to 
the VEP quantum liquid. In it /i" = ^" + Sfi where /Lt° = lima;_>.o fi is one half the Mott-Hubbard gap of Eq. (jASp of 
Appendix A and the shift 6fi generated by the hole trapping effects is given in Eq. (IBip of Appendix B. 

As discussed in Ref."'^, in the Fermi energy expression Ep — fi + 5Ep{(j)) the chemical potential fi arises from the 
isotropic c fermion energy dispersion and 5Ep = |esi('fJsi)| — \^v-ei{qBsi)\ stems from the anisotropic si fermion 
energy dispersion. The extra term |51si((fJ^j^)| sa |J7|| sin2(/)| appearing for x G {xc,x^,) in the expression of the 
anisotropic Fermi energy term SEp given in Eq. (j62p does not change the physics discussed in Section III-A. In turn, 
for X € {xo,Xc) and T = the equality SEp = |A„_e;(gJ^-^)| — \Asi{q^gi)\ holds. The anisotropic Fermi energy 
term SEp plays an important role in the VEP quantum-liquid physics. Its maximum magnitude Tiiax{6Ep} = |A| 
determines and equals that of the anti-nodal one-electron gap SEp'^ — ina,x{SEp} — |A|. From the use of the 
critical-temperature expressions given in Eq. (l43l) and (|45|) and 5Ep[(f)) expression provided in Eq. (|62|) . SEp^ can 
be expressed as a function of Tc/T"''"'^ as, 
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SE^^ = max{dEp} = \A\{x) = |A|(a;op)(l - sgn{x - x^p}! ^1 - T,(x)/T,"--|) ; \A\{xop) =lc^- (63) 

C. The general virtual-electron strong effective coupling pairing mechanism 

According to our approach, c fermion strong effective coupling is that whose breaking upon one-electron excitations 
leads to sharp features in the one-electron spectral function. It is a necessary but not sufficient condition for the 
occurrence of coherent pairing in the initial ground state. The absolute value q^ = |g''| of the hole momenta q'^ and 
—q^ of two c fermions with strong effective coupling belongs to a well-defined range q^ € {Qpc^Qec)- Here g^^ is the 
maximum hole momentum given below for which there is c fermion strong effective coupling. We denote by Qg^ and 
Ql]. the corresponding c and si band momentum domains for which the c and si fermions, respectively, participate 
in virtual-electron pairs with strong effective coupling. Q^^ and Q^]. correspond to a set of c — sc and si — sc lines, 
respectively. As further discussed below, within the present scheme: 

i) Only strong effective coupling, i.e. that whose breaking upon one-electron excitations leads to sharp features in 
the one-electron spectral function, can lead to coherent virtual-electron pairing. 

ii) For each c band hole momentum absolute value g'* in a well-defined range q^ G (qpciqec) there is a c band 
approximately circular c— sc line of radius g'' centered at — tt. Such a range refers to the c fcrmions of hole momenta 
q^ and —q^ whose effective coupling is strong. There is no one-to-one correspondence between the c fermion pairs and 
the si fermions of strongly coupled virtual-electron pairs. As justified in this section, the strong effective coupling of 
c fermions of hole momenta q'^ and —q^ at or near a, c— sc line of radius q^ results from interactions with si fermions 
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of momentum q at or near a corresponding si — sc line centered at the si band zero momentum. The four nodal 
momenta belonging to the latter line have the same absolute value g^^, — q^^^{q^), which is uniquely determined by 
the radius magnitude q^ of the corresponding c — sc Hne. Here g^^, S (q^c^lBsi) where g^^ = g^ and g^^ = q^^^ 
correspond to q^ = g^^ and q^ = q%^, respectively. For g^^ < q^^^ and thus q^ > q'p^ the si — sc hne is constituted by 
four disconnected line arcs centered at zero momentum and crossing perpendicularly the four nodal lines, respectively, 
whose auxiliary momenta refer to Fermi angles 4> = 7r/4, 37r/4, 57r/4, 77r/4. The angular width of these four line arcs 
is a decreasing function of g^, which vanishes at g'' = g^^. As found below, in that limit these four line arcs reduce 
to four nodal momenta whose absolute value is g^. The function g^^ = Qardl'^) ^^'^ the magnitudes of q^^ and 
g^ are given below. That the number of si fermions that mediate the strong effective coupling of c fermion pairs 
of hole momenta g'' and —q^ decreases upon increasing g'\ vanishing for g'' > g^^, is consistent with for the range 
q^ € (^Fc lee) the strength of such effective coupling being a decreasing function of the c fermion absolute energy, 

kc(g-"')| = ([g^]^~[g^J^)/2r<. 

In turn, weak effective pairing coupling between c fermions is that whose breaking upon one-electron excitations 
leads to broad flat features. Weak effective pairing coupling never contributes to coherent virtual-electron pairing. 
Strong effective pairing coupling may or may not lead to such a coherent pairing. 

In the following we justify and supplement the information reported in point (ii). The ground-state ec-pairing line 
is the c— sc line of largest radius q^ = g^^. It separates in the c momentum band the c fermions with strong effective 
pairing coupling from those with only weak effective pairing coupling, respectively. (The designation ec-pairing line 
follows from it separating c fermions with two different types of effective coupling [ec].) The hole momentum domain 
Qgj, refers to the set of c — sc lines. It is limited by the c Fermi line and ec-pairing line, respectively. As confirmed 
below in Section IV-E, the energy scale, 

2m* 7e ' "■ ' 

is the maximum magnitude of both the energy bandwidth of the superconducting-ground-state sea of c fermions con- 
tributing to coherent pairing and the energy bandwidth corresponding to the hole momentum domain Q^^ considered 
here. 

Alike the smaller c — sc lines enclosed by it, for the interaction range f//4i e (uq, ui) the ec-pairing line is nearly 
circular for hole concentrations x e {xc,Xci). It is a fairly good approximation to consider that it remains circular 
for X e (xcijX*). The energy range |ec(g''')| = |ec(~g'')| G (0, IVec) of c fermions having strong effective coupling 
in the ground state corresponds to the uniquely defined hole momentum range g'' S (qpciQec) reported in (ii). For 
these c fermions, the hole momenta q^^ belonging to the ec-pairing line are those of largest absolute value g^^. For 
the interaction range U/At S {uq,ui) and x € {xc,x*) it is approximately given by, 

(65) 

An interesting property is that for x € {xcX^.) the c band hole momentum area corresponding to c fermions with 
strong effective coupling in the ground state is independent of x and given by, 

S!e = ni[q':f - [q'^f) = ^^° . (66) 

[X^ Xc) t 

It refers to the hole momentum area of the domain Qec- Its expression involves only the superconducting dome hole 
concentration width (x* — Xc) and the basic energy scales Aq and t. 

In Appendix D it is found that for U/At g (uo,ui) and hole concentrations in the range x € (xc,x*) the energy 
bandwidth |ec(g''')| of c fermions of momenta g^ and — g^ having strong effective pairing coupling due to residual 
interactions with si fermions of momentum q G Qll obeys the inequality, 

kc(g")|<W^ecfl-^f^) -Wec{l-\cos2cP\); | A,i (g") | « | A| | COS 2<^| for g e Q^J . (67) 




For x S {xc,x^,) the inequality Wee < W^ holds. (Here VFI" is the energy bandwidth of the ground-state c momentum 
band filled by c fermions.) Hence only si fermions whose auxiliary momenta point in the nodal directions have 
interactions with strongly coupled c fermions whose energy |ec(g'')| belongs to the whole range |ec(9'')| € (0, Wee)- 
The inequality (j67p is obeyed by the following hole momentum range, 



g'G(4c,'z'); 9' = «,?Fc; «, = W i + ("L - i) (^ - ^^) ^ ^/i + («L - i) (i - I cos 20|) , (68) 
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Equivalent to that inequality are the following inequalities obeyed by the magnitude |Asi(g)| of the si fermion pairing 
energy per spinon and Fermi angle defined in Eq. (jA6l) . 



|A.r(.1|<|A|(l-Mi^); | cos2^| < (l - ^0) , (69) 

respectively. The second inequality given here involving the Fermi angle is in turn equivalent to restricting the 
range 4> e (0, 7r/2) to an arc of angular width 2(j>arc centered at (/) = 7r/4, 

e (7r/4 - <j,arc, 7r/4 + <j,arc) ', <Parc = ^ arcsiu f ^££L_^_i_ \ ^ _ ^^^^j^^ [f^] ^ (0, ^/4) • (70) 

Here e^'^{q'^) is the following c fermion energy dispersion, 

e-(g'') ^ Wee + eeiq'') ; < ef{q'') < Wee ■ (71) 

It is the energy of a c fermion of hole momentum g '' measured from the energy level of the ec-pairing line. The energy 
(|7ip is largest at the c Fermi line and vanishes at the ec-pairing line, 

Citc)=0; C{<iFc)=Wee. (72) 

The shape of the ec-pairing line is for x e {xe,x■^,) fully determined by the form of the c band energy dispersion 
^Ti'i'^) as follows, 

ciee e ec - pairing hue ^^=^ ^fi^ee) = . (73) 

For simplicity we limit our analysis to the quadrant referring to the Fermi angle range 4> G (0, 7r/2) and corresponding 
si— sc line arc whose angular range is that given in Eq. (|70p . G (7r/4— ^arc, t^ /'^-'r't'arc)- To each c band approximately 
circular c — sc line of radius q^ centered at — tt corresponds such a si band si — sc line arc whose momenta q^^^^ are 
approximately given by, 

q^ 

qie{4>) = Ai^{4>) qarc{4>) ; qarc{4>)~^§^qBsl{4>)\ garc(0) = garc(7r/2 - (/)) ; (/. G (7r/4 - <?!)a^„ 7r/4 -^ (/)arc) ■ (74) 

Here qarc{4>) ~ [Qarc/^Bsil lBsi{4>) is the absolute value of the auxiliary momentum qarc{4>\ ^fi i^ the matrix of Eq. 
(O, and the si — sc hue arc nodal momentum absolute value g^^, — qareil^) i^ given below. As stated in (ii), for each 
c band c— sc line of radius q*^ in the range q^ G {Qfc 'iec) there is for (p £ (0, 7r/2) exactly one of such si band si — sc 
line arcs. The line arc limiting angular widths 2(parc — and 2<j>arc = 7r/2 refer to q'^ — g^^, and q^ — qp^, respectively. 
The angle range </> € (7r/4 — (pare, 7i'/4 + 4>arc) of each si — sc line arc runs symmetrical around (j) — 7r/4. The physical 
importance of the si — sc line arcs follows from the corresponding sharp spectral features line arcs studied below, 
which emerge in the one-electron removal weight distribution and are observed in experiments on the cuprates^"^. 

We note that for € (0, 7r/2) the si band auxiliary momenta belong to the quadrant for which their two components 
are negative^. For instance, for momenta near the si boundary line and x G {xci,Xc2) this follows from the relation 
of the si band momentum angle 0si = ^-I-tt of Eq. ((TT]) to the Fermi angle </>. Indeed the range (psi S (tt, 37r/2) refers 
to (/) e (0, 7r/2). Consistently with Eq. ([71)) . the si boundary line is for (p G (0, 7r/2) the si — sc line arc whose nodal 
momentum absolute value is given by g^^ — q^^i- It corresponds to the approximately circular c Fermi line whose 
radius reads q^ = q^^: One confirms below that g^^ = qg^i for g'* = q^^. Hence the auxiliary momentum angle 
(psi — (j) + T^ can be generalized to all si — sc line arcs. Each of such line arcs then refers to an auxiliary momentum 
angle range 0si e (57r/4 - (pare, 57r/4 -|- (pare) corresponding to G (7r/4 - (pare, 7r/4 -I- (pare)- 

As mentioned above, there is no one-to-one correspondence between the c fermion pairs and corresponding si 
fermions participating in virtual-electron pairs with strong effective coupling. Concerning general c— si fermion inter- 
actions, there are in average two c fermions for each si fermion. However, only part of these c fermions and si fermions 
participate in strongly coupled virtual-electron pairs. There are four classes of strongly coupled virtual-electron pairs 
associated with the four disconnected si — sc line arcs, respectively. For any strongly coupled virtual-electron pair 
configuration whose auxiliary momentum Cartesian components are [qarc cos (psi , qare sin (psi] where (psi = (P + tt there 
are in the other three disconnected si — sc line arcs three strongly coupled virtual-electron pair configurations, respec- 
tively, with the same energy and momentum absolute value qarc- The si band auxiliary momenta of the four strongly 
coupled virtual electron pair configurations under consideration have Cartesian components ± [qare cos (psi, qarc sin (psi] 
and ±[— garcSini/)^!, garcCOS0si]. The set of four si — sc line arcs that such four momenta belong to are associated 
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with four sub-domains Qf^^'^^ and Qf^^''^^, respectively, of the si band momentum domain Qll of strongly coupled 
virtual-electron pair configurations. 

The point is that the strongly coupled virtual electron pair configurations whose auxiliary momentum has Carte- 
sian components ±[qarc cos 4>si,qarc sin 4>si] couple to different c fermion pairs than those of Cartesian components 
±[— Q'arc sin 0si,qarc cos (/isi]. Hence one can divide the c band domain Q^^ into two sub-domains that in average have 
half of the c fermion pairs each. We call them Q'^^ and Qjc^ respectively. The strong effective coupling of c fermions 
of hole momenta belonging to Q'^^ and Q~^ results from interactions with si fermions whose momenta belong to 
Qtc^'^^ and Qfc^'^^, respectively. Furthermore, symmetry arguments concerning the relationship of the two strongly 
coupled virtual-electron pairs with the same energy, momentum absolute value, and momentum direction but different 
senses suggest that the number of si fermions of each of the two si band sub-domains Qf^'^'"'^ (and Qj^' '*^) ^^ *^^^ 
half that of c fermions in the c band domain Q+f (and Qec)- Therefore, for each of the four sub-domains Qfc^''^^ and 
Oec'~'^^ the si band momentum area S^^ corresponding to si fermions of virtual-electron pairs with strong effective 
coupling in the ground state is for x G (xc, a;*) one fourth that of Eq. ([66)) . 



(X* - Xc)t 

It follows that the momentum area of the whole domain Qll equals that of Q^^, and reads 451^ — 4Ao/(a;* — Xc) t. It 
refers to the four sets of concentric si — sc line arcs centered at 0. The present results are valid for U/4t e (uq, ui). 

The simultaneous validity of the equalities (|66p for the whole c band domain Q'^^ and ([75]) for each of the four 
sub-domains of Q^l requires the hole concentration to obey the inequality x < x^^^ = 1 — [2Ao]/[7r^(a;* — x^i\. 
Hence x* must be smaller than or equal to x™^^ . This requires that Ao/i < [7r^/2](l — a;*)(a;* — x^. For x^, — 0.5 
and x* = 0.27 this gives approximately Ao/t < 0.79. The fulfillment of this inequality assures that for x S (xc,a;*) 
the total numbers of c fermions and si fermions are larger than those of c fermions and si fermions, respectively, 
participating in strongly coupled virtual-electron pairs. For t//4t = 1.525 and the t magnitude t « 195 meV found 
in Reffi^ to be appropriate to the four representative hole doped cuprates other than LSCO for which Xc = 0.5 and 
x* = 0.27 one finds Aq ~ 84meV, so that Ag/t ~ 0.28. In turn, the LSCO cation-randomness effects considered in 
that section lessen the Aq magnitude to Ao ~ 42meV, so that Ag/i ~ 0.14 for that random alloy. 

That in average the number of si fermions whose momenta q belong to the domain Q^;!. is larger than that of c 
fermion pairs of hole momenta q^ and —q^ belonging to Q%^ is consistent with the dependences on g'' of the si — sc 
and c ~ sc lines length. The c — sc line length sa 2Trq^ increases upon increasing q'^ . In contrast, that of the si — sc 
line, 

L,^{qZJ^4h^{q^^J■, hiiq^J ^ 2(1 - sm2cj^arc)\/[qi^,? - [q^J' = 2 ^ff ~ '^°;-^ J[qi^,V - [<J^ , (76) 



decreases. Here lsi{qarc) is the approximate length of each of the four si — sc line arcs, g^ is the minimum absolute 
value of the nodal si — sc-line-arc momenta q/^^ = q^^^{7T/4) and their auxiliary momenta q^ = garc(7r/4) of Eq. ([71| . 
and q^gi and g^^ are the absolute values of Eq. (lEll) of Appendix E of the nodal and anti-nodal si boundary-line 
momenta qi^i and q^f/^ , respectively. According to the definition of Ref.^, their corresponding auxiliary momenta 
?Bsi ''^^'^ ^Bsi' respectively, have for the quadrant for which e (0, 7r/2) so that (pgi £ (vr, 37r/2) and thus qoxi !i 



and qox2 < the Cartesian components given in Eq. (jASI) of Appendix A. The x dependence of qg^i and q^^i is 
discussed in Appendix E. 

Hence the number of si fermions available to supply the energy needed for c fermion strong effective pairing coupling 
decreases upon increasing the c— sc line radius q^. For q^ equal or close to q^^ there are plenty si fermions to supply 
the energy needed for the c fermion strong effective coupling. In the opposite limit of q'' tending to (jg^ the c — sc 
line length lirq^^ is maximum. However, the length Lsi{q^^^) given in Eq. (j76p tends to zero, as each si — sc line arc 
becomes a single discrete momentum value. Therefore, then only two sets of two si fermions whose auxiliary momenta 
have Cartesian components ±[g^/-\/2, g^/-\/2] and ±[— g^/-\/2, q^^/^/2] are available to supply the energy needed for 
the ground-state strong effective coupling of c fermions of hole momenta belonging to Q^^f and Q~'^, respectively. 
Consistently, there is no strong effective pairing coupling for g'' > g^^ and g^.^, < g^. 

The momentum area S^l of Eq. (|75p is approximately given by. 



Ibs 



5fc= / dqZjsliq':rc)^WBsl~C)^Jki7l?-kBsl?- (77) 

Here lsi{qarc) is the approximate length of each si — sc line arc given in Eq. ()76p . From the use of the two S^^ 
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expressions provided in Eqs. ([75]) and ([77]) . respectively, one straightforwardly arrives to, 

€ - iL, , ^° ■ (78) 

On combining this expression with the auxiliary si — sc-line momentum expression of Eq. ([74]) one confirms that, 

Qarc((/') = qBsii<P) for gf^^ = q^^^ , e (0,7r/2) , 

= C(0) for qarc = 'Ze^'c , </- = 7r/4 . (79) 

Here q^ is the auxiliary momentum of the si — sc- line momentum qf^'^ = A'^-^ q^ of smallest absolute value (|78p . 

The absolute value g^^ € ilecilBsi) '^^ the nodal momentum belonging to a given si — sc line arc is obtained by 
replacing in expression ([78| the energy scale Aq by the energy spectrum \Aec{q'^)\ defined in Eq. (|D1|) of Appendix 
D. As given in Eq. (jD2p of that Appendix, its magnitude varies from |Aec(g)| = for q'^ = qp^ to |Aec(g'')| = Aq for 
q^ = q^^. This is consistent with c fermions with strong effective coupling and hole momenta at or near the c Fermi 
line interacting with si fermions of momenta at or near the si boundary line. In turn, c fermions with strong effective 
coupling and hole momenta at or near the ec pairing line interact with the si fermion whose nodal momentum has 
absolute value g^ given in Eq. ([78]). In the limit q^^ — >■ g^ the si — sc line arc of minimum absolute nodal momentum 
q^ has vanishing length and for e (0,7r/2) and thus (/)si S (tt, 37r/2) reduces to a single discrete nodal momentum 

ffdN _ Ad ffN 
lee — -^sl lee ■ 

Replacement in expression ([78]) of the energy scale Aq by the energy spectrum |Aec(<f'')| then leads to. 



Qare'^lBsi ' '^ " - = g^^i - — ^ ^ ^^^ ^^ TT ; g-^ G (g^„ 9,^J . (80) 

These expressions are equivalent to, 

Qare « QbsI - (1 " sm 2<j>,,.eM,, ~ gfj . (81) 

The angle (pare of Eq. (1701) can be expressed in terms of the nodal momentum g^^ of Eq. ([50]) as follows, 

^are = I arcsin ( ^f^^^^) e (0, 7r/4) . (82) 

As given in Eq. ([8T]) . the absolute value g^^ provided in Eq. ([80]) of the nodal momentum belonging to the si — sc 
line arc can be expressed as a function of the angle 4>arc- 

Consistently with the limiting behaviors reported in Eq. (|79p. the si — sc line arc considered here becomes for 
q^rc — Qee ^nd thus for q^ = g^^ a single discrete nodal momentum ql!;J^ whose absolute value is given in Eq. ([78]) . 
Hence q^re — Qee^ ■ Ii^ the opposite limit reached at q^rc — QbsI and then g'' = q^^ the si — sc line arc becomes 
the si boundary line of angular range cj) G (0,7r/2), so that q^rc = 'IbsI- Between these two limits the angular range 
4> S (7r/4 — 4>arc, 7r/4 + (pare) of the sl — SC line arcs increases from a single angle (j) = 7''/4 for g^^ = g^ and (Jarc = Qec^ 
to a maximum angular width p e (0,7r/2) for g^^ = g^^^ and (Jarc = QbsI- This corresponds to auxihary momentum 
angular ranges increasing from a single angle (pgi — bn/A to a maximum angular width (p^i G (tt, 37r/2), respectively. 

The inequality ([67]) is equivalent to restricting the energy dispersions f-%'^{q'^) of Eq. ([7T]) and ec{q'^) to the following 
ranges, 

C{q^) e {Vgf{q), W^ec) ; e^iq'^) G (-[M^ec - Kc^^'(9)], 0) , (83) 

respectively. Here the energy [Wee — Vec (*?)] i^ positive or vanishing. For a c fermion of a c strongly coupled pair 
whose energy e^'^(g'') is measured from the ec-pairing line the energy scale, 

Kc-'^(9) = -^|Asi(9)|; Kc^^(g-)e(0,I^,,), (84) 

with limiting values, 

Vei^iO = V:ff{qEst) = V!if{qi!['') = Wee , 
plays the role of an effective potential energy. 
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Since the momentum g of a strongly coupled virtual-electron pair configuration equals that of the corresponding si 
fcrmion, q is at or near a si — sc line arc. It follows from Eq. (p7)) that the energy Ey-ei{q) given in Eq. ([56)1 of the 
corresponding virtual-electron pair obeys the inequality, 

E^,-eiiq) < i?i(0) + |esi(<ratt)l ; ^i(0) = 2I^ec(l " | cos20|) . (85) 

Here q^^^ is a si — sc- line arc momentum q^^.^ pointing in the directions defined by the angle (p — t^/^ i (pare- Hence 
breaking virtual-electron pairs of momentum q — q^^.^ belonging to a si — sc line arc under one-electron excitations 
leads to sharp spectral features provided that their energy obeys the inequalities of Eq. (j85|) . Moreover, only such 
virtual-electron pairs may have phase coherence. Note that the Fermi angle range (p S (7r/4 — (pare, 7''/4 -I- (pare) of a 
si — sc line arc has been constructed to inherently the inequalities of Eqs. (j67[) and (j85p and the energy ranges of Eq. 
(|83| being obeyed. 

For U/At g {uq,ui) and x € {xa,Xc2) where xa ~ x^./2 the maximum magnitude of the energy |esi(g^rJ)| is often 
below the resolution of experiments on hole doped cuprates^. In that case the inequality ([55]) approximately reads, 

E,-el{(iarc) < El^l^) = 2T4^ec(l - | COs2,/)|) . (86) 

D. Virtual-electron pairing breaking under electron removal: Signatures of the virtual-electron pairs in the 

one-electron spectral-weight distribution 

Taking into account the ranges of Eqs. (1571) . ([55t . and ([55)1 and using the general energy functional of Eq. ([Ml) , one 
finds that the maximum energy of the one-electron removal processes that lead to sharp spectral features is reached 
as follows: Upon breaking a virtual-electron pair, a c fermion pair whose c fermions have hole momenta tf'' ' and —([^ ' 
at or near the c — sc line corresponding to the energy \ec{±(j'^ ')\ — Wee (1 — | cos20|) is broken. Simultaneously, a 
spinon pair of a si fermion whose momentum (f = (f^^^ is at or near the si — sc line arc whose nodal momentum 
absolute value q^^ is given by Eq. (j80p with q'^ = q'' ' = |(f'"| is broken as well. The maximum energy condition 
imposes that such a si — sc line arc momentum corresponds to a minimum 7r/4 — (pare or maximum 7r/4 + (pare Fermi 
angle (p magnitude. Furthermore, one of the two involved c fermions recombines with one spinon within the removed 
electron. The spinon left behind then leads to the emergence of one si band holei whose momentum cj" — g^^^ is that 
of the broken si fermion. Finally, the second c fermion goes over to the c Fermi line. The resulting energy spectrum 
is straightforwardly obtained by use of the hole-momentum-distribution-function deviations, 

SNJ^iq'') = [S^H^_^H , + S^H^^H , - S^H^^^^4 ; SN^M = %,-..^^ , (87) 

in the energy functional of Eq. ([5^ . It is given by the energy scale Ei{(p) of Eqs. (1551) and ([55)) . Such an energy 
scale appears on the right-hand side of the inequality E^-eiiQare) ^ Ei{(p) also provided in the former equation. It 
equals the maximum magnitude of the virtual-electron pair energy Ey-eiiQare) of strongly coupled virtual-electron 
pairs of momentum q = q^^^. It follows that the one-electron spectrum Ei{(p) corresponds to a hue in the (0, w) 
plane associated with the boundary separating the one-electron sharp features from broad incoherent features. Such 
a boundary refers to the equality E^-ei{<larc) ~ Ei{(p) associated with virtual-electron pairs whose si — sc-line arc 
momentum has (p limiting magnitudes given by 7r/4 ± (pare- In turn, the above inequality E^-eii<lare) — Ei{(j)) of Eq. 
(|67p corresponds to a {(p,uj) plane domain bounded by it and associated with virtual-electron pairs whose si — sc-line 
arc momentum has (p magnitudes belonging to the whole corresponding range (p £ (7''/4 — (pare, 7''/4 + (pare)- Indeed, 
breaking virtual-electron pairs of momentum (f ss (j^^,^ upon one-electron removal excitations leads to sharp spectral 
features provided that their effective coupling is strong. This implies that their energy obeys such an inequality. It 
corresponds to a {(p, uj) plane domain constituted by a set of sharp one-electron removal spectral features lines. 

There is an one-to-one correspondence between each sharp one-electron removal spectral feature line arc (1-el-sharp- 
feature line arc) and a virtual-electron-pair si band si — sc line arc. Specifically, for x € {xa, Xc2) and U/At e (uq, ui) 
the energy scale |esi(q'jf^g)| is small and below the experimental resolution, so that the energy and hole momentum of 
the 1-el-sharp-feature line arcs read. 



E{k'') « ^— ^ £,^el{2(j>are) ^ 2Wee{l - siTi2^are) , 

IT m* 

k >^ = (i^^^ ~ tare , ^ (vr/4 - (Pare, 7r/4 + (Pare) ', k"" ^ k^ "^ - 1^ = - ^^(^Ze^ ~ ^Ce , (88) 

larc IBsI 



respectively. Importantly, the energy scale EiW) k, E{q^^^) exactly equals the energy Ey-^i of the strongly coupled 
virtual-electron pair broken under the one-electron removal excitation. Except for the small energy scale |esi((?j 
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ignored here, strongly coupled virtual-electron pairs associated with the same 1-el-sharp-feature line arc have the same 
energy. Indeed, Ey-e.i{<larc) depends only on the nodal momentum absolute value g^^ of the corresponding si — sc 
line arc. It may alternatively be expressed as a function of the angular width 2(j}arc of the line arc under consideration. 

Note that different 1-el-sharp-feature line arcs have different energy magnitudes E{q^^^) = Ey-ei{qarc)- Hence 
they refer to broken strongly coupled virtual-electron pairs of different energy. It follows that the energy E(q^^^) = 
Ey-e-iiQarc) of each 1-el-sharp-feature line arc provides a direct experimental signature of the corresponding broken 
strongly coupled virtual-electron pairs. Its energy exactly equals that of such pairs. The £^-ei{'24>arc) expression 
provided in Eq. (|88p refers to the energy 2|ec(g'')| expressed in terms of the corresponding si — sc-line arc angular 
width 2(t)arc- The coefficient (1 — sm2(f)arc) equals the factor (1 — | cos2(/)|) of the Ei expression given in Eq. ([62]) for 
2(f) = (^/2 ± 2c^arc) and 20arc e (0, 7r/2). 

The 1-el-sharp-feature line arc excitation hole momentum k^ — k+if also given in Eq. (j88p is centered at momentum 
— TT. In that expression an expression for the corresponding momentum k valid only for the nodal direction is also 
provided. Each of the four arcs of the 1-el-sharp-feature line has exactly the same angular width 2(j)arc as the 
corresponding si — sc line four arcs. Moreover, the Fermi angle (f) range provided in Eq. (|88p of each 1-el-sharp- 
feature line arc has been constructed to inherently the inequalities of Eqs. (p7|) . ([55)) . and ([55]) being obeyed. Its 
expression provided in that equation can be expressed in terms of the momentum q^,,^ of the hole that emerges in the 
si band and final hole momentum qp^ of the c fermion transferred from the c band hole momentum tf'' or —q'^ to the 
c Fermi line. The auxiliary momentum qarc — [^fi]~^ Qarc ^as been constructed to inherently pointing in the same 
direction as the hole momentum —k^. This implies that the absolute value k^ of the excitation nodal momentum 
k^ = k^^ — if is approximately given by k^ « kp [qarc/^Bsi\- -^ comparative analysis of the hole Fermi momentum 
kp {4>) and nodal Fermi momentum kp expressions of Eq. ([TT]) and 1-el-sharp-feature line arc hole momentum k ^ {(f)) 
and nodal momentum k^ expressions of Eq. (l88l) reveals that for g^^ = g^^.^ and thus qarc = <7Ssi the sharp 
one-electron spectral features line centered at —7? is the Fermi line. In turn, for g^^ = g^ and thus qarc — Qec where 
(f) — 7r/4 the arc of such a line considered here becomes a single discrete momentum value. It corresponds to a nodal 
momentum that we denote by /cj^ and reads k^ « ~[^f /^Ssil ^ec ■ Here (fj^ is the auxiliary nodal si — sc-line 
momentum of Eq. fTS)) whose absolute value is given in Eq. ([75)1 . 

The studies of Refi^ check the validity of the strongly coupled virtual-electron pairing mechanism introduced in 
this paper in experiments on hole-doped cuprates. As discussed in that reference, the 1-el-sharp-feature line arcs, 
whose energy provides a direct signature of the strongly coupled virtual-electron pairs, are observed in the experiments 
on LSCO of Refi^. The studies of Rcf.^"^ confirm that the boundary that separates the broad one-electron removal 
spectral features from the 1-el-sharp-feature line arcs is indeed the theoretical line Ey-ei{<iarc) ~ 2Wec(l — | cos 2(f)\) of 
the virtual-electron energy allowed by the inequality (j86p . It refers to the momenta associated with the minimum and 
maximum Fermi angles 7r/4 — (f)arc and 7r/4 -t- (f>arci respectively, of each 1-el-sharp-feature line arc. The theoretical 
expression Ey-^i (Qarc) ~ 2Wec(l — | cos 2(f>\) has exactly the same form and, for the theoretical parameters appropriate 
to LSCO, the same 2Wec magnitude as the corresponding experimental empirical formula given in Eq. (1) of Ref.^i. 
That formula refers indeed to the boundary separating the one-electron sharp features from broad incoherent features 
in that hole-doped cuprate superconductor. Such an overall agreement seems to confirm that the strongly coupled 
virtual-electron pairs are observed in the hole-doped cuprate superconductors, consistently with the pairing mechanism 
proposed here for the VEP quantum liquid. 

E. Zero-temperature superfluid density 

Concerning charge excitations and currents associated with the phases 9cp the zero-momentum c fermion pairs 
participating in phase-coherent virtual-electron pair configurations behave independently of the corresponding si 
fermions of momentum (j. Hence within such phenomena the phases 9cp are associated with the zero-momentum 
c fermion pairs of the phase-coherent-virtual-electron-pair superconducting state macroscopic condensate. In the 
following we call them coherent c fermion pairs. For such a state the energy cost of a phase twist is for small 
< {x — Xc) ^ 1 and zero temperature approximately given by /3cp(0)|V^cpP/8. Here pcpiT) is the superfluid density 
at temperature T of coherent c fermion pairs. For < (x — Xc) ^ 1 the fluctuations of the phases 0j.o and thus of the 
c fermion-pair phases 9cp are strong. Hence the critical temperature Tc is governed by such fluctuations. As a result 
for < (a; — Xc) ^ 1 the transition taking place from the superconducting state to the quantum vortex liquid is for 
the quasi-2D VEP quantum liquid a Berezinskii-Kosterlitz-Thouless like transition^iiS. One finds according to the 
properties of such a transition that under the suppression effects the following relation holds, 

^fcsTc^^^S^lZkl^p (T,). (89) 

IT m* 
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Here UcpiTc) is the density limT->.r<: ncp{T) of coherently paired c fermions for < (Tc — T) ^ 1. (For < {T—Tc) ^ 1 
such a density vanishes.) ncp{Tc)/2 is the corresponding density of the macroscopic condensate coherent c fermion 
pairs. In units of lattice spacing a one, the c fermion mass can be written in terms of the coefhcient Tc of Eq. (|Aip 
of Appendix A as tti* = T? j^Vct. Here fi is the Planck constant. (Often we write such a mass in units of both a and 
Ti one, so that it reads m* == l/2rci.) From the T^ expression given in Eq. (j33]) one finds T^ ~ g[)/S.Q/2kB — >■ as 
(x — Xc) — >■ 0. Indeed, in that limit the suppression coefficient of Eq. ((35)) reads 7^ = 1. For < (x — Xc) <C 1 the use 
of the relations given in Eq. ([55]) then leads to, 



Pcp(O) ~ 50 



4Ao 



(90) 



Here cjq is the amplitude given in Eq. (H^ . So that the relations re ~ 2rs « 2e^"'*""/'^ and a;* = 2rs/7r « 
Tc/tt w p/TrJe^''*""/'^ apply, expression (lUUl) and most expressions given below are derived for the approximate range 
U/At e (woWi). It follows from Eq. ([TO]) that for < (x — Xc) ^ 1 the superconducting-ground-state c band 
momentum area that corresponds to the c fermions contributing to coherent pairing is. 
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(91) 



The hole-momentum domain Q'^ of the c fermions that contribute to coherent pairing refers to a set of c — sc lines 
whose radius q^ belongs to the range q^ e {q'fc^ ^cv)- Here q^p is the radius of the zero-temperature coherent unpaired 
c fermion line whose hole momenta are denoted by qz,- It is the coherent c — sc line of largest radius q^p. For 
< (x — Xc) ^ 1 it is nearly circular and separates in the c momentum band the c fermions participating in coherent 
pairing from those that do not participate. The hole-momentum domain Q^p is bounded by the c Fermi line and 
coherent unpaired c fermion line, respectively. 

For < (a; — Xc) <C 1 the coherent unpaired c fermion line radius q'^ approximately reads, 




(92) 



For U/At £ {uo,ui) this expression is a very good approximation for the range x € {xc,Xci) and a reasonably good 
approximation for the range x € {xci,x^,). 

The energy bandwidth Wcp « [{qcp)"^ — iQFc)'^]/'^''^c of the corresponding superconducting-ground-state sea of c 
fermions contributing to coherent pairing reads. 



Wcp = go 8A0 ; We, 



max Wc-^ 
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(93) 



This confirms the validity of the maximum energy bandwidth expression provided in Eq. (j64p . The parameter jc 
appearing here is given in Eq. (|38l) . The expression of Wcp provided in this equation is a good approximation for 
a range of hole concentrations x S {xc,Xcp) for which Wcp < Wcc- The maximum magnitude W^c of the energy 
bandwidth Wcp is reached at a hole concentration x — Xcp given below. Only for x below Xcp there is formation of 
coherent c fermion pairing upon increasing x. Such a formation occurs while the short-range spin correlations are 
strong enough to supply the energy needed for it. For the hole concentration range x € {xcp, x*) the energy bandwidth 
Wcp is independent of x and reads Wcp = Wee- Furthermore, it vanishes for x > x*. The singular behavior occurring 
at X = X* marks a sharp quantum phase transition to a disordered state without short-range spin order and thus 
without long-range superconducting order. 

For the approximate range U/At G (uo,ui) generalization of the zero-temperature superfluid density expression 
(|M| found above for < (x — Xc) ^ 1 to the hole-concentration range x € (Xc,x*) for which there is long-range 
superconducting order leads to. 
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Here the theta function 0(z) is such that Q{z) = 1 for z > and 8(z) = for z < 0. The superfluid density Pcp 
vanishes both for x < Xc and x > x*. Hence it has a singular behavior at x = x*. It marks the above mentioned 
sharp quantum phase transition to a disordered state without short-range spin order and thus without long-range 
superconducting order for x > x.^. 
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Why does the superfluid density ([M)) not depend of the suppression coefficient 7ci? Since there are no suppression 
effects both for < {x — Xc) ^ 1 and < (cc* — cc) <C 1 when 7d = 1 and such effects are strongest at x = Xop, to 
answer such a question it is useful to consider the ratio, 
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The energy scale 2|f2||T=o is the maximum magnitude of the superconducting virtual-electron pairing energy . The 
dependence on the suppression coefficient 7™™ = (1 — ctd) of the ratio [2|ri|/2|A|]|T=o,a;=2:op justifies physically why 
the zero-temperature superfluid density pcp remains unaltered under the suppression effects and given by Eqs. ([Ml) . 
If the ratio [2|r2|/2|A|]|r=o,x=Sop was independent of 7™*", the superfluid density would be suppressed. In that case 
one would expect that pcp oc 7™*" at cc = Xop- However, while the energy available to the short-range correlations 
2|A||t=o to supply pair formation remains unaltered, the energy cost 2|r2||T=o of that pair formation decreases by a 
factor given exactly by the suppression coefficient 7™™. This effect cancels the decreasing of pcp, whose magnitude 
remains independent of 7™*". Such an analysis is straightforwardly generalized to the whole range x € (xc,Xa.). 

It follows from expression (j94p that the corresponding superconducting-ground-state energy bandwidth of the sea 
of c fermions whose pairs contribute to the superfluid density is for hole concentrations in the range x G (xcX*), 
vanishing spin density m = 0, and zero temperature given by, 



W,p = 



90 Q{xcp ~ x) + -— 9{x - Xcp) 



8A0 . (96) 
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The in-plane penetration depth Xab can be expressed in terms of the superfluid density of Eq. (|94|) as follows. 
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Aab = -^^T^^/tT-^— ; T2- = —l— Pcp- (97) 



Here (i|| is the plane separation of the system described by the quasi-2D Hamiltonian ([5]), — e denotes the electronic 
charge, and po is the vacuum permeability. As confirmed in Ref.^^, for the parameter values appropriate to the 
five representative hole-doped cuprate superconductors the penetration depth Xab is much larger than the coherence 
length ^ of Eq. (|46l) . 

Finally, the x dependence of the zero-temperature superfluid density pcp of Eq. (j94p can for approximately U/At € 
{uo,ui) be obtained from solution of the following rate differential equation of supcrfluid-density increase upon 
increasing the hole concentration x, under suitable physical boundary conditions, 

dpcp _ ^ec „, _ , I _ I _ Wec_ _ I _ I _ 

■ ^\Xcp Xj , Pcp\xcp<x<x^ — Pcp\x—Xt.p — .9 ' Pcp\x—Xf. — Pcp\x—x^ — U . 



dx 2h?r rcp\Xap<^X<^X, rcp|X-Xcp ^^2 

F. General superconducting virtual-electron pairing energy 

Here we generalize the expression of the superconducting virtual-electron pairing energy per electron |ri5i(q)| given 
in Eq. (j6ip for momenta g « QbsI ^^ '-'^ near the si boundary line to momenta q — q^rc belonging to any other phase- 
coherent virtual-electron pair configuration si — sc line arc, as defined below. The whole set of such lines generates 
the momentum range qG Q^^ of the phase-coherent virtual-electron pair configurations. The general function |risi(g)| 
given here is that involved in the si band energy dispersion esi{q) expression of Eq. ()55|) . That dispersion appears in 
the first-order terms of the general energy functional provided in Eq. (IM|) . 

The general coherent virtual-electron pairing energy per electron |ilsi(Q)| refers to pairs of c fermions whose hole 
momenta belong to any c— sc line between the c Fermi line and the coherent unpaired c fermion line. The radius q^ of 
such a c — sc line belongs to the range g'* G {qpc^ Qcp)- Only for it is the coherent virtual-electron pairing energy per 
electron introduced here finite. The energy bandwidth Wcp of such c fermions provided in Eq. (1961) plays the same 
role for the coherent unpaired c fermion line as the energy bandwidth Wee given in Eq. (IM|) for the ec-pairing line. 
As found below, for x € {xcp, x^) these two lines are the same line. 

For the range of hole concentrations x € (xc-,Xa.) the shape of the coherent unpaired c fermion line is at zero 
temperature defined by the following relation, 

q^p G coherent unpaired c hue <^=> ^cpiqcp) = . (99) 

Here the c fermion energy dispersion ecp{q'^), 

tcp{q^) = Wep + ec(g '') , < t,p{q'') < W,p ; t,p{q!^j,) = ; e^piq^^}) = W^p , (100) 
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has the vanishing energy level at the coherent unpaired c fermion line. For C//4t e {uq, ui) and the hole-concentration 
range x € {xc,x^,) the coherent unpaired c fermion line radius is approximately given by Eq. (|92l) . 

The corresponding si band momentum domain Q^l^ of the phase-coherent virtual electron pair configurations refers 
to the set of coherent si — sc line arcs associated with coherent c — sc lines contained in the c band hole- momentum 
domain Q'^p- For instance, for the si — sc line arc corresponding to the Fermi angle range 4> S (0, 7r/2) whose si band 
quadrant has auxiliary momentum angle range (jisi G (tt, 37r/2) and components qoxi < and qox2 < it is limited 
by the si boundary line and the coherent si — sc line arc with minimum nodal momentum absolute value q^ . The 
latter is given by expression (pO| at q'^ — q^^. This leads to, 

for X G {xc,Xcp) and q^ = g^ for x G {xcp,x^,) where g^ is provided in Eq. ([75]). The si boundary momenta qg^i 
and q^^i appearing here are given in Eq. (IE1[) of Appendix E. 

For the hole concentration range x € {xc,Xcp) the coherent si — sc line arc angle width 2(j)arc belongs to the range 
20arc e (7r/2 - 20J^P^,7r/2 + 2(/)^Pj. Here cj)"/^^ is the angle (t)arc of Eq. ^ at 9^^ = g^ and thus g'' = g^p. It reads, 

Crc = o arcsin -^^^ 2^ = - arcsin I 1 e (0, 7r/4) ; a; S (xc, Xcp) , 

^ \1Bs1~ lecj ^ \ X* J 

— 0; X e {xcjj,x^) . (102) 

Its limiting values are ip'^^^ = 7r/4 at x — Xc and (/)^^^ = for a; e (xcp, x*). 

The c fermions of hole momenta q^ and — g'' at or near the c Fermi line participate in phase-coherent virtual- 
electron pair configurations of momentum q « (Ibsi ^^ o^ near the si boundary line. The strong effective coupling 
of such c fermions results from residual interactions with the corresponding si fermions of momentum q « (IbsI ^^ 
or near the si boundary line. In turn, in the opposite limit of c fermion hole momenta g'' and —q'^ belonging to 
the coherent unpaired c fermion line such a residual interactions are with si fermions of momentum g belonging to 



the coherent si — sc line arc whose nodal momentum has minimum absolute value q^ . For the hole concentration 

line arcs and its nodal momentum has absolute value g^ given in Eq. (jlOip . In turn, for x G {xcp,x*) the angular 
width 2(parc vanishes and q^ equals the absolute- value nodal momentum g^ given in Eq. ()78p . For the latter x range 
the coherent si — sc line arc of minimum nodal momentum absolute value g^ then reduces to the discrete nodal 



7^. 
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range x € (xc^Xcp) that line has the minimum angular range (j) G (7r/4 — ^J^f^^, 7r/4 -I- (j)'^^^,) of the coherent si — sc 
,^„1 ™^™„„+,,™ u^„ „K„^i,,+„ ,,„!,,„ yN gjven in Eq. (jlOip . In turn, for x G {xcp, 

dcp 

of minimum nodal momentum absolute value i,^p 
momentum g'= qeJ^ ■ Hence for x G {xcp, x^.) the coherent unpaired c fermion line is the ec-pairing c fermion line. For 
intermediate g'' and — g''* hole momenta at or near the coherent c— sc lines between the c Fermi line and the coherent 
unpaired c fermion line the strong effective coupling associated with coherent pairing is due to residual interactions 
of the c fermions under consideration with si fermions of momentum at or near the corresponding coherent si — sc 
line arcs whose nodal momentum absolute value belongs to the range g^^ G {qcp:QBsi)- 

Based on the results of Sections IV-C and IV-E, we are now ready to introduce the superconducting pairing energy 
spectrum |J7si(g)| contributing to the si band energy dispersion esi(g) expression (|55p . For a virtual-electron pair 
of momentum g ~ q^^.^ belonging to a given coherent si — sc line arc the general superconducting virtual-electron 
pairing energy per electron reads. 



1^.1(^0)1 - e {e^piq'^) ~ V!l^{qtS) & {Wcp - e,p{q'')) 7. ffoV2 ^f (gt.J ^-^(g") +;-^( g') 
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(103) 



Here q^ is the radius of the corresponding c— sc line. It is uniquely related to the nodal momentum absolute value 
g^j, of the coherent si — sc line arc under consideration by the function g^^ = QarAl^) defined in Eq. (|80|. Moreover, 
since coherent pairing requires strong effective coupling, the energy range of Eq. ([S7)) together with the expression 
€cp{q^) = Wcp + £c{q^) of Eq. (|100l) imply that the energy scale V^J^ {q) appearing in this pairing energy expression 
reads, 

V^"iqL) = ; W.p< W.. (1 l^^i(g'-)l 
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The pairing energy given in Eq. (|103p is expressed in terms of the energy ecp{<i^) + ^cp{—q^) of the c fermion pairs 
of hole momentum at or near the coherent c — sc line of radius q''. Those interact with the si fermion of momentum 
q = q^^^ at or near the corresponding coherent si — sc line arc. This occurs within the phase-coherent virtual-electron 
pair configuration of the same momentum. The use of the function g^^, — qarAq^) of Eq. (|80| leads to the following 
equivalent yet simpler expression for the general pairing energy per electron defined by Eqs. (jl03p and (jl04p . 

I^si(ejl =7^50^2 ^i(ej f 'r''x l «7.ff f 'r''x l Ao|sin20|. (105) 



qssi Vcp J \ qBsi qcp 



Here, 



Kt(q) = |V,HA.i(g)|| = ^ G.i(g) , g = g^c ; G.i(g;t) « I sin20| , (106) 

where the function Gsi{q) is defined in Refji. The value Gsi{qarc) ~ |sin2(/)| is a good approximation for si band 
momenta belonging to coherent si — sc line arcs. 

The energy |r2si((f^rc)l has its maximum magnitude for momenta q^^^ « 'fJsn when it involves c fermions of hole 
momenta at or near the c Fermi line such that ecp{iq^) = Wcp- The strong effective coupling of such fermions results 
from interactions with si fermions of momenta at or near the si boundary line. It vanishes in the opposite limit 
referring to c fermions of hole momenta at or near the coherent unpaired c fermion line for which ecp{iq ) — 0. 

Both the superconducting pairing energy expressions (jl03p and (jl05p have physical significance only for si band 
momenta q ~ q^^^ at or near a coherent si — sc line arc. Interestingly, the energy |r2si((f^^^)| has a much simpler 
expression (jlOSp in terms of the momenta of the virtual-electron pair coherent si — sc line arc than as given in Eq. 
(I103p . Indeed, the coherent si — sc line arcs have been constructed to inherently the inequalities (|104p being obeyed. 

On varying the virtual-electron pair coherent si — sc line arc nodal momentum absolute value within the range 
Qarc S (9cpi qssi) ^^'^ ^^us that of the c fermions hole momenta within q'^ € (^fc 9cp) the pairing energy 2|r2si(gjf^j,)| 
varies within the corresponding range 2|f2si(<f^^^)| S (0, 2|r2|| sin20|). It has the following Hmiting behaviors, 

max2|r!,i(<7;tJ| =2|f],i(g|i^)||T=o = 2|r!||T=o; 2\n,,iq^p)\ ^ ; 2\n,,iqi^,)\ ^ 2\n\\sm2^\ . (107) 

V. CONCLUDING REMARKS 

In conventional superconductivity of isotropic 3D many-electron problems the objects that pair are the quasipar- 
ticles of the Fermi-liquid theory—. The low-energy eigenstates are described by occupancy configurations of such 
quasiparticles so that their momenta are good quantum numbers and their interactions are residual. That property 
simplifies enormously the description of the many-electron physics^S. In turn, in the VEP quantum liquid the charge 
c fermions and spin- neutral two-spinon si fermions are for x € (ccc, x*) the constituents of the phase-coherent virtual- 
electron pair configurations. The momentum occupancies of such quantum objects generate the energy eigenstates of 
the Hubbard model on the square lattice in the one- and two-electron subspace^. For the VEP quantum liquid studied 
in this paper such c and si band momenta are close to good quantum numbers. The c — si fermion interactions play 
a central role in the pairing mechanism of that quantum liquid, consistently with the evidence that unconventional 
superconductivity is in different classes of systems mediated by magnetic fluctuation o^^'^^ . 

Our studies consider the fluctuations of the phases associated with the long-range antiferromagnetic, short-range 
spin, and long-range superconducting orders occurring at zero temperature for x G (0,xo), x S (xo,x«), and x S 
(xcX*), respectively. For x G {xo,Xc) the short-range spin order coexists with Anderson insulating behavior brought 
about by intrinsic disorder. For hole concentrations in the range x G (xc,a;*) the short-range spin order and long- 
range superconducting order coexist. Fortunately, for that x range the effects of intrinsic disorder or superfluid-density 
anisotropy are very weak. This justifies the success of our oversimplified description of such effects in terms of a single 
suppression coefRcient -y™*" ^ (l—ad). The results of this paper focus on the physics associated with the corresponding 
hole concentration range x € (xc,x*). (The effects of intrinsic disorder are much stronger for x € {0,Xc), as shortly 
discussed in Appendix B.) 

The exotic d-wave superconducting phase of the VEP quantum liquid involves virtual-electron pairs whose phases 
are coherent. Those are constituted by one zero-momentum c fermion pair and one spin-singlet two-spinon si fermion 
of momentum q G Qcl,- Hence a virtual electron pair involves the two charges — e of its c fermions and the spin-singlet 
configuration of the two spin-1/2 spinous of its composite si fermion. The phases of the virtual-electron pairs read 
— 6o + 9i. Here do are overall centre-of-mass phases and the phases 9i are related to the internal pairing degrees of 
freedom. The corresponding macroscopic condensate of zero-momentum c fermion pairs is associated with the phase 
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coherence occurring for x € (xc, a;*). The fluctuations of the phases 6q and 6i become large for < (2: — Xc) ^ 1 and 
< (x* — x) ^ 1, respectively. For the intermediate interaction range U/At £ {uo,ui) the critical hole concentrations 
are given by Xc ~ 0.05 and x* G (0.23, 0.28). For it such fluctuations are controlled by the weak 3D uniaxial anisotropy 
effects and electronic correlations through the mass ratios 1/e^ = M/m* and re — m'^/m* « 2rs « 2e~''*""''^, 
respectively. The direct relation to the problem investigated in Ref<^ has simplified the study of the fluctuations of 
the phases 9q and 61. 

The virtual-electron pairing energy is supplied to the zero-momentum c fermion pair by the short-range spin 
correlations associated with the spin-singlet spinon pairing. This occurs through the c - si fermion interactions 
within each virtual-electron pair. The interrelated spinon and c fermion pairings involved in the exotic coherent 
virtual-electron pairing are associated with the coexisting two-gap scenario consistent with the unusual properties of 
the representative hole doped cuprate o^'^^i'^^'^^ : A superconducting energy scale 2|51| and pseudogap 2|A|, over the 
whole dome x G (xcX*), as illustrated in Fig. 1 of Refji^. The energy parameter 2|A| is the maximum magnitude 
reached at (p = 0,tt/2 of the spinon pairing energy 2\Asi{qggi)\ « 2|A|| cos20| appearing in Eq. (|62|) . The energy 
scale 2|f2| is the maximum magnitude reached at (j) — 7r/4 of the superconducting virtual-electron pairing energy 

2\Qsi{qBsi)\ ~ 2\n\\sm2(j)\ of Eq. (|107p . The latter pairing energy is associated with the effective coupling of 
the corresponding zero-momentum coherent c fermion pairs. Such zero-momentum c fermion pairs couple to charge 
probes independently of the si fermions, which couple to spin probes. It follows that concerning charge excitations 
the zero-momentum coherent c fermion pairs behave as independent objects relative to the virtual-electron pairs. In 
turn, the virtual-electron pairs exist in virtual intermediate states generated by one-electron and spin excitations, 
which break such pairs. 

The pseudogap state occurs both for hole concentrations x e (xo,Xc) at temperatures below the pseudogap tem- 
perature T* and for a temperature-dependent hole concentration range centered at a; = Xop for finite temperatures in 
the range T S {Tc,T*). For such a state the si fermion spinon pairing energy remains finite and there are c fermion 
pairing correlations, yet there is no phase coherence. Normal ground states may be reached by applying a magnetic 
field H aligned perpendicular to the planes. Our results are inconclusive on whether for H ~ the ground state 
of the e^ = m*^/M = Hubbard model on the square lattice is superconducting. They seem to indicate that some 
small 3D uniaxial anisotropy is needed for the emergence of superconductivity. Indeed, at constant U /At values the 
hole-concentration width {x■^, — Xc) of the superconducting dome decreases upon decreasing the small 3D uniaxial 
anisotropy parameter e^ = m*elM . 

Concerning previous related studies on the large-?/ Hubbard model and t — J model and Heisenberg model on 
a square lattice involving for instance Jordan- Wigner transformations^^ or the slave particle formalisro^i^Sr— , here 
the single occupancy constraint is naturally implemented for all U /At finite values. Indeed, the spin-1/2 spinous 
refer to the rotated electrons of the singly occupied sites of the energy eigenstates configuration a^'^^ . In the one- 
and two-electron subspace where the VEP quantum liquid is defined only the charge c fermions and spin-neutral 
two-spinon si fermions play an active role. The main difference relative to the above related schemes is that for 
them the spinless fermions arise from individual spin-1/2 spins or spinous. In contrast, according to the studies of 
Re£i the si fermions of the representation used in the studies of this paper emerge from an extended Jordan- Wigner 
transformation performed on spin-neutral two-spinon composite si bond particles. A property of physical importance 
for the unusual scattering properties of the VEP quantum liquid investigated in Ref.^"^ is that the si band is full for 
m = ground states and contains a single hole for one-electron excitations. 

When expressed in terms of the rotated-electron operators, the VEP quantum liquid microscopic Hamiltonian 
has basic similarities to that considered in Refi^ in terms of electron operators. The advantage of our scheme is 
that it applies to intermediate U/At values. Indeed, the results of Refs.^^''^'^ strongly suggest that the physics of 
several classes of hole-doped cuprates with critical hole concentrations Xc ~ 0.05 and x* ~ 0.27S"— corresponds to 
U/At « u* — 1.525 G (uoi^i)- The effective superconductivity approach for such cuprates introduced in this paper 
is consistent with their physics being closely related to the doping of a Mott insulator'"' . There are previous studies 
where spin-charge interactions also lead to superconducting effective pair couplingi^i^"— . Our VEP quantum liquid 
acts on the one- and two-electron subspace^. For it rotated-electron double occupancy is exactly zero. This implies 
that the mechanisms behind its long-range superconducting order are different from those of the 50(5) theory^S. 

Finally, in RefsJ^i^^ evidence is provided that the use of the theoretical VEP quantum liquid scheme introduced 
in this paper leads to quantitative agreement with experiments concerning the universal properties of the cuprate 
superconductors 
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Appendix A: Basic information on the c and si fermion description 



In this Appendix some basic information on the c and si fermion description used in the studies of RefSfii^^ is 
provided. The critical hole concentration x* expression given in Eq. (|47p refers to the interaction range U /At > mq ~ 
1.302 and involves the spin ratio r^. That ratio can be defined for the whole range of U /At values and increases 
smoothly upon increasing U/At. The charge and spin ratios are for the approximate range U/At € (uq^ui) given by 
Tc « 2rs « 2e~'**"°/'^. Here ui « 1.600. For U/At > such ratios have the following limiting behaviors, 
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For the range U/At S {uq,ui) the critical hole concentration x* is given by x* « 2rs/7r. For U/At > uq its limiting 
magnitudes are, 

X, « Ig-ituo/U ^ U/Ate (wo, Ml), 
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An important energy scale is the the a; == and m = Mott-Hubbard gap 2/i'^ given bj^, 
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Another energy scale that plays a key role in our studies is that of Eq. ([6]), lim^j^o |A| = Aq. It has the following 
approximate limiting behaviors^, 

Ao « 16t e-^\/4*7^ , U/At < 1 , 
= max {Aq} « t/-K , U/At = mq j 

« e(i-4*"°/^)[i/7r][l-([//4i-Mo)e"''^^^^'"^"°'], f//4i € (uq, Ui) , 

« 7r[4i]V[/, C//4t> 19, (A4) 

where uq « 1.3, w* = 1.525, and ui ~ 1.6. The U/At dependence reported here for approximately U/At e (uo,ui) 
corresponds to an interpolation function used to connect the following limiting behaviors valid for < {U/At — 
uo)/{ui - Uq) < 1 and U/At sa u* = 1.525i, 
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Let us consider a Brillouin zone centered at the momentum — tt. 
Fermi momentum kp as follows, 



The hole Fermi momentum kp is related to the 



kp — kp + n — kp {(f>) 6*0 ; (j) = arctan —f^ 
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Here the expression of the Fermi angle 4>, which defines the direction of the hole Fermi momentum kp, also holds for 
(j) S (0,27r). The angle (/)an appearing in Eq. (|A6[) is small for the hole concentration range x G {xh,x^.) for which 
the Fermi line is particle like. The precise magnitude of the hole concentration Xh > Xc2 remains unknown. The 
studies of Refiii^ consider that it belongs to the range Xh S {xc2 ,x^,). The hole Fermi momentum fc^ can be expressed 
in terms of the c Fermi line hole momenta and si boundary momenta as given in Eq. (jlOp . For x G (xci,Xc2) the 
corresponding angles (/)^ and ^f^ are provided in Eq. pT|) . 
For one-electron excited states, the si band momentum, 



qssi = [Ai^ 



"^ ' VJsi =qBsi{(t>)e^^ 



(A7) 



is called auxiliary momentum of the si boundary-line momentum q^^i- Both for x < Xci and x > Xc2 the form of the 
matrix [^^j^]"^ is not known. For approximately x G {xci,Xc2) it reads [^fj^]"^ « Ap where the F rotation matrix 



Ap is given in Eq. p5|) . 

Among the si boundary-line momenta qg^i , it is useful to consider those whose corresponding auxiliary si boundary- 
line momenta qssi point in the nodal and anti-nodal directions. The nodal and anti-nodal si boundary-line momenta 



-dJV 



and qg^^ , respectively, are defined as those whose corresponding auxiliary momenta q^^i and q^ 



instance for the quadrant such that go^i ^ and qQx2 < the following Cartesian components. 



'-^fi have for 



9f.i = 



qli/V2 
qLi/V2 



-an 



AN 
^Bsl 





AN 

qssi 



(A8) 



Here q^^i and g^^ are the absolute values of both the auxiliary momenta q^gi and q^si ^^^^ corresponding momenta 
(fjj^ and gj^^, respectively. 

For m = ground states and their two-electron excited states the si boundary-line momentum reads q^gi — qBsi — 
qBsi(,4') e^-i-Tr for X e (xc, a;*). Thus it equals the auxiliary momentum. For the sake of generality, we often use the 
notation q^^i for the si boundary-line momenta of the si fermion occupancy configurations of these states. For such 
states the hole c Fermi momentum is independent of the doublicity d introduced in Ref.^. It is given by. 



qk 



qp, 



+ TT ^ qp^{(l)) 64,^ ; (t>c = (f>- 



(A9) 



The c and si energy dispersions appearing in the general energy functional of Eq. (I54p depend on the Cartesian 
components of the c band hole momentum and si band momentum through the elementary functions edq) and esi(<z), 
respectively. Those are known in some limitsi. The c fermion energy dispersion edq^) is for U /At > 0, tti = 0, and 
X G (xcX*) given by. 



=(9") = ^liq'l 



El' 

1=1.2 



Ml) 



edq'p' 



)]• 



(AlO) 



For X G {xc, x■^,) and U/At > uq the c Fermi line is approximately circular and the c fermion energy dispersion (|A10p 
and the chemical potential /i are given bj*^, 



^^(^-..)^_M!f_M^. \,^(^q'^)\^W^%=.,; ^'^^' M«M0 + W^.^ 



(All) 



Provided that ^0 is replaced by the energy scale /i° given in Eq. (|Bip of Appendix B, the /i expression given here is 
a good approximation for the VEP quantum liquid at a; G {Xc,x^,) and U/At G (mo,wi). 

The si fermion energy dispersion e'^i{q) and pairing energy per spinon |Asi(g)| appearing in Eq. (|55|) are given by. 



4('?)-4"('fo); |A,i(q)| = |AL(go)| 



(A12) 



Here qo obeys Eq. dH), e"J{qo) = J2z=i,2[^siiqox,) - esi{qBsixi)], |A!|;^(qo)| = gi ^0 F'J^{cfo) , i^Ji((fo) = [\esi{qoxi) 
esi{qox2)\]/Wsi, and the energy bandwidth Wsi decreases for increasing values of U/At. 
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The elementary function esi{q) is such that Asi{qg^) — 0, |esi(g^^) — esi(0)| = Wsi, and Fl.-^{qBsi) ~ |cos20|. 
Hence at T = the maximum magnitude 2|A| of the spinon pairing energy 2|A5i(q)| is reached at q = qj^fi^ , 

2|A||T=o=5i2Ao = 2|A,i(q'|,^i^)|; x G (x„ a;,) , m = . (A13) 

The shapes of the c Fermi hne and si boundary Hue are fully determined by the form of the auxiliary energy 
dispersions e^{q''') and €^i{q) given in Eqs. (lAlOp and (|A12|) . respectively, as follows, 

qpc e hole c Fermi line <^=^ e°{qFc) = , 
(fssi e si boundary line <=> e°^(qj^^)=0. (A14) 

The c and si fermion group velocities derived from the energy dispersions of Eqs. ([55)) . (jAlOp . (IA11[) . and (IA12[) 
read, 

Vc{q^)^y,-.e,{q'^); Kf(q):=-V,HA,i(<r)|; F,i(qO - V,^e.i(g) ; Vjiiq) ^ ^ .^ ^lio) ■ (A15) 

We call their unit vectors e'0^j(^, e^o ^-j, e^Aj--), and e^^i^^h^, respectively. 
Finally, the c fermion elementary charge current is given by^, 

Jc(g'') = -eac/^^c(g''); JcC^fJ = -e%e^.+.; a^; = ^ « ^ , f//4i > uo « 1.302 . (A16) 

Here Vc{q^) is the c fermion velocity of Eq. (jAlsp . 771^ a renormalized transport mass, and 4>c — (f>- 

Appendix B: The hole-trapping effects 

In this Appendix the strong effects of intrinsic disorder on the square-lattice quantum-liquid physics of Ref.— in 
the hole concentration range x S {0,Xc) are briefly discussed. In contrast to the range x S {Xc,Xi,) referring to the 
problem studied in this paper, the effects reported here change such a physics qualitatively. 

As discussed in the following, intrinsic disorder is strongest in the range x G (0,a;o). It has qualitatively different 
effects in that range and in the range x G {xo,Xc), respectively. Here xq = x^/[l + (27ra;*)^ (t/Ao)] is the hole 
concentration of Eq. ((8)) at which the Fermi-energy anisotropy coefhcient 770 = jAI/W^'' given in that equation reads 
770 = 1 for the square-lattice quantum liquid. It is given by xq « 0.013 and Xq « 0.024 for x^ = 0.27, t — 295 meV, 
and the magnitudes Aq — 42meV and Ao = 84meV found in Refi^^ for LSCO and the remaining representative 
systems, respectively. For the square-lattice quantum liquid of Ref;^ there occurs a quantum phase transition from 
a Mott-Hubbard insulator with long-range antifcrromagnctic order at the hole concentration a; = to a short-range 
incommensurate-spiral spin ordered state for < a; ^ 1^^ . An energy scale that could be used to characterize such 
a transition is the c fermion unfilled sea bandwidth W^ of Eq. (jAlip of Appendix A. It is associated with the hole 
kinetic energy. For the square-lattice quantum liquid it vanishes at a; = and is finite for finite hole concentrations. 

Within the effective and simplified description of the problem considered here, for hole concentrations in the range 
X € (0, Xc) the holes are trapped in the vicinity of randomly distributed impurities, so that their kinetic energy vanishes 
and the c fermion hole energy bandwidth W^^ is not well defined. Nevertheless, the critical hole concentration xq of the 
Neel-state - spin-glass-state sharp quantum-phase transition is that at which the ratio jAj/W^'' of the corresponding 
square-lattice quantum liquid is given by one. Such a ratio is for the Fermi-energy anisotropy coefficient 770 — |A|/Wj, 
of Eq. ([5]) introduced in Ref.^. In turn, the absolute value |esi((75| of the si fermion energy dispersion of Eq. (JA11[) of 
Appendix A changes from |esi(<7Bii)| ~ '-' ^^r si boundary-line auxiliary momenta pointing in the nodal directions to 
|esi(gj^^) — |A| for such momenta pointing in the anti-nodal directions. Therefore, the energy parameter |A| gives 
a measure of the si boundary-line anisotropy. 

Consistently with |A| providing a measure of the si boundary-line anisotropy, for the square-lattice quantum liquid 
at hole concentrations below xq for which |A| > W^} the corresponding anisotropy effects dominate over those of 
the isotropic c Fermi line. It follows from Eq. (fT0|) that the Fermi line involves contributions from both that line 
and the si boundary line. Hence for the square-lattice quantum liquid the Fermi line is strongly anisotropic for hole 
concentrations below a:o- This weakens the hole kinetic-energy effects associated with the energy scale W^}. The 
hole-trapping effects profit from such weakening of the hole kinetic-energy effects: They give rise to the extension of 
the long-range antiferromagnetic order from a; = to a; G (0, a;o) for the VEP quantum liquid. As a result, the critical 
concentrations xq and Xc are those at which the energy scales |A| and W^ emerge for the latter quantum problem, 
being ill defined for hole concentrations below xq and below Xc, respectively. In contrast, for the square-lattice 
quantum liquid they are well defined for the whole range a; G (0,a;*). 
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TABLE I: Magnitudes of the Fermi-velocity anisotropy coefficient tja and Fermi-energy anisotropy coefficient r/o of Eq. ^ for 
several hole concentrations x, U/At « 1.525, and thus Ao ~ 0.285 1. It is considered that the Fermi-velocity anisotropy is small 
when rjA < 2x0 ~ 0.05. Consistently, the crossover hole concentration Xci = 1/8 approximately is that at which tja ~ 2x0 for 
U/4:t £ (uo, 111). The magnitudes given here refer both to the square-lattice quantum liquid and the four representative systems 
other than LSCO. 



The hole-trapping effects lead to two related phenomena: (i) the extension of the long-range antiferromagnetic 
order of the Mott-Hubbard insulator from a; = to the range x G (0, a:o); (ii) The coexistence for x S {xq, Xc) of the 
short-range spin order emerging at x = a;o with Anderson insulating behavior. The latter behavior is brought about by 
intrinsic disorder. For hole concentrations in the range x € {xq, Xc) such a disorder and associate hole-trapping effects 
are not strong enough to remove the short-range spin order. Finally, the only effect of the hole-trapping effects on the 
long-range superconducting order of the quantum liquid with weak 3D uniaxial anisotropy is shifting its emergence 
by Sxc — Xq. This changes the magnitude of the critical hole concentration Xc from Xc = Gi for that quantum liquid 
to Xc = Gi + Xq for the VEP quantum liquid. In turn, the magnitude of the hole concentration cc* remains unaltered. 

For the square-lattice quantum liquid the hole concentrations xq and Xd — 1/8 of Eq. ([8]) do not mark sharp 
quantum phase transitions. For intermediate values of U/At they refer instead to crossovers between hole-concentration 
ranges where the interplay of Fermi-line anisotropy and electronic correlations leads to a different physics. Here that 
remains true for Xci = 1/8. In turn, the hole-tapping effects render xq a critical hole concentration. It marks the 
sharp quantum phase transition between the Mott-Hubbard insulator with long-range antiferromagnetic order and 
the Anderson insulator with short-range spiral incommensurate spin order. 

The dependence (and independence) on the Fermi angle (j) of the si fermion velocity V^^ of Eq. © (and c 
fermion velocity Vpc also given in that equation) implies an anisotropic (and isotropic) character for the si boundary 
line (and c Fermi line.) The Fermi angle (j) given in Eq. (jA6|) of Appendix A defines the direction of the hole 
Fermi momentum kp also provided in that equation. According to the analysis of Ref.i, the Fermi line anisotropy 
is measured both by the Fermi-velocity anisotropy coefficient tja = maxjrAJ = maxly^^-^j/VFc and Fermi-energy 
anisotropy coefficient 770 = max{5Ep}/Wc. For the square-lattice and VEP quantum liquids, those are well defined 
for the ranges x G (0, x*) and x G (XcX^,), respectively. The velocity Vg^i and anisotropic part of the Fermi energy 
SEp of Eq. (1201) achieve their maximum magnitudes for si band auxiliary momenta pointing in the nodal and anti- 
nodal directions, respectively. Hence the coefficients 77A and ijo contain complementary yet different information. For 
approximately U/At S {uq,ui) they are given in Eq. ([5]). The hole concentration xq is defined in that equation as 
that at which tjq — 1 for the square-lattice quantum liquid. For such an intermediate f//4i range it is considered that 
the Fermi-velocity anisotropy is small for hole concentrations larger than Xd — 1/8. 

For U/At « 1.525 the expressions provided in Eq. ([5]) lead to xq sa 0.024 for the square-lattice quantum liquid 
for which Ag/i ~ 0.285 and four representative systems other than LSCO and to xq « 0.013 for LSCO. Indeed, 
according to the analysis of RefJ^, Ag/i ~ 0.142 for the latter system. The magnitudes of the coefficients ?7a and 
r]Q are for U/At « 1.525 and several values of x given in Table |T] for these four representative systems and in Table [H] 
for LSCO. As given in the latter table, the anisotropy coefficients are smaller for LSCO. However, intrinsic disorder 
and randomness are larger for that compound. The Fermi- velocity anisotropy is then considered small at a smaller 
r]/\ = 2xo ~ 0.03 value than that t/a = 2x0 ~ 0.05 of the other four representative systems. Indeed, due to the 
LSCO larger intrinsic disorder and randomness, an anisotropy associated with tja ~ 0.05 has more impact on the 
physics than for the remaining four representative cuprates. However, in both cases the hole concentration at which 
r]A = 2xo sa 0.03 and r]/\ = 2a;o ~ 0.05, respectively, is approximately Xci = f /8. 

As confirmed by the data of the Tables U and [TTl the coefficient 770 decreases upon increasing x slower than r/A . For 
the VEP (and square-lattice) quantum liquid the Fermi-line anisotropy is strongest at x = Xc (and for x — > 0), when 
?7A and 770 are largest (and r/A,Vo ~^ 00). Such an anisotropy vanishes for x — > x* , when r/A, ??o — >■ and the Fermi line 
becomes fully isotropic. Following the x dependence of 77A and 770, it is considered that the Fermi line is anisotropic for 
hole concentrations in the range x G (xc, Xci). For the x range x G (xci, Xc2) it has some Fermi-energy anisotropy, yet 
the Fermi velocity is nearly isotropic. The Fermi line is nearly isotropic for the hole concentration range x £ (Xc2, x*). 
That for x — >■ x* the Fermi line becomes as isotropic as that of an isotropic Fermi liquid is consistent with the physics 
tending to that of a Fermi liquid as that critical hole concentration is approached. 

Both for the x range x € (0, Xc) and hole concentrations obeying the inequality < (x — Xc) ^ 1 the hole-trapping 
effects change the chemical potential /i « /io + W^ of Eq. (|Alip of Appendix A. For the range x S (0, Xc) it is pinned 
and given approximately by // ~ /i*^. In turn, for < (x — Xc) <C 1 it is shifted by —Sfi « —Xc^Tr/m* — —AnrcXct. The 
emergence of its dependence on x is then shifted from the hole concentration x = to x = Xc. The suppression of 
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TABLE II: Magnitudes of the same coefficients as in Table |l] for several hole concentrations x, U/it « 1.525, and Aq « 0.142 i 
for LSCO. Again, it is considered that the Fermi-velocity anisotropy is small when r^A < 2x0. In the present case 2a;o ~ 0.03 
also implies that x ~ Xd = 1/8 for U/At £ {uq, ui). 



the chemical potential shift is for the hole concentration range x G {0,Xc) due to its pinning by the intrinsic-disorder 
impurity potential. For U/4t > uq ~ 1.302 it is approximately given by. 



^ « / + [W^ - Sfi]9{x - Xc) w / + — — ^^^ 9{x -Xc), 0<x <Xc and < (a; - Xc) < 1 , 



M 



M 



+ WJ} ; fl° — fi^ — Sfi; Sfi = AnrcXct , x € (xc, x*) . 



(Bl) 



Here the second chemical potential expression refers to the hole-concentration range x G {xc,Xt.) mostly considered 
in the studies of this paper. For it the only modification of the chemical potential /i « //q + W^ of Eq. (jAlip of 
Appendix A is the replacement of the x independent energy term /io by pfl — [/i" — AnrcXct], the x dependent term 
W^} remaining unaltered. In turn, the first /x expression of Eq. (JB1[) is consistent with for hole concentrations below 
Xc the energy scale W^ being ill defined for the VEP quantum liquid. For the particular case of LSCO, in addition 
to suppressing the magnitude of the energy scale Aq by a factor two, the cation-random effects considered in Ref,— 
cause a further pinning of the chemical potential up to approximately x « 0.14 and lessen the magnitude of W^- 
Hence the formula (|Bip does not apply to that random alloy. 



Appendix C: Hamiltonian terms that control the fluctuations of the c fermion-pair phases 



Here it is shown that the Hamiltonian terms (j29p are approximately equivalent to those given in Eq. 
to reach that goal we start by expressing (|29|) in terms of the phases 9j , 



3T]). In order 



Tjhonds \ ^ A\ 



F Ao/1 /I bi 



J' j'"b-const] 



,sl,dj,g 



(h.c. 



(CI) 



Long-range superconducting order implies the occurrence of coherent c fermion pairing, so that (/I^ ^ fl^^ ^ is finite. 
We thus use a mean-field approximation for which the Hamiltonian ()29|) is replaced by. 



H. 



bonds 






E 

,j"lJ-const] 



Ao(4,,4, .,>4,, 



,sl,d,l,g 



(h.c.). 



(C2) 



According to Eq. (|25p . local c fermion pairs with real-space coordinates r*,/ and rjn equal to those of a given two-site 
bond refer to the same local rotated-electron pair. Hence within our approach for pairs with the same centre of mass 



']/2 



']/2 



' d.l 



and belonging to the same family and thus having the same values for the 



indices d and I (where d = 1,2 is not the doublicity) the following relation holds. 



2h* 



(fl fl ) 

K/i.c/i >l "^' 



E '^'^'^ 

9=0 



2hg\'=l. 



(C3) 



Here the index g was added to specify the two-site link type. Hence the expectation values amplitudes are controlled 
by the coefficient h* of the corresponding spin two-site bond in the si fermion operator defined in Eq. (|26p. Such a 
relation follows from the two c fermions and the corresponding two-site bond stemming from the same spin-singlet 
rotated-electron pair. The expectation values amplitudes decrease upon increasing the distance between the two c 
fermions of a pair. As discussed in Section III-C, the physics behind the relation (jC3|) is that the generator of the spin 
degrees of freedom of the overall occupancy configuration generated by each of the N/2 spin-singlet rotated-electron 
pair operators of the j summation of Eq. (|23|) is a si bond-particle operator. 
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Indeed, the use of the relation (|C3p and si bond-particle operator expression (I^Bl) in the effective Hamiltonian of 
Eq. (|C2p allows after some straightforward algebra to express it in terms of a si bond-particle operator as follows, 

H'efr-i: E ^''^^(4.c4,c)^i.i+(h.c.). (C4) 

The use of the transformation ([3]) and replacement of <t>j^si by the phase (f>^ ^-^ of Eq. ([5^ allows expressing (jC4p in 
terms of the si ferniion operator /I, ^j^ as follows, 

^.o„.. = E E ''"'"^l^fl-cfl,.) 4..1 + (h-c.) . (C5) 

i=iO-'j") ' ' 

Finally, note that (jCSP is the effective Hamiltonian obtained by replacing in Eq. ([31]) the operator /I ^ /I ^ by 
(/I ^ /I ^). This is alike for the effective Hamiltonians (|C1|) and (|C2|) . Since (f29| and (|C1|) are the same Hamiltonian 
terms, one concludes that those are approximately equivalent to (pT|) . 

Appendix D: The energy range of the c - si fermion interactions behind c fermion strong effective coupling 

Here we derive the c fermion energy range of the c - si fermion interactions within a virtual-electron pair configura- 
tion. Our results refer to the interactions that lead to c fermion strong effective coupling. That range is expressed as a 
function of the c fermion hole momenta zLq'^ and si fermion momentum q. We start by considering the square-lattice 
quantum liquid perturbed by weak 3D uniaxial anisotropy without hole-trapping effects. For it there is short-range 
spin order for < .t ^ 1 and the critical hole concentration Xc above which there is long-range superconducting order 
reads Xc ~ Gi rather than Xc ~ Gi-|-xo. Thereafter, we extrapolate our analysis to the range a: € (xc,a;*) of the VEP 
quantum liquid. 

The absolute maximum magnitude of the ratio \il.si{q^r^)\/AQ where |risi(g^^^)| is the pairing energy per electron 
of Eq. (jlOSp gives for 7d = 1 the maximum rate concerning the energy that the short-range spin correlations can 
supply to c fermion strong effective coupling through the c - si fermion interactions within a virtual-electron pair 
configuration. Indeed, in expression (|105|) the suppression coefficient 7^ accounts for the part of that energy used in 
phase-coherent pairing. For 7d = 1 such a ratio reaches its absolute maximum magnitude \ilsi{qg^)\/ Ao — 7c/4 at 
zero temperature, q^^^ — gjj^, and x — Xop- Importantly, as confirmed by Eq. (j64p the ratio Ao/Wec = 7c/4 exactly 
equals that maximum magnitude. 

In the absence of hole trapping effects there is c fermion strong effective coupling for < x < a;* , yet it leads to 
phase-coherent pairing only for x > Xc- The above rate ratio |51si((fJj^)|/Ao — Ao/Wec — 7c/4 plays an important 
role in the c - si fermion interactions. Indeed the c fermion energy range of these interactions contributing to strong 
effective coupling is for < x <C 1 controlled by a function given by, 

|Aec(<r'')| = |^kc(<f'')|; Ceiq"^) = -Wj-^^^ . (Dl) 

Wee ^0 



It is such that, 



|A,,(rf'^)| = |Ae.(-g")|e(0,Ao); |Ae.(<7l;j|=0; |A,e(gl)| - Ao . (D2) 



The momentum tf^'j, appearing here belongs to the zero-temperature ec-pairing line defined in Eq. (j73p . We recall 
that the energy scale Wee of Eq. ((64)) is both the maximum energy bandwidth corresponding to the hole momentum 
domain Q^^ of c fermions with strong effective coupling and the maximum magnitude of the energy bandwidth of the 
superconducting-ground-state sea of c fermions contributing to phase-coherent pairing. Alike the maximum pairing 
energy 2\nsi{qg^)\ is smaller than the energy 2Ao associated with the short-range spin correlations, also the energy 
\Aee{'i'^)\ + |Aec(— g'*) Js Smaller than the c fermion energy |ec(g'')| -I- \ec{—q'^)\- Indeed, for 7^; = 1 the maximum 
magnitude of the ratio 2|0,i(g|i^)|/2Ao exactly equals [\Aee{q'')\ + \Aee{-q'')]/[\ee{q'')\ + \ec{-q'')\]. 

In the present < a; <C 1 limit the energy scale 2Ao of Eq. (jA4p of Appendix A is both the maximum magnitude of 
the si fermion spinon-pairing energy and the maximum magnitude of the c fermion energy scale |Aec('7*'')| + |Aec(— g'*)! 
of Eq. (|D2p . For < a; <C 1 it controls the energy bandwidth of the following ranges, 

2|A,i(gO|e(0,2Ao); |Aec(g")| + |Aee(-<f")| G (0,2Ao) . (D3) 
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Consistently with the ranges (jD3[) and the above analysis, the energy scale 2Ao is for < a; ^ 1 the maximum energy 
of the short-range spin correlations. It is the source of the energy supplied through the c - si fermion interactions 
within each virtual-electron pair configuration, to ensure c fermion strong effective coupling. The occurrence of such a 
strong effective coupling is a necessary condition for virtual-electron pairing phase coherence. This role of the energy 
scale 2Ao is consistent with: i) 2Ao = lima;_j.o 2|A| being the maximum magnitude of the si fermion spinon pairing 
energy associated with the short-range spin correlations; ii) The presence of Aq in the Hamiltonian terms of Eq. (|29p . 
which control the elementary processes associated the c - si fermion interactions. 

Only at zero temperature and for a; — > all phases 6ji line up. Only then gi — |(e*^^-^)| « 1 and the pseudogap 
energy 2|A| = gi 2Ao of Eq. (|36p reaches its maximum magnitude 2Ao. Hence for finite hole concentrations and/or 
temperatures the pseudogap energy 2| A| = gi 2Ao gives the suppressed magnitude of the energy scale 2Ao due to the 
fluctuations of the phases 9ji. That pseudogap is the order parameter of the short-range spin correlations. As a result 
of 2Ao being the maximum energy of the short-range spin correlations, at small hole concentrations < a; <C 1 the 
overall quantity |Aec('f'')| + |Aec(— <f'')|-l-2|Asi((7)| must be equal or smaller than 2 Aq. The ranges provided in Eq. (jD3[) 
correspond to two limiting cases of such an overall inequality: i) The range 2|Asi(g)| G (0, 2Ao) holds for c band hole 
momenta ±q'' at the c Fermi line, so that |Aec(g'')| = \Aec{-q'')\ = 0; ii) The range |Aec((f'')| + |Aec(-g'')| € (0, 2Ao) 
holds for si band momenta (f pointing in the nodal directions, so that 2|Asi((7)| = 0. 

Hence for < x ^ 1 the residual interactions of c fermions with general hole momenta q^ and — g'' and energy 
^c{q^) = ^c{—q^) = —W^ec [|Aec((f'')|/Ao] with a si fermion of general momentum q* and spinon-pairing energy 
2|Asi((j)| contribute to the strong effective pairing coupling of the former two objects provided that, 

< |Aec((f")| + |Aee(-g")| + 2\Asl{q)\ < 2Ao . (D4) 

This inequality is equivalent to, 

M.-)|<M^e.(l-^fP). (D5) 

Generalization of the < a; ^ 1 inequality (jDSp to finite hole concentrations below x* and C//4i g (uo,u^) 
involves the replacement of the energy parameter Aq = liuix^o |A| by the general spin energy scale |A| in the ratio 
[|Asi(q)|/Ao]. For finite hole concentrations below x^. it then becomes the inequality provided in Eq. (|67|) . Due to 
the hole trapping effects discussed in Appendix B, which are strongest for x € (0,a;o), for the VEP quantum liquid 
the latter inequality is valid for x G (xo,a;*) rather than for x e {0,x^). 

Strong effective coupling is required for the occurrence of phase-coherent virtual-electron pairing. However, strong 
effective coupling also occurs in the pseudogap state. That 2\il\ = does not affect though the validity of the 
inequalities provided in Eqs. (|67p and (|D5I) . It defines the c fermion energy range of the c - si fermion interactions 
that lead to strong effective coupling independently on whether it is associated with phase-coherent virtual-electron 
pairing or not. 

Appendix E: The si boundary line nodal and anti-nodal momenta 

The controlled approximations used to derive the following expressions of the si boundary-line absolute-value 
momenta q^gi and q^^i rely on the change of the si boundary line shape and size of the square-lattice quantum 
liquid of Ref.— from a square of edge magnitude V^n and length 1/2477 for x — >■ to a circle of radius ^^(1 — x)27r 
and length (1 — x)2tt'^ for x — >■ 1. For U/At G {uq,ui) the latter shape is a good approximation for x > Xcs where 
Xc3 ~ (xci -I- x») « 0.40 is defined below and Xci — 1/8. The obtained approximate expressions read, 

QbsI ~ 7f*^^~^)' 2;G(Xc,Xci), 



V28 



= ^(1 -Xc3)27r, xe(xi,Xc3) 



AN 
1b si 



\J{1 - x)27r, xe(xc3,l); 

X Kz \Xc-i •^clj 1 

C2 V(l - x)2t^ I xe(Xcl,Xc3), 

y/(T^x)2^, xe(x,3,l), (El) 
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respectively. Here, 



Ci = Ca\x=x^i — -z 



1 — -— tanh(-v/7r) 



'^2 — 



4 V^^s-a^ci 



256 /- CiV77r __ j^N 



The coefficient Ca = 'ZBfi/V2'?Bsi ^^^ then the following approximate limiting behaviors, 



r^, - ggfi ^ l-i^tanh(v^7ra;/a;,i) 




V2 



xe(a;,3,l). (E3) 



As discussed in Ref.-'^, for approximately C//4t > uq and hole concentrations in the range x G (0,a;ci) the shape of 
the si boundary line of the square-lattice quantum liquid is independent of U /At. This justifies the independence of 
U /At for that x range of the q^gi and g^^ expressions given in Eq. (lEip for x e (xc, Xd). Indeed, the results of that 
reference apply to the latter x range. (The behavior q^^i ~ tt — \r2/Kx reported in Refii is valid for 2: <C 1 rather than 
for X < a;ci— ; Note also that for U/At > uq the expression q^fi ~ 7r(l — [l/27r] tanh(-\/87ra;)) of Eq. (|Eip applies; In 
turn, the crossover hole concentration magnitude Xci = 1/8 refers only to the range U/At G (uo,ui), for which that 
expression reads q^f^ ~ 7r(l — [l/27r] tanh(-y7ra;/xci)), as given in Eq. (jEip .) 

For the square-lattice quantum liquid of Ref.- the coefficient Ca = QB^i/V^qEsi ^^ ^Q- l|E3p exactly reads Ca = 1 
in the limit x —^ 0. In the range x € (0,a;ci) for that liquid and x € (xcXd) for the VEP quantum liquid the 
overall decreasing of g^^ and q^^i is very similar, so that Ca — Ci ~ 0.97 at a; = Xci = 1/8. At that x value the 
magnitude of q^^i equals the radius \/(l — Xc3)2tt of the nearly circular si boundary line at a larger hole concentration 
X = Xc3 ~ {Xci + x^). In the range x € (x^jXcs) the q^^i magnitude decreases very little. This justifies why within 
our approximation q^^i ~ [7r/v2](l — Xci) — \/{l — Xc3)2n for x G {xci,Xc3), as given in Eq. (JE1[) . Since for x > Xcs 
it is a reasonable good approximation to consider that the si boundary line is a circle of radius ■\/(l — x)2tt, it is 
assumed that qg^i ~ C2 ^(1 — a;)27r for the range x G {xcijXcs), C2 reaching the value C2 = 1 at a; = Xc3- For 
the latter x range the magnitudes of both q^^i and C2 are decreasing functions of x. While the exact form of the 
decreasing function C2 = C2{x) remains an open problem, that C2 ~ ■y/7r(l — Xci)/2 and C2 = 1 for x — Xd and 
X = Xc3, respectively, is expected to be a good approximation. As given in Eq. (jE2p . here we assume that the 
coefficient C2 decreases linearly as C2 c>c (X3 — x) where 0:3 is provided in that equation. That is the simplest curve 
connecting the above C2 magnitudes at x ~ Xd and x = Xcs, respectively. 
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